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V 

Classical and thermodynamic stability of black holes 



Ricardo Jorge Ferreira Monteiro 
Thesis summary 

We consider the stability of black holes within both classical general relativity and 
the semiclassical thermodynamic description. In particular, we study linearised pertur- 
bations and their contributions to the gravitational partition function. Exploring the 
connection between classical and thermodynamic stability, we find classical instabilities 
and new families of vacuum black holes. 

We start by studying negative modes of black hole partition functions, which repre- 
sent pathologies in the one-loop quantum corrections. In particular, we extend this study 
to charged black holes (Reissner- Nordstrom) , using a method based on gauge-invariant 
perturbations, and to rotating black holes (Kerr-AdS), where a numerical technique is 
employed. In the both cases, we find a negative mode in the region where local thermo- 
dynamic stability fails, as expected. 

We then present the first examples of linearised classical instabilities of vacuum 
asymptotically flat black holes. We analyse numerically perturbations of Myers-Perry 
solutions, both in the single spin and in the equal spins (odd D) cases. For sufficiently 
high rotation, in the so-called ultraspinning regime, new negative modes of the partition 
function may arise whose threshold marks both the onset of a classical instability of the 
black hole (not just of the associated black branes) and the bifurcation to a new family 
of black hole solutions. In the case of singly-spinning solutions, we find the threshold 
stationary modes signalling the instabilities, confirming a conjecture by Emparan and 
Myers. In the case of solutions with equal spins, we are able to find perturbations that 
grow exponentially in time in D = 9 (we believe that this extends to higher odd D). 
Furthermore, the new family of solutions bifurcating at the onset of the instability should 
have a single rotational symmetry, saturating the rigidity theorem. 
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Chapter 1 

Black holes 



Black holes are arguably the most interesting objects in theoretical physics. Understanding 
their dynamics forces us to fit together two widely accepted theories of Nature: general 
relativity (Einstein's classical theory of gravity) and quantum mechanics, a goal that has 
eluded theoretical efforts so far, despite encouraging successes. Black hole thermodynamics 
is at the crossroad between the classical and the quantum pictures. In this thesis, we will 
study the behaviour of black holes as seen both from classical general relativity and from 
their quantum thermodynamic properties. 

This Chapter includes a basic introduction to black holes, a summary of their ther- 
modynamic properties, and a review of black hole solutions in spacctimcs of different 
dimensionality. Only later will we focus on our main subject of research: black hole 
stability. 

1.1 Introduction 

In Newtonian gravity, a massive body can have an escape velocity greater than the ve- 
locity of light. The analogue of such a "black planet" in general relativity is a black 
hole. However, the analogy does not go far since black holes are intrinsically relativistic 
objects highlighting two basic features of Einstein's theory: causal horizons and spacetime 

5 
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singularities. 

According to general relativity, the gravitational force is caused by the curvature of 
spacetime. Spacetime consists of three spatial dimensions and a time dimension put to- 
gether in a geometrical way. (In higher-dimensional gravity, additional spatial dimensions 
are considered.) It is a pseudo-Riemannian manifold {M.,gab) of Lorentzian signature 
whose dynamical metric gab obeys the Einstein field equations. Let Rah be the Ricci 
curvature and R the scalar curvature. The Einstein equations are 



The parameter A, the cosmological constant, denotes the contribution from vacuum energy, 
and, on the right-hand side, Tab is the energy-momentum tensor of all the matter fields. It 
is clear that the dynamics of spacetime is affected by any field content. On the other hand, 
any field is affected by gravity since it lives on a curved spacetime, with massive particles 
moving along timelike geodesies and massless particles moving along null geodesies. 

Solid evidence for general relativity has been provided by the corrections to Newto- 
nian dynamics in the Solar system, and the indirect detection of gravitational waves from 
binary pulsars [l]. 

A crucial distinction from Newtonian gravity is that the "action-at-a-distance" is 
substituted by a built-in causality structure in Einstein's theory. The initial- value formu- 
lation of general relativity splits the 10 Einstein equations into 6 evolution equations and 
4 constraint equations |2|. The latter constrain what is suitable initial data on a Cauchy 
spacelike surface. The former consist of a system of hyperbolic quasilinear equations that 
evolves the initial data in time. The resulting causality structure looks locally like the 
light-cone structure of special relativity. However, since the spacetime is dynamical, there 
is the possibility that the Cauchy evolution of smooth geometry and matter data on a 
spacelike surface may lead to a singularity due to a cathastrophic event, such as the gravi- 
tational collapse of a massive body or a high energy collision. By singularity, we mean that 
the spacetime is geodesically incomplete for timelike or null geodesies, and that geometric 



Rab 



Rgab + ^gab 
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invariants constructed with the metric curvature may diverge. The causal propagation of 
such a spacetime pathology could challenge the physical significance of general relativity, 
in spite of the experimental successes. 

Singularities indeed arise in Einstein's theory as shown by the theorems of Penrose 
and Hawking [3]. Their standard interpretation is that general relativity breaks down for 
curvatures of the order of the Planck scale, e.g. RabcdR"'^^'^ ~ c^/hG where h is Planck's 
constant, giving way to a quantum description of spacetime. Since general relativity is non- 
renormalisable when treated as a quantum field theory of gravitons, a more fundamental 
quantum theory of spacetime is required, and general relativity should be viewed as an 
effective low energy theory. The physical significance of classical general relativity relies on 
a protection mechanism from singularities known as cosmic censorship |4|. This hypothesis 
(in the "weak" formulation) states that no naked singularities exist, i.e. singularities 
arising from realistic matter (except perhaps at the Big Bang) can only form behind a 
surface, called an event horizon, from within which no information may reach observers 
outside that surface. 

A black hole is the region contained inside an event horizon. More precisely, for 
spacetimes with an asymptotic conformal structure, a black hole is the region of space- 
time that does not lie in the causal past of future null infinity, and its boundary in the 
full spacetime A4 is called the future event horizon T-L^ (to distinguish it from the past 
event horizon 'H'^ , the boundary of communication of past null infinity, present in time- 
symmetric solutions but absent for dynamically formed black holes). The singularity 
theorems mentioned above then imply, under certain assumptions for the matter content, 
that there is a singularity in the black hole. 

The formation of black holes through the gravitational collapse of massive objects has 
been studied analytically and numerically (see reviews e.g. in [5] and [6]), and recently the 
formation through high energy collisions has also been addressed numerically |7|. Using 
semi-realistic matter, an event horizon forms in agreement with cosmic censorship. In 
fact, no matter is necessary for the formation of a black hole, as the focusing of incoming 
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gravitational waves may be sufficient [s] . The results support the expectation that a black 
hole will form whenever a given amount of energy is contained in a sufficiently small 
region of space, as proposed by the hoop conjecture |8j. What is truly remarkable is 
that the spacetime (at least in the four-dimensional asymptotically flat case) settles down 
to a unique stationary black hole solution, regardless of the details of the initial matter 
distribution. 

The uniqueness theorems (see [9| for a review) show that four-dimensional vacuum 
black holes which are asymptotically flat, stationary and have a connected horizon are 
described by the Kerr metric [10| . This solution depends on two parameters only: the mass 
M and the angular momentum J of the black hole, which are defined asymptotically as 
conserved charges of the spacetime. They satisfy the extremality bound \J\ < GM'^, whose 
saturation gives a degenerate horizon (zero surface gravity) , while otherwise the spacetime 
is nakedly singular. The Kerr-Newman black hole generalises the Kerr solution to 
include electric and magnetic charges, the uniqueness for a connected horizon still holding 
in the Einstein-Maxwell theory. The connectedness assumption can be dropped in the 
static limit, except in the charged case with degenerate horizons where there exist multi- 
black hole solutions of the Majumdar-Papapetrou type [12| . In the case of solutions with 
a cosmological constant, uniqueness has not been proven for asymptotically de Sitter (dS, 
A > 0) or anti-de Sitter (AdS, A < 0) spacetimes but the only localised black hole known 



in four dimensions is the generalisation of the Kerr case 13 . We shall later discuss how 
this simple picture changes drastically in higher-dimensional gravity, where there is a great 
variety of rotating black hole solutions. 

Observational evidence for the existence of black holes in the universe is indirect 
only, since no signal comes out of them according to general relativity. X-ray signals do 
come from accretion disks around black hole candidates, sometimes in binary systems. 
The geodesic motion of nearby bodies can also indicate the presence of a black hole. The 
size, mass and angular momentum of the candidates is inferred from that data, and the 
evidence for black holes is abundant. A supermassive black hole is thought to exist at 
the centre of most galaxies |14j , and even evidence for near-extremal black holes in binary 



1.2. QUANTUM DESCRIPTION 9 

systems has been reported, for instance |J| > 0.98 GM^ in |15|. Quantum evaporation of 
black holes, to be discussed in the next Section, opens the possibility of direct evidence, 
but the effects are too small for large astrophysical black holes. Very small astrophysical 
black holes, which would provide a definite signature, have not been detected. 

Hereafter, we shall use natural units c = G = h = kB = i, where fcs is Boltzmann's 
constant, implicit in the discussion of thermodynamics. 



1.2 Quantum description 
1.2.1 Black hole thermodynamics 

The laws of black hole mechanics have a close analogy with the common laws of thermo- 



dynamics |16 - 18 . This observation and the discovery of Hawking radiation 19 , 20 paved 



the way for progress in the quantum understanding of spacetime. 

A small variation in the mass M of a rotating black hole with angular momentum 
J, and say electric charge Q, satisfies 

dM = —dA + ndJ + (t>dQ , (1.2) 
87r 

where k is the surface gravity, A is the area of the event horizon, Q is the angular velocity 
of the horizon and (f) is the electric potential on the horizon. The analogy with the common 
first law of thermodynamics is clear since M is the conserved charge associated with the 
time-translation symmetry, i.e. the energy. J and Q are identified with "particle numbers" , 
and Q and </> are identified with "chemical potentials" , which are constant on the horizon, 
as required by equilibrium. We are left with the identifications of the temperature and 
the entropy of the black hole, T = ok/Stt and S = A/ a, respectively, where a is a positive 
constant. The zeroth law of thermodynamics is the statement that the surface gravity k 
is constant on the horizon. The second law is the area law of classical general relativity, 
stating that the event horizon area A never decreases [2l|. When matter is considered. 
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both the black hole entropy and the total entropy - black hole plus matter outside - 
must be non-decreasing classically. (Hawking radiation, a quantum effect, may cause the 
horizon area A to decrease, but this is compensated by the entropy of the radiation so 
that the total entropy never decreases.) To complete the four laws of thermodynamics, 
the third law says that it is not possible to make n vanish through a finite-time physical 
process. 

Although M, A, J and Q play the role of the extensive thermodynamic variables 
in the common applications of the first law, they have different scaling dimensions. For 
asymptotically flat solutions in four spacetime dimensions, the first law and the scaling of 
those quantities imply the Smarr relation [22] , 

M = 2(^^+ nj^ +<j)Q. (1.3) 

The relation fails for solutions with a cosmological constant, since an independent length 
scale is introduced. 

The laws of black hole mechanics result from the structure of the Einstein equations 



and from the fact that event horizons of stationary black holes are Killing horizons 17 
A Killing horizon is generated by a Killing vector K which is null on the horizon. For 
a four-dimensional black hole, K = dt + , where dt and 9^ are the stationarity and 
axisymmetry Killing vectors, respectively. We will later discuss specific aspects of higher 
dimensions. The first law has actually been shown to hold for any covariant Lagrangian 
theory of gravity (e.g. arising as a higher-derivative curvature correction to general relativ- 
ity) if the Bekenstein-Hawking entropy is substituted by the so-called Wald entropy 23p4 



The corresponding second law, however, remains an open problem. Let us also point out 
the result of Jacobson |25| stating that, while the first law is a consequence of the equations 
of motion, the converse is also true if we formulate the first law conveniently. Indeed, the 
Einstein equations can be derived by requiring that energy- momentum fluxes 5 Q across all 
local Rindler horizons (causal horizons of uniformly accelerated observers) through each 
spacetime point are reversible, i.e. 6Q = KdA/Sir. 
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The considerations above are purely classical. The temperature of a black hole, 
which emits nothing classically, can only be understood quantum-mechanically. The break- 
through by Hawking |19','20l was to use quantum field theory on a curved background to 
show that black holes behave like black bodies in the usual thermodynamic sense, emiting 
radiation with a thermal spectrum. The Hawking temperature is given by T = k/Itt, i.e. 
a = 4. The classical limit /i — ?• is clear if we do not use natural units: ksT = Hk/Ittc — )• 0. 
A stationary observer detects a thermal bath with temperature T/V, where V = \/— 
is the redshift factor = dt). The radiation is the effect of acceleration, characteristic of 
flat space too: Rindler observers, i.e. observers with uniform acceleration a, detect Unruh 
radiation with T = a/lir [26| . 

One of the great puzzles in this picture is what happens to the information of 
an object that enters the black hole. Since the radiation emitted by the black hole is 
thermalised, this information seems to be irremediably lost, violating unitarity ^27J. The 
rate of evaporation, i.e. mass loss through Hawking radiation, is given by the Stefan- 
Boltzmann law, dM/dt ~ -ylT"^ ~ -M'^, since A ~ and T ^ M ^. Evaporation 
is thus unimportant for large astrophysical black holes, but it dominates the behaviour of 
very small black holes, and it would provide a distinct observational signature in case they 
would be produced in accelerator experiments. In the final stages of evaporation, as the 
black hole approaches the Planck size, a full quantum description is unavoidable. Such a 
description must solve both the information paradox and the singularity problem. 

The discussion above implies that the entropy of a black hole, known as the Beken- 
stein-Hawking entropy, is given by 5 = A/ A, which strongly suggests that the quantum 
degrees of freedom of the black hole are effectively distributed over a surface, rather than 
a volume. This crucial observation is the basis of the holographic principle, proposed by 't 
Hooft |28| and Susskind [29|. This principle says that quantum gravity in a given volume 
should be described by a theory on the boundary of that volume, analogously to a common 
planar hologram that generates a 3D image. The number of degrees of freedom is thus 
drastically smaller than the naive expectation. 
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The Euclidean path integral approach to black hole thermodynamics, on which this 
thesis heavily relies, will be reviewed in Chapter [3] 

1.2.2 String theory and the AdS/CFT correspondence 

Only a theory of quantum gravity can identify the microscopic degrees of freedom which 
give rise to the Bekenstein-Hawking entropy. Moreover, according to the holographic 
principle, this theory should admit two equivalent formulations: one as a bulk theory, in 
which semiclassical general relativity arises explicitly as a low energy limit, and one as a 
boundary theory, in which the description of gravity is not explicit but the distribution of 
degrees of freedom is more conventional. 

String theory, a promising candidate for a theory of quantum gravity, has made 
progress on these two challenges coming from black hole thermodynamics: (i) it has pro- 
vided a counting of microscopic states for specific classes of black holes leading to the 
Bekenstein-Hawking entropy, first reported in |30|, and (ii) it has provided a concrete re- 
alisation of the holographic principle, which is the anti-de Sitter / conformal field theory 



(AdS/CFT) correspondence 31 32]. We shall make a brief comment although this is not 



a required background for the research described in this thesis. 

The idea of string theory is that the universe is composed of small strings, which 
resemble point particles at low energies, and that the worldsheet of a string (the (1+1)- 
surface spanned by its time evolution) is quantised. The consistency of the quantum field 
theory of the worldsheet has drastic consequences: the introduction of supersymmetry (to 
avoid the existence of tachyons and naturally include both bosons and fermions in the 
spectrum) and, more remarkably, the existence of ten spacetime dimensions (to ensure 
Lorentz invariance). Agreement with the commonly perceived four dimensions requires 
that six spatial dimensions are compactified on a very small scale. To offer in return, the 
theory: contains a graviton-type string state in its spectrum; requires that, at low energies 
and curvatures, the spacetime background on which the worldsheet moves satisfies the su- 
pergravity equations of motion (to ensure conformal invariance of the worldsheet theory) ; 
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and has a natural high energy cutoff scale (the string length) since interactions are not 
point-like, which solves the non-renormalisability problem of general relativity. The dras- 
tic consequences mentioned above can also be seen in a different light: supersymmetric 
field theories have appealing properties, and the compactification of the extra-dimensions 
provides a geometric framework to understand how the various parameters of the Standard 
Model of particle physics can arise from a simpler fundamental theory. There are five dif- 
ferent types of consistent string theories (types I, IIA, IIB, heterotic SO{32) and heterotic 
£"8 X Es), some of which have been shown to be related by dualities making them equiv- 
alent. They have all been conjectured to arise as special cases of an eleven-dimensional 
theory called M-theory. 

The quantum understanding of black holes provided by string theory comes mainly 
from the introduction of Dp-branes, which are extended objects with p spatial dimensions 
on which the endpoints of open strings are restricted to move. A gauge theory arises 
as the low energy description of these open strings. From a dual gravity perspective, 
Dp-branes are extremal black branes (i.e. black holes with infinitely extended horizons 
along p directions). The wrapping of D-branes on compact spaces can reproduce black 
holes in lower dimensions, whose quantum properties can then be derived from those 
of the D-branes. Strominger and Vafa [30] counted the degeneracy of D-brane states 
corresponding to microstates of a five-dimensional class of extremal black holes, leading to 
the first microscopic derivation of the Bekenstein-Hawking entropy. We should point out 
that, while D-branes are intrinsically string theory objects, it has been argued that the 
computation of the entropy may actually not rely on the string theory input, but rather 



on the symmetry of the near-horizon geometry or of the low energy wave equation [33 , 34 
(see (35] for potential issues). 



An outstanding development brought about by the understanding of D-branes is the 
AdS/CFT correspondence proposed by Maldacena [sT]. The original conjecture states that 
type IIB string theory with AdS^ x boundary conditions is dual to A/" = 4 SU{N) super- 
Yang-Mills theory defined on R x S^. This is a realisation of the holographic principle, 
since we have a gravity theory (string theory) in a five-dimensional spacetime {AdS^, 
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after compactification on S^) proposed to be equivalent to a four-dimensional field theory 
without gravity. In particular, the semiclassical supergravity limit on the gravity side 
corresponds to the strongly coupled limit on the field theory side. Black hole physics 
in AdS can describe strongly coupled gauge theories at finite temperature (the Hawking 
temperature of the black hole) and vice-versa. The conjecture is expected to extend to 
other cases of string theory or M-theory with AdS (or "almost AdS" ) boundary conditions, 
the dual field theory being defined on the AdS boundary. 



1.3 Higher-dimensional solutions 

In this Section, we motivate the study of higher-dimensional gravity and review briefly the 
literature on the existence and properties of black hole solutions, leaving the important 
issue of stability for the next Chapter. Our main focus will be on vacuum solutions to 
general relativity, with particular attention to Myers-Perry black holes, which will be 
analysed in the last part this thesis. See [36] for a comprehensive review, although some 
developments are more recent. 

String theory is the original and most important motivation for the study of higher- 
dimensional gravity. Black holes in spacetimes with up to ten dimensions are part of 
string theory (or eleven dimensions, for M-theory). We mentioned in the last Section that 
one of the theory's major successes, the microscopic derivation of the Bekenstein-Hawking 
entropy, was obtained for a class of five-dimensional black holes jSOl. The AdS/CFT cor- 



respondence 31,32] and its applications further motivate this study. The correspondence 
implies that the dynamics of Z)-dimensional field theories at finite temperature and that 
of (1) + l)-dimensional black holes are equivalent. Gravity solutions have been constructed 
d la carte to describe strong-coupling features of field theories, in an effort to learn some- 



thing about quantum chromodynamics |37] and condensed matter physics 38 . Other 
scenarios, which try to explain the hierarchy problem (the weakness of 4D gravity when 
compared to the Standard Model interactions) by introducing extra-dimensions, possibly 
in string theory embeddings, open the possibility that black hole production in particle 
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collisions is within future experimental reach, as reviewed in ^39j . 

Higher-dimensional gravity is also important in its own right. A new insight is 
gained by studying gravity using the spacetime dimensionality D as a parameter. Indeed, 
many properties of black hole solutions are specific to four dimensions. 

In Newtonian gravity, the attractive gravitational force is supressed with the radial 
distance as r~^^^^\ while the repulsive centrifugal force is supressed as (for given 
mass and angular momentum) independently of D since it acts on a plane of rotation. 
This is why stable planetary orbits are not possible for D > 4. One might expect that the 
physics of black holes, especially when there is rotation, will also be different. In fact, there 
will be black holes with different horizon topology and also disconnected horizons. Black 
holes with an arbitrarily large angular momentum for a given mass can also be found. 
These are closely related to the existence of solutions with extended horizons (absent in 
four dimensions without cosmological constant). 

Another obvious contrast is the number of rotation planes. The Cartan subgroup, 
i.e. the maximal Abelian subgroup or maximal torus, of SO{3) is ^7(1), which means 
that there is a single independent rotation plane in D = 4. The associated rotational 
Killing vector = dfj, defines the angular momentum as an asymptotic conserved charge, 
J ~ fq2 *d$, . However, the Cartan subgroup of SO{D — 1) is f7(l)"', so that there 
are n = [{D — l)/2j (where [J stands for the smallest integer part) independent rotation 
planes in D — 1 spatial dimensions. These allow for the definition of n independent angular 
momenta, Jj ~ Is^~^ ■ (Notice that only one of the ^j's must be a Killing vector of 
the spacetime according to the rigidity theorem, to be discussed below.) 

As we mentioned in Section |1.1[ the uniqueness of the Kerr black hole among sta- 
tionary asymptotically fiat vacuum black holes with a connected horizon is a remarkable 
feature oi D = A. In higher dimensions, the generalisation of the Schwarzschild solution, 
obtained by Tangherlini [40], is also the unique static vacuum black hole f4T]. Uniqueness 
extends to the higher-dimensional Reissner- Nordstrom black hole for the Einstein-Maxwell 
theory in the static non-degenerate case |42|, while in the degenerate case there exist the 
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higher-dimensional Majumdar-Papapetrou solutions, as in four dimensions. It is the in- 
clusion of rotation that changes this simple picture. 



1.3.1 Myers-Perry solutions 

The Myers-Perry solution |43| is the natural generalisation of the Kerr black hole to higher 
dimensions, and (so far) the only exactly known rotating solution in D > 5. It is a rather 
non-trivial generalisation, due to the possibility of rotation in different planes. Fortunately, 
as the Kerr solution, it can be written in the so-called Kerr-Schild form, 

9^lu = 'n^u + H{x^)k^ky , (1.4) 

where r/^i^ is the Minkowski spacetime metric and is a null vector with respect to g^^. 
The fact that k^ is then null also with respect to rj^u allows for a simplification of the 
Einstein equations similar to linearisation. 

We present the solution here in more detail because it will be useful in this thesis. 
Let us take D = 2n + 1 + e, where n = [(D — l)/2j is the number of angular momenta 
and € = (D — l){mod 2} is {D odd) or 1 (D even). There are n azimuthal angles (pi 
and n + e direction cosines /Xj obeying = 1. In Boyer-Lindquist coordinates, the 

metric can be written as 

i=l i=l 

+ ^ 1.5 

V — 2m 



where the last term is present only for even D, and 



n+e 9 n 



U{r, = ^ n(^' + «') ' ^(^) = n(^' + • (1-6) 

i=l * j=l i=l 

The event horizon is located at r = r_|_, where r_|_ is the largest root of V{r) — 2m = 0. 
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The Kerr black hole case {D = 4) is included as n = e = 1. 

The solution is parametrised by n + 1 length scales: the horizon radius r+ and the 



n rotation parameters Oj. The horizon area is 



where Ad- 



+ i=i 



2 7r(^-i)/2 
n{D-l)/2] 



:i.7) 



is the volume of a unit-radius {D — 2)-sphere. The surface gravity k and the angular 
velocities on the horizon are given by 



K = r_|_ 



E 



1 



2 - e 



rl + af ■ 



:ii 



The asymptotic charges, the mass M and the angular momenta Jj, which uniquely specify 
a solution, are given by 



[D - 2)Ad-2 

M = m . 

Svr 



Ji 



D-2 



aiM . 



:i.9) 



Extremality, for which the temperature T = k/2tt vanishes, occurs for 



Er 

i=l 



2 - e 



+ {MJi/Af 



1 D odd 
1/2 D even 



(1.10) 



In D = 4, this introduces the bound J < on the angular momentum, beyond which 
the Kerr solution is nakedly singular. Does such a bound exist in higher dimensions? Let 
us pick up a direction through the origin in the space {Ji} for a fixed mass. We can readily 



see from (1.10) that a direction with regular extremal limit requires Ji ^ Vi. Otherwise: 



D even : There is no bound when one or more of the angular momenta vanish. 



since it is impossible to satisfy (1.10), which leaves all the remaining angular momenta 
unbounded for a given mass. 

• D odd : There is no bound when two or more of the angular momenta vanish, 
leaving all the remaining angular momenta unbounded for a given mass. However, there 
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D = 5: 



D = Q: 




D = 7: 



D = 8: 





Figure 1.1: Shape of the parameter space of angular momenta for fixed mass in Myers- 
Perry black holes for D = 5, 6, 7, 8. 



is a bound corresponding to (1.10) when a single angular momentum vanishes. The satu- 



ration of the bound gives a nakedly singular solution, rather than an extremal black hole, 
since it is required that A ^ (no event horizon). 

Figure [LT] presents some examples. In D = 5, there are two angular momenta and, 
if one of them vanishes, there is a bound on the other for a given mass, the saturation 
of which gives a naked singularity (the four corners in Figure [lT] ). In L> = 6, there are 
two angular momenta and, if one of them vanishes, there is no bound on the other. In 
D = 7, there are three angular momenta and two cases occur: if only one vanishes, there 
is a singular bound (the intersections of the extremality curve with the Ji,2,3 = planes 
in Figure [lT]); but if two vanish, there is no bound (along the axes). In Z? = 8, there are 
three angular momenta and, if any of them vanish, there is no bound on the other two. 



The fact that there are solutions with arbitrarily high angular momentum for a 
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given mass allows for different regimes, e.g. in the D > 5 singly-spinning we can have 

D-2 D-2 

\J\/Mo-3 ^ 1 but also \J\/Mo-s ^ i. This hints at a rich phenomenology, as we shall 
confirm. 

The isometry group of the Myers-Perry family is M x [/(l)" (up to discrete factors), 
where M corresponds to the time translations and C/(l)" to the rotational symmetries. 
In Chapters [7] and [Sj we will consider particular cases of the Myers-Perry family with 
enhanced symmetry, namely solutions with a single spin (say Ji ^ 0, Ji = \fi > 1), which 
have isometry group M x [/(I) x SO{D — 3), and solutions with equal spins (Jj = J 
which have isometry group R x [/ (n) . 

The Myers-Perry family has been extended to asymptotically (A)dS spacetimes in 
(4^ {D = 5) and (45l , |46| {D > 5). It is the only exactly known family of localised black 
holes with that asymptotic behaviour. 

1.3.2 Plethora of black holes 

There are no solutions with extended event horizons in four-dimensional general relativity 
unless a cosmological constant is present A reason for this is that no asymptotically flat 
black holes exist in D < 4. Notice that the direct product of two Ricci-flat manifolds is also 
a Ricci-flat manifold, so that solutions with extended horizons are trivially constructed 
from lower dimensional black holes. Consider the metric of a (D + A^)-dimensional vacuum 
black brane which uniformly extends a L)- dimensional asymptotically flat vacuum black 
hole with metric gab along M^, 

^^brane — 9ahdx°'dx^ + dz ■ dz . (l-H) 

The simplest example is the Schwarzschild string, SchwarzD x M. Black branes need 
not be uniform and need not be obtained from a direct product. We will discuss how 



^Although black hole spacetimes in D = 4 which asymptote to AdS in all spatial directions are restricted 
to have spherical event horizon topology [47| , solutions with extended cylindrical, planar and hyperbolic 
horizons exist 48 -51 . 
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families of non-uniform branes bifurcate from uniform branes due to linear instabilities 
with stationary threshold. 

The existence of black strings is suggestive. Take a segment of a black string and 
bend it to make it circular like a ring, introducing also rotation on the plane of the circle so 
that the centrifugal force may balance the gravitational pull. Emparan and Reall showed 
in [52] that there is indeed an asymptotically flat black ring solution in D = 5 (generalised 
to have two angular momenta by Pomeransky and Sen'kov [53j). The topology of the black 
ring event horizon is x 5^. Topological restrictions, such as Hawking's proof |54| that 
the horizon topology must be spherical in Z) = 4, are much weaker in higher dimensions 
(see e.g. [55]). 



The existence of black rings shows that uniqueness is lost generically for rotating 
black holes. There is a region in the parameter space of = 5 black holes, near the corners 
of the Myers-Perry parameter space of Figure where three solutions with connected 
horizons exist: the Myers-Perry black hole and two black ring solutions. Let us concentrate 



on the case with a single spin, shown in Figure 1.2, There are two black ring solutions 
meeting at a cusp in the graph: the fat ring has smaller entropy and finite J range, and 
the thin ring can have arbitrarily large J. The existence of these two solutions shows that 
uniqueness fails even for a given horizon topology. In the plot, the Myers- Perry black hole 
and the fat ring meet at a zero-area point where both cases give the same nakedly singular 



solution. Figure 1.3 presents lower-dimensional pictures of a fat ring and a thin ring. 

One may also wonder whether stationary vacuum solutions with disconnected hori- 
zon components can exist, which is believed to be impossible in = 4. Take a black ring 
(BR), put it around a Myers- Perry (MP) black hole, and then make the black ring rotate 



faster to balance the gravitational pull, in a Saturn- like configuration as in Figure 1.3 
Elvang and Figueras (56) constructed the singly-spinning black Saturn solution, which has 
a disconnected event horizon: an S^ component and an S^ x S^ component. There are re- 
gions of the parameter space where the solution can have the same asymptotic charges as a 
Myers-Perry black hole or a black ring. However, the solution possesses a continuous type 
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A 




M3/2 



Figure 1.2: Entropy S versus angular momentum J for fixed mass M, for the singly- 
spinning black rings (red) and the D = 5 Myers-Perry black hole (blue). Notice the 
zoomed detail in the top-right corner, where the fat and the thin ring solutions coincide 
at the cusp. 





Figure 1.3: Lower-dimensional pictures of the event horizons of a fat ring, a thin ring and 
a black Saturn (left to right). 

of non-uniqueness given the asymptotic charge^ masses and angular momenta can be 
assigned to each horizon component by Komar integrals, and they can vary while keeping 
the asymptotic mass M = Mmp' + Af'BR' and angular momentum J = J'mp' + -^'BR' fixed. 
The case J'mp' = — <^'BR' shows that non-static J = solutions are possible. Furthermore, 
there are interesting frame dragging effects between the two components. Solutions have 
also been constructed which represent two concentric black rings rotating on the same 



plane (di-rings) |58| or on perpendicular planes (bi-rings) [59 60 



These solutions were obtained with powerful solution generating techniques for D 



Continuous non-uniqueness was first obtained for black rings with a dipole cliarge 57 
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5, reviewed in |36|j, which assume the existence of two rotational Kilhng symmetries. There 
is uniqueness in this class, given the asymptotic charges and the so-called rod structure that 
generates the solution [61]. The failure to construct an analogous scheme in D > 5 is the 
reason why Myers-Perry solutions are the only rotating solutions known exactly. However, 
thin black rings with large rotation in D > 5 have been constructed using approximate 
methods |62| . These methods were later put into a more systematic framework, where 
black holes with quasi-extended horizons are given approximate solutions called blackfolds 



[63 - 65 . The approach allows for the construction of a variety of solutions with different 
horizon topology. This is important since, as we mentioned, topological restrictions are 
much weaker in higher dimensions. 

An important recent development is related to the rigidity theorem, generalised 
from D = 4 [3] to higher dimensions in [55 66 , which requires stationary (non-extremal) 
black holes to have at least one rotational symmetry. Until recently, not a single black 
hole had been found with only one rotational symmetry, though such solutions had been 
conjectured to exist [67]. The blackfold approach gave the first example, as it allowed for 
the approximate construction of helical black rings in any D > 4 which possess a single 
rotational symmetry [65] . The perturbative numerical work of |6^, presented in Chapter [sj 
of this thesis, provided the first evidence of such solutions with spherical horizon topology, 
in Z? = 9 (and all higher odd dimensions, we expect). A single rotational symmetry was 
also obtained in all even D > 8 for near-horizon geometries which give an infinite class of 
distinct horizon topologies and are consistent with the requirements of asymptotically fiat 



or AdS extremal black holes 69 



The goal is to uncover the phase diagram analogous to Figure 1.2 but for all solutions 
in any D > 4. Notice how the fat ring and the Myers-Perry solution coincide at one point 
in that Figure. The black Saturn in thermodynamic equilibrium (horizon components 
with the same temperature and angular velocity) appears in this plot in a similar way to 
the black ring, also joining the Myers- Perry solution at the singular point; see Figure 2 



of 70 . In D > 5, this singularity is "resolved" since the angular momentum is unbounded 



for the singly-spinning Myers- Perry black hole (recall Figure 1.1). The expectation [62[|71 
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is that the connection between phases at the singular bound in D = 5 will correspond in 
D > 5 to bifurcations of regular solutions which interpolate between Myers-Perry black 
holes and phases such as the black ring and the black Saturn. We shall elaborate on this in 
Chapter [7| which is based on the numerical work of |72|. The results imply non- uniqueness 
in L> > 5 given spherical horizon topology. 

It would be very interesting to see how the picture described here would change 
with the inclusion of a cosmological constant, especially a negative one having in mind the 
AdS/CFT correspondence. We mentioned that the Myers-Perry solution is the only one to 
have been extended to asymptotically (A)dS spacetimes. However, black rings have been 
constructed approximately in AdS [73] and there is no reason why a variety of solutions 
should not exist in this case too. 



1.3.3 Thermodynamics 

The existence and ~ as we shall see - the stability of black hole solutions depend critically 
on the dimensionality of spacetime. However, the laws of black hole thermodynamics are 
a consequence of the structure of Einstein's equations in any dimensions. 

For charged black holes, the first law in higher dimensions is 

dM = ^ + J^S^^djf") + ct^'f^Q'f^ , (1.12) 

where we have a sum over different planes of rotation - index i -, a sum over charges - 
index i - and a sum over different event horizon components - index (a) . This extension of 



( 1.2 ) can be derived with the standard procedure [I7j, by considering that each component 
is a Killing horizon of K^"'' = dt + . Notice that the individual angular momenta 

Jj^"^ and charges Q^"^ are not asymptotically defined conserved quantities in the case of 
disconnected horizon components. Even in the case of a connected horizon, a dipole charge 



Q*'^ has no net asymptotic contribution, despite being included in the first law 57,74 . 
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The Smarr relation, valid for asymptotically flat spacetimes, reads 

iD-3)M = {D-2) (^^L^^ + + ae <pfUQf^ , (1.13) 

where a£ is the scaling dimension of the charge (it differs for dipoles and conserved charges) . 
Thermodynamic equilibrium between components of the horizon requires 

k('^) = «, n^^ = ni, (l>P = <f>i, (1.14) 

which reduces the space of solutions, e.g. the contributions J^^"^ from different horizon 
components to the total conserved angular momenta Jj are fixed. The first law becomes 
simply 

dM = -^dA + Q,idJi + (t)idQ^. (1.15) 



Chapter 2 



Stability of black holes 



The construction of a stationary black hole spacetime satisfying the equations of motion 
is not enough to appreciate the physical significance of the solution. In a "realistic" 
scenario, we have to know whether or not the spacetime is robust for small perturbations 

of the geometry and matter fields. If not, a probe, say a particle moving on the black 
hole background, may cause a disruptive backreaction; and the possibility of dynamically 
forming such a black hole through a physical process, such as gravitational collapse, is put 
in doubt. 

Nevertheless, unstable solutions are not devoid of physical significance. Instabili- 
ties have characteristic timcscalcs, and a solution may be "sufRcicntly stable" for shorter 
timcscale phenomena. Furthermore, the manifestation of an instability is itself of great 
interest. The onset may indicate a phase transition through a bifurcation, as we shall 
see. Information on the timescales is very useful especially in the cases where different 
instabilities are present, since it may give hints on what the final state is, given that the 
full time evolution is often beyond our technology. 

The complete description of black holes - stability included - must take into account 
their quantum thermodynamic properties. It is well known that the Schwarzschild black 
hole is classically stable at the linear mode level and yet it has a negative specific heat. 
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which signals a (local) thermodynamic instability. Standard thermodynamic arguments 
can also identify (global) phase transitions between different spacetimes. 

In this Chapter, we will discuss both the classical and the thermodynamic stability 
of black holes. It is one of the goals of this thesis to explore their connection. 

2.1 Classical stability 

We are concerned here with the stability of stationary black hole spacetimes as solutions to 
a system of differential equations of motion, namely the Einstein equations, or extensions. 
The literature often consists of linear stability studies. Just as happened for uniqueness, 
stability does not always hold for asymptotically flat vacuum black holes in D > 4 when 
rotation is considered. In fact, we shall argue that uniqueness and stability are intimately 
related. Is there uniqueness of classically stable black holes? This remains a major open 
problem. 

2.1.1 Statement 

Classical stability is based on the initial value problem of general relativity |2|. Brushing 
aside the intricate technical aspects, the stability problem can be put as follows. Take 
a Cauchy surface S of a stationary solution {M,gab) and perturb the initial data C 
on S. That is, instead of the initial data C that develops into the stationary solution 
from S, {Ms, 9ab)j take smooth initial data C which is "close" to C and satisfies the 
Einstein constraint equations. Now consider the future development {Ms^Oab) of C. If, 
for any "small" perturbation (C close enough to C), the future development {My;, dab) is 
"close" to - or actually "approaches" - {M.T:,gab), then the stationary solution {M.,gab) is 
stable. For instance, starting with the Kerr spacetime, the perturbed solution (A^s^^ab) 
could become a Kerr solution at late times, but with slightly different mass and angular 
momentum. Otherwise, the stationary spacetime is unstable, and the time evolution in 
{AiT,,'gab) will lead "away" from (A^s^ffab) for a class of perturbations. The inclusion of 
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matter fields presents no problem of principle, the initial value problem being set for 
the solution (M., gab,^s)- 

For black hole spacetimes, only the development of the data on E between the future 
event horizon and spatial infinity is important, as the evolution inside the horizon 
cannot affect the evolution outside. We are therefore interested only in the stability of the 
black hole exterior. 

The non-linearity of the Einstein equations makes the stability problem very chal- 
lenging. A rigorous proof of non-linear stability for small perturbations has been achieved 
only for the simplest of spacetimes, Minkowski space, by Christodoulou and Klainer- 
man 



75 . For black holes, apart from numerical approaches, only the linear problem has 
been studied. In this case, one considers perturbations hat = gab — dab of a black hole 
background Qab, subject to boundary conditions to ensure the regularity of the perturbed 
spacetime Tjab- The perturbations satisfy the linearised Einstein equations, 

{/^Lh)ab - VaVf,h% + 2V(,V^/ife), +gab{- V^V^h.^ + □/I'e + ^""^cd) 

+ i2A -R)hab = IdnSTab, (2.1) 



where is the Lichnerowiz operator, defined as 

{ALh)ab = -Wchab - '^Rab'^h^ + '^Rlah)c , (2-2) 

and 5Tab is the linear perturbation of the energy-momentum tensor. The gauge ambiguity 
in hab must be dealt with, either by fixing the gauge or by considering gauge-invariant 
quantities. Stability requires that these perturbations remain bounded in their time evo- 
lution, the bound depending on the initial data. This usually means that an appropriate 
positive definite functional J-[h] (t) of the perturbation, defined on a time slice outside 
the horizon, obeys J-[h]{t) < CJ-[h]{tQ), where the initial data is defined on Sj^, t > to, 
and C is a constant depending on the parameters of the background geometry. It may 
also be possible to show "pointwise" boundedness, e.g. for the perturbation as a function 
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of the radius in a spherically symmetric spacetime, as opposed to a functional over a time 
slice. When considering the implications to the non-linear problem, stronger conditions of 
time decay - and not just boundedness - are typically required at the linear level (this is 
because some integrals over time should also be bounded). Notice, however, that the solu- 
tion approached at late times by the time decay of perturbations will have small variations 
in its parameters, such as mass and angular momenta, with respect to the initial solution. 
The time decay of the scalar wave equation is an important first step in understanding 
the decay of metric perturbations. 

The symmetries of the black hole background can simplify the problem consider- 
ably. For stationary backgrounds Qab, a linear perturbation is a superposition of Fourier 
modes h^^^ ~ e~^^^, and stability implies that, when analysing each mode separately, only 
modes with Im(a;) < are allowed. Spatial symmetries of the background also allow for 
decompositions, e.g. into spherical harmonics when there is spherical symmetry. Mode 
analyses are much simpler and very effective, as we shall see in this thesis. In the case 
that no unstable modes are found, the approach falls short of the mathematical rigour 
required to establish linear stabilityj^ However, the detection of an instability is clearer, 
since the initial data may be fine-tuned to reproduce a particular unstable mode. 



2.1.2 Review 

Let us start with the static case (for a more detailed review, see |77| ) . In D = 4, the 
linear mode stability of the Schwarzschild black hole was established by [78| - [8T] . Spherical 
symmetry allows for a decomposition of linear perturbations into an even-parity class 
(constructed from scalar harmonics of S^) and an odd-parity class (constructed from 



vector harmonics of S" ), as shown by Regge and Wheeler 82 . Each class of perturbations 



is governed by a master equation: the Regge- Wheeler equation in the odd-parity case [82] 



and the Zerilli equation in the even-parity case 83 . After a Fourier decomposition of the 



^Because an infinite sum of modes wliicli decay in time may not decay in time. Furtliermore, tliere are 
examples, sucli as coloured black holes [76], where regular initial data can be constructed from unstable 
modes which are irregular per se, and thus usually discarded. 
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time dependence, each becomes a Schrodinger-type ODE allowing for no unstable regular 



modes. A rigorous approach avoiding that decomposition is put forward in 80,81 , where 
the boundedness of linear perturbations is shown. Analyses of Schrodinger-type master 
equations have established the mode stability of the non-degenerate Reissner-Nordstrom 
black hole [841, using the Newman-Penrose formalism |85|, and the Schwarzschild-(A)dS 



black hole 86 , 87 . In the latter case, rigorous studies of the boundedness and decay of 



solutions to the scalar wave equation are available, see e.g. the review 88 for dS and [89] 
for AdS, where a scalar field mass obeying the Breitenlohner-Preedman bound 90 is 
considered 13 

In higher dimensions, the procedure has been applied with success to establish the 
mode stability of Schwarzschild black holes for arbitrary D 92||93 . Stability extends to the 



D = 5 Reissner-Nordstrom case, but no such proof is available for D > 5. Schwarzschild- 
dS black holes have been shown to be stable for D < 6. While Reissner-Nordstrom-dS 
black holes are stable for D < 5, an instability was found for sufficiently large charge in 
D > 6 [94|95| . In the asymptotically AdS case, not even the stability of Schwarzschild- AdS 



black holes is established for D > 4. 

Now we consider rotating asymptotically flat black holes, starting with D = 4. For 
the Kerr spacetime, which is codimension-2, Teukolsky obtained the remarkable result 
that the linearised Einstein equations can also be reduced to a single separable mas- 
ter equation 96 , using the Newman-Penrose formalism [85jj^ Whiting showed that the 



Teukolsky equation does not allow for unstable Fourier modes |98j. Such a separability 
of the equations has not been achieved in the Kerr-Newman case, which remains an open 
problem. Ref. [99| reviews the rigorous study of the scalar wave equation for slowly- 
rotating (|J| <^ M^) small-charge (|Q| <^ M) Kerr-Newman black holes, avoiding the 
use of separability properties. That approach has provided proofs of the boundedness 



^In asymptotically AdS spacetimes, the presence of a Cauchy horizon demands more care in setting 
boundary conditions at spatial infinity, in order to have a well-defined initial value problem [90[|91| . Neg- 
ative mass-squared scalar fields can have a well-defined dynamics if the mass obeys the Breitenlohner- 
Freedman bound (90|. 

^The Teukolsky equation considers radiative modes, thus excluding stationary modes, which consist 



only of trivial mass and angular momentum variations, as shown by Wald 97 . 
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and decay of solutions to the scalar wave equation, and has also strengthened the results 



of [80,81 for the Schwarzschild case. Let us also point out that numerical studies of black 
hole formation (e.g. |6,'7]) provide important evidence for stability in the full non-linear 
regime. Astrophysical black holes should be unique and stable. 

The picture of stability for higher-dimensional rotating black holes is very different. 
The first discovery of an instability in higher dimensions concerned solutions which, despite 
being static, possess an extended horizon. Gregory and Laflamme showed numerically that 



black branes of the type ( 1.11 ) constructed with a D-dimensional Schwarzschild black hole 



are unstable (lOOj . Consider a transverse and traceless (TT) linear perturbation of the 



metric (1.11) of the form 



'^■^brane = gah{x)dx°' dx^ + dz ■ dz — )• '^■^brane + e^^^ hah{x)dx°' dx^ , (2.3) 

The condition is equivalent to hat being TT with respect to the D-dimensional black hole 
metric Qab, 

h\ = 0, V''hab = 0. (2.4) 



The linearised Einstein equations (2.1) reduce to 



{ALh)ab = -k'hab, (2.5) 

where is the Lichnerowicz operator for the black hole background, and k"^ = k ■ k. 
Hence perturbations with non-zero k correspond to negative modes of A/,. The boundary 
conditions are that hab should be regular on the future horizon and vanishing at infinity. 
Gregory and Lafiamme studied time dependent modes hab e'"* which preserve the spher- 
ical symmetry of the Schwarzschild Qab backgroundj^ They found exponentially growing 



solutions to (2.5) for A; < fc*, as shown in Figure 2.1 The instability pinches the branes 



along the extra-directions z, breaking the translational symmetry. The final point of the 



*The equations (2.51 reduce to a single ODE. In a Fourier decomposition hab oc e~'"*, the hermiticity 
of the ODE operator guarantees a real eigenvalue cj^, so that F = io; is purely real (and positive) for an 
instability. 
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Figure 2.1: Profile T{k) for the Gregory-Laflamme unstable mode of a Schwarzschild brane 
(-C'schwarz = 4), where Vg stands for the Schwarzschild radius. 



time evolution remains an open problem 101 . The stability properties may be improved 



with the inclusion of charge (see the review [102j ), according to the Gubser-Mitra conjec- 
ture, to be discussed in the next Section. 



Notice that, while the mode {k, T) = (0, 0) is a gauge mode [103 , the mode {k, T) = 
(A;*,0) is a static linear perturbation signaling the existence of a one-parameter family 
of static non-uniform branes. The existence of such a family of black strings was shown 
first in |104j , at the perturbative level, and then in [105 , with a non-linear numerical 
analysis. Along the family of non- uniform solutions, the pinches become greater and 



greater, supposedly until the string splits into an array of black holes; see [102 ,106 for 



reviews and [107 for more recent work. The idea is represented in Figure 2.2 The same 



reasoning should apply to the rotating non-uniform strings discovered in 108 and [68] , 
also associated with instabilities as verified numerically in [68[ . 

The fact that the event horizon is infinite along the z directions but finite along the 
black hole directions is the cause of the instability!^ Now, in higher dimensions, black hole 
horizons can be characterised by very different length scales. Say we have a horizon scale 
Lx along certain directions x and a horizon scale Ly <^ Lx along other directions y. A 
linear perturbation along the x directions with a length scale A satisfying Ly < X < Lx 



^Ref. |l09| proposed an analogy with the Rayleigh-Plateau instabihty of fluid tubes, which was made 
more precise in the context of the fluid-gravity correspondence 110[|111 . 
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Figure 2.2: Family of solutions from a uniform black string to an array of black holes. 



may behave like a Gregory-Laflamme unstable mode. Recall that the horizon of a Kerr 
black hole is a deformed sphere, but that there is a limit on the deformation which is set 
by extremality, | J|/Af^ < 1. We already mentioned that such a bound is absent for singly- 
spinning Myers Perry black holes in Z) > 5: we have the regime \ J\/Md-3 ^ 1, for which 
the horizon is approximately spherical, but also the regime |J|/Mo-3 ^ 1, for which 
some directions along the horizon spread, becoming quasi-extended. In the latter regime, 
the black hole is expected to be unstable, because the metric near the axis of rotation 
takes the form of a black-brane metric and should thus possess a Gregory-Laflamme-type 
Similarly, very thin black rings, satisfying |J|/Af2 ^ 1, have a horizon 



instability 



71 



with topology x S'^ for which the radius is much larger than the radius of the 
squashed 5^. Locally, a segment of the very thin ring looks like a boosted black string 



and is expected to be unstable to Gregory-Laflamme-type modes 52 



In Chapters [T] and [sj based on j35 72 , we will analyse this type of instability nu- 
merically for Myers- Perry black holes with a single spin and with equal spins. Our results 



for a single spin are consistent with the expectations of 62 , 71 regarding connections be- 



tween black hole families in D > 5. These connections arise in an analogous way to the 
bifurcation of non-uniform branes from uniform branes, i.e. through the stationary thres- 



hold mode (/c,r) = (A;*,0) of an instability, as in Figure 2.2 In the case of Myers- Perry 
black holes with equal spins, we find a similar instability in D = 9 near the extremal- 
ity bound. No clear geometrical understanding of the instability, by analogy with the 
Gregory-Laflamme case, has been proposed for equal spins due to the presence of the 
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extremality bound. However, a definition of an ultraspinning regime will be proposed, 
based on a connection with thermodynamics, and conjectured to be a necessary condition 
for this type of instability. The thermodynamic argument is only expected to hold for 
perturbations which preserve the stationarity dt and axisymmetry Qid^p- Killing vectors. 

Also very recently, Shibata and Yoshino detected an instability of rapidly-rotating 
singly-spinning Myers-Perry black holes, using a non-linear numerical method [112 113 



The instability occurs for non-axisymmetric deformations. Although the approach is non- 
linear, the instability probably holds at the linear level, but it should not have a stationary 
and axisymmetric threshold mode. Hence the thermodynamic argument for our ultraspin- 
ning conjecture does not apply. Indeed, no D = 5 Myers-Perry black hole falls into our 



ultraspinning class. Ref. 114 found no evidence of instability for D = 5 Myers- Perry 



black holes with equal spins. 

We mentioned above that very thin black rings should be unstable to Gregory- 



Laflamme-type modes along the ring circumference. Ref. 115 points out the possibility 
that this affects the entire thin ring branch, since all thin rings can fit several dangerous 
modes in their circumference. Furthermore, evidence is provided for the instability of fat 



rings for radial perturbations. Fat rings had been argued to be unstable in 1116 , based 



on the thermodynamic turning-point methodj^ The question remains whether there is 

3 

a portion of the thin ring branch with J/M2 « 0(1) which is stable. More generally, 
among the plethora of higher-dimensional asymptotically flat vacuum black holes, are 
Myers-Perry solutions with moderate rotation the only stable ones? 

To conclude, let us consider rotating asymptotically (A)dS black holes. In dS, the 
full stability problem remains open even in D = 4. In the Kerr- AdS/Myers- Perry- AdS 
case, however, it is known that an instability for non-axisymmetric perturbations occurs 
when 3i > 1, where £ is the AdS curvature radius. The instability is caused by 

superradiant modes which are amplified near the horizon and reflected back to it by the 



Such a connection between classical stability and thermodynamics is of a different nature than the one 
we will explore in this thesis, which relies on local thermodynamic stability in the spirit of the Gubser-Mitra 
conjecture. 
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AdS gravitational potential [117 -121 . When the rotation is slow, \Qi\£ < 1 Vi, these black 
holes should be stable [117 . Ref. 89 provides a rigorous study showing the boundedness 
of solutions to the scalar wave equation for slowly rotating Kerr-AdS black holes, when 
the scalar field mass obeys the Breitenlohner-Freedman bound |90|. 



2.2 Thermodynamic stability 



The thermodynamic description of black holes, outlined in Section 1.2 accounts for their 
quantum properties at the semiclassical level. Stability with respect to quantum effects 
can be studied, within the semiclassical approximation, using the well-established criteria 
of thermodynamic stability. 

Thermodynamic stability tells us how a system in thermodynamic equilibrium re- 
sponds to fluctuations of energy, temperature and other thermodynamic parameters. The 
conditions for stability are derived from the second law of thermodynamics. If they are 
violated, which often happens for black holes, the system is not in a preferred configu- 
ration. We distinguish between global and local stability. Because black holes, as other 
self-gravitating systems, are not extensive some care is required in interpreting these con- 
ditions. 

For definiteness, we will focus on vacuum black holes in thermodynamic equilibrium, 
for which the first law reads 

dM = TdS + VtidJi . (2.6) 



2.2.1 Global stability 

Global stability is concerned with the phase of a system corresponding to the global 
maximum of the total entropy. 

Consider a system in equilibrium with a thermodynamic reservoir at temperature T 
and angular velocities Jlj. If we allow for mass and angular momenta exchanges (AM 7^ 0, 
AJj 7^ 0) between the system and the reservoir, the relevant ensemble to describe the 
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system is the grand- canonical. The preferred phase of the system is the one that minimises 
the Gibbs free energy, 

G = M -TS -niJ^, (2.7) 

given T and ilj. Therefore, for two possible phases - 1 and 2 - satisfying Gi{T,0,i) > 
G2{T,Qi), phase 1 is unstable. Starting with phase 1, thermodynamic fluctuations will 
eventually lead to a phase transition to phase 2. To see this from the second law, the crucial 
point is that the reservoir is assumed to be very large so that no exchange affects its internal 
equilibrium. Hence, exchanges are reversible for the reservoir: AMfcs = TASj^cs+^i^Jires- 
Conservation of mass and angular momentum requires AMres = —AM and AJjres = 
— AJj. The second law then favours fluctuations obeying 

AS + ASros>0 ^ A5 - ^^'^ -J^'^^Ji > ^ AG<0. (2.8) 



If we restrict the exchanges such that AJj = 0, but AM ^ 0, then the relevant 
ensemble is the canonical. The preferred phase is the one that minimises the Helmholtz 
free energy, 

F = M-TS, (2.9) 

given T and Jj, since the second law implies that AF < is favoured. 

An isolated system (AM = AJj = 0) is described by the micro- canonical ensemble. 
The second law is direct: AS* > 0, the preferred phase of the system being the one that 
maximises the entropy S given M and Jj. 

In higher-dimensions, there are several black hole phases so that we may expect to 



find phase transitions due to the global instabilities. For instance. Figure 1.2 represents 
three phases in a certain parameter range. However, as we will show later, following [68] , 
asymptotically flat vacuum black holes are locally thermodynamically unstable. The most 
celebrated example of a global phase transition between locally stable phases occurs in 
asymptotically AdS spacetimes: it is the Hawking- Page phase transition between thermal 
AdS and a large asymptotically AdS black hole jT22j . According to the AdS/CFT corre- 
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sponce such a semiclassical phase transition between asymptoticahy AdS spacetimes has 
a corresponding "confinement /deconfinement" phase transition in the free energy of the 
dual field theory, and that is indeed what Witten found [ 123j . 



2.2.2 Local stability 

Local stability is concerned with whether a certain phase in equilibrium is a local maximum 
of the total entropy. This depends on how the system responds to small fluctuations of 
its thermodynamic parameters. The phase is locally stable if the response counteracts the 
effect of the fluctuation leading back to equilibrium. The standard analysis is based on 
linear response functions, such as specific heats. 

Let us consider the general condition of local thermodynamic stability. As happened 
for global stability, it is a consequence of the second law of thermodynamics. Take a 
black hole in equilibrium with a reservoir at temperature T and angular velocities rij. 



as described by the grand- canonical ensemble. The inequality (2.8) gives the preferred 
evolution of the system. The important point is that, if a phase is locally stable, then 
small fluctuations must be entropically suppressed. That is, we change the sign of the 



middle inequality (2.8) so that small fluctuations {AM, AJj} obey 



, „ AM - QiAJi 

AS — ^ ^<0. (2.10) 



The fluctuations take the black hole temporarily out of equilibrium with the reservoir. For 
small fluctuations AM and A Jj, the change in the entropy of the black hole is 

dS 1 d'^S 

AS = — Ax^ + - Ax^Axy + . . . , x^ = (M,Ji). (2.11) 

Notice that, although the thermodynamic fluctuations are off-equilibrium, the coefficients 
in the expansion above are still determined by the equilibrium phase, because one may 
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take = Xfj^iT, Qi). Using the first law of tliermodynamics, 

1 . n 



dS = -dM-^ dJi , (2.12) 



tlie inequality (2.10) becomes 



AM - QiAJi 1 d'^S 



AS — ^ - = -——-Ax^Ax, + ...<0. (2.13) 

T 2 ox,,ox^ 



For small arbitrary fluctuations Ax^, this is true if and only if 
d^S{x^) 



dxfj,dxu 



Xfj, = (M, Jj) , is positive definite. (2-14) 



This is the condition for local thermodynamic stability. If it is not satisfied, the black hole 
will drop out of its equilibrium phase. 



We can express these conditions in several different ways. In Appendix 2. A 



we 



prove that each of the following statements (the list is not exhaustive) is equivalent to 



condition (2.14) 



• — - — , = {T,Qi) , is positive definite. (2-15) 

• ^J^^ , Xfj, = (5, Ji) , is positive definite. (2.16) 

ox ^oxi/ 

d'^W(yx) 

• , Vfi = (/3, -l3^i) , IS positive definite. (2-17) 
We introduced here /3 = 1 /T and the Legendre transform of the entropy, 

W = S - I3M + fS^iJi = -I3G . (2.18) 



The formulation (2.15) is more convenient to relate the stability condition to the usual 
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linear response functions, like the specific heat. Explicitly, 
where 



(2.19) 



Cn = T(^) (2.20) 



dT J 



n 



is the specific heat at constant angular velocities (all fixed). The isothermal differential 
moment of inertia tensor is 

There is also the vector 



where the second equality, given by the symmetry of the Hessian matrix, corresponds to 
a Maxwell relation. The specific heat at constant angular momenta Cj, defined as 

C.=t(^)^. (2.23) 

satisfies the identity: 

Cj = Cn-T{e-%Vir]j, (2.24) 

which follows simply from the chain and cyclic rules for partial derivatives. We show in 



Appendix |2. A that the usual stability conditions, 

Cj > and eij is positive definite, (2.25) 



amount to the statement (2.15). The identity (2.24) then implies that Cq > Cj. 



Let us also obtain these conditions starting from the partition function. We can 
write the grand-canonical partition function Z = e~^'^ as 

Z{y^) = j {Ildx,) p{x„) e-y^^^ , (2.26) 
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where p{xa-) is the density of states with mass or angular momenta x^- The mean value 
of these variables in the ensemble is given by 



dlnZ 



(2.27) 



If the partition function is well-defined, the statistical variance of a linear combination of 
the variables must be positive: 



{{C^X^ - {CyXu)f) = c^Ci,{{x^Xu) - {x^){xy)) > , 



(2.28) 



where are real coefficients, one of which, at least, is non-zero. This clearly implies that 
the matrix 

{x,x^} - {X,} (X.) - - - (2.29) 

is positive definite, which is simply the statement (2.17). We will introduce in Chapter [s] 
the gravitational partition function of Euclidean quantum gravity, and confirm that it is 
ill-defined when local thermodynamic stability fails. 

Consider now that the fluctuations are restricted to AJj = 0. The conditions are 
relaxed so that Cj > is enough to ensure stability for exchanges AM ^ with the 



reservoir. Notice, for instance, that the matrices in (2.14) and (2.16) give the appropriate 
condition if restricted to the {0, 0} component: 



(2.30) 



The canonical partition function being well-defined requires, analogously to (2.29), 



a^(-/3F(T,J,)) 



T^Cj > 0. 



(2.31) 



What about the local thermodynamic stability of a black hole as an isolated system, 
i.e. in the micro-canonical ensemble? The common derivation of the stability condition 
does not make reference to a reservoir, which we assumed above, but only to the internal 
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equilibrium of the system, where there is no net variation of mass and angular momentum. 
However, those systems are extensive. Suppose that a system, which we take to be in 
internal equilibrium in the phase S{M,Ji), is extensive. Then it can be divided into 
two separate subsystems with equivalent properties, each with entropy S{M/2, Ji/2) = 
S{M, Ji)/2. The suppression of fluctuations within the total system requires 

S{M + AM, Ji + AJi) + S{M - AM, Ji - AJi) < 2 S{M, Ji) , (2.32) 



which for small fluctuations implies (2.14), without reference to a reservoir. 



This fails for gravitational systems because gravity is a long-range interaction. We 
cannot assign a mass to different parts of the spacetime since the energy/mass is non-local. 
It is well known that the Schwarzschild solution is classically stable (at the linear level), 
and yet its specific heat is negativ^ C = — SvrM^ since T = l/SvrM. In fact, we will show 



in Chapter 6l based on 68 , that all asymptotically flat vacuum black holes are unstable 



according to the condition of local thermodynamic stability. 

For extensive systems, the thermodynamic limit gives a direct connection between 
classical and thermodynamic stabilityj^ This clearly fails for black holes, as shown by the 
Schwarzschild case. However, there is also a well-known relation, the Gubser-Mitra conjec- 
ture 125 126 , which states that hlack branes with a non-compact translational symmetry 
are classically stable if and only if they are locally thermodynamically stable. The reason 
for this is discussed in Chapter [6] for vacuum black branes of the type ( |1.11[ ). However, 
we can readily notice that these branes are extensive along the i* directions, since we may 
take a partition of the spacetime by periodically identifying each of the z coordinates, 
z = z + L. The thermodynamic limit in the z directions is exact because U / L^ (where U 
can be the entropy, the mass, the angular momenta or the charge) is independent of L. 



^Negative specific heats are common in self-gravitating systems 124 , e.g. a virialised self-gravitating 
ideal gas with Np particles in four dimensions has energy E = i?kin + i^pot = — i^kin = — '^NpT. 

*This is the gas limit U /V — >■ u{Np/V,T) for large volume V and number of particles Np, where U is 
an extensive quantity, say the energy. The statistical variance in an ensemble becomes negligible in the 
thermodynamic limit making the ensembles equivalent, e.g. for the energy in the canonical or the grand- 
canonical ensembles; {{E^) — {E)^)/{E)^ — > 0. A quantum statistical system then behaves classically for 
quasi-stationary long-wavelength processes. 
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How can the internal equilibrium conditions be derived for black holes? The works 



[35 72 , on which Chapters |6j^ are based, give a partial answer to this question using 
the gravitational partition function of Euclidean quantum gravity. New pathologies in the 
partition function correspond to instabilities which are also classical. They correspond 



to the failure of (2.29) being positive definite, not for fluctuations of the mass and angu- 
lar momenta, but for additional degrees of freedom that are not captured by the usual 
thermodynamic description. 



2. A Appendix: Equivalent statements of local stability 



In this Appendix, we will connect some different ways of stating the condition for local 



thermodynamic stability. Let us start by showing the equivalence of the statements (2.14) 



(2.15), (2.16) and (2.17). We denote 



X/, = (M, Ji), yf, = {T,ni), Xf, = {S,Ji), = {l3,-(3Qi) ^ 



(2.33) 



The first law of thermodynamics (2.6) can be written as 



dS = yfj.dXfj, , or dG = -x^dy^ 



(2.34) 



The Hessian matrices are 



d'^S{xx) djju dy^ dy^ 



32/ 



d'Gjyx) 
dViidyu 



which implies that 



dx^ dx^ dya 
dxy dxa dxu 
dy^ dy^ dx^ 



(5xy • dyy)^u , 

-((^xy)^"^ • 5xx)^j, , 



(2.35) 



(2.36) 
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Since 







-13 (9xx)i.^ , 



(2.37) 



we find that 



S = TU^-G-^-U. 



(2.38) 



Now, the matrix (— S) is positive definite - therefore satisfying the statement (2.14) - if 
and only if 

Vw/O v^-{-S)-v>^. (2.39) 



Prom Eq. (2.38) 



. (_s) .v = T (W)T • (-G)-i • (W) 



(2.40) 



Since |det(U)| = l/T, we conclude that (2.39) is equivalent to 



(2.41) 



so that (— G) is positive definite. Its inverse must also be positive definite, which is 



the statement (2.15). We can easily see that the statements (2.16) and (2.17) are also 



equivalent: since 



dM = y^dx^ , dW = -x^dy^ , 



(2.42) 



we have 



(M) 



d^M{xx) dy^ 



(W) 



d^W(yx) 
dy^idyu 



dxf^dxy dx 



dxy 
dy,i 



(2.43) 



from (I2I35I). 



Let us now express these statements in terms of the usual linear response functions. 
Recall the standard linear algebra result: 
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Theorem: Let L be a symmetric real matrix denoted by 



-'fj.i/ 



A Bi 



(2.44) 



where the matrix Cij is invertible. Then (i) the following determinant identity holds: 



det(L) = {A- {C-')kiBkBi) det(Q,) ; 



(2.45) 



and (ii) L is positive definite if and only if 



A — {C )kiBkBi > and Cij is positive definite . 



(2.46) 



Proof: Notice that L = V"^ L V, where 



A - {C-^)kiBuBi 



and V, 



{C ^)ikBk Si 



(2.47) 



The statement (i) holds since det{V) = 1, and the right-hand side of (2.45) is det(L). 



The statement (ii) follows from steps analogous to Eqs. (2.39)-(2.41 ), showing that L is 
positive definite if and only if L is positive definite. 



If we apply this theorem to the matrix (2.19), for which 



Cij = e^j and A - {C-^)kiBkBi = P[C^- {e-^)ki Vk m) =^Cj, (2.48) 



we see that the conditions (2.25) are equivalent to (2.15) 



To conclude, notice that these conditions can also be obtained using the Helmholtz 
free energy -F(T, Jj). Local stability holds if F is convex with respect to the temperature 
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concave with respect to the angular momenta: 



Part II 

Black hole partition functions 
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Chapter 3 



Thermodynamic negative modes 



The thermodynamic description of black holes reviewed earlier holds at the semiclassical 
level. Quantum effects are incorporated by using thermal quantum field theory on a 
classical spacetime background. Therefore, a semiclassical effective theory of quantum 
gravity - and not necessarily a complete theory - should be able to identify the low energy 
degrees of freedom of thermodynamic instabilities. This is indeed the case for Euclidean 
quantum gravity. 

We saw in the last Chapter that a local thermodynamic instability must be signalled 
by a pathology in the partition function. For Euclidean quantum gravity, the pathology 
consists of a divergent mode in the one-loop quantum corrections of the gravitational 
partition function. It is commonly referred to as a negative mode, since it has the wrong 
sign for convergence. 



3.1 Introduction 



Shortly after the advent of black hole thermodynamics, the Euclidean path integral meth- 
ods of quantum field theory at finite temperature were extended to semiclassical quantum 
gravity |127 , 128 . Gibbons and Hawking jl29| proposed a construction for the partition 
functions of black holes. These were defined as path integrals with given boundary con- 
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ditions, which correspond to fixing the temperature, through the periodicity in imaginary 
time, and possibly other quantities such as angular velocities. A Euclideanised black hole 
solution is seen as a saddle point of the path integral - an instanton - its action being 
related to the thermodynamic free energy. 

The semiclassical approach to the path integral allows for more than that. It is 
possible to go beyond the instanton approximation, corresponding to the classical black 
hole, and analyse small thermal (quantum) fluctuations around it. This was made possible 
after a better understanding of the path integral for gravitational perturbations, through 



the works of Gibbons, Hawking and Perry [130 131 . They found that conformal pertur- 
bations of the metric, which always decrease the Euclidean action and seem to render the 
path integral divergent, are an unphysical artifact that can be eliminated by choosing a 
suitable integration contour and by applying a standard gauge-fixing procedure. 

The first application of those methods was the work of Gross, Perry and YafFe 
[132 . They found that the Schwarzschild instanton possesses a non-conformal radial 
negative mode in the path integral for perturbations, which renders the partition func- 
tion ill-defined. This was expected because the Hawking temperature formula for the 
Schwarzschild black hole, T = l/SvrM, corresponds to a negative specific heat, so that 
there is no local thermodynamic stability. The authors further interpreted the instability 
as the possibility of spontaneous nucleation of black holes in hot flat spacej^ The physics 



of black hole nucleation was clarified when York 135 considered the partition function 
with boundary conditions at finite radius. Such a cavity, with fixed temperature on the 
wall, allows for two black hole solutions. The smaller radius solution is unstable and tends 
to the usual Schwarzschild case when the radius is taken to infinity. The larger radius 
solution is stable and its nucleation is thermodynamically allowed, since its free energy is 
inferior to that of the hot flat space in the cavity. 

A non-trivial test of the correspondence between local thermodynamical stability and 



^In de Sitter space, the cosmological horizon gives rise to the Gibbons-Hawking temperature 133 , and 
a black hole in thermal equilibrium must match this temperature. Ginsparg and Perry found that the 
corresponding instanton possesses a negative mode 134 . 
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a well-defined gravitational partition function was given by Prestidge [136 . He analysed 
the Schwarzschild-AdS instanton, soon after the AdS/CFT correspondence was proposed, 
and found numerically that the negative mode disappears when the specific heat of the 
black hole becomes positive. The curvature of AdS simulates the finite cavity in the 
asymptotically flat case. The correspondence was put on a firmer footing when Reall 
[137 showed, for a certain class of black holes, that a negative specific heat implies the 
existence of a negative mode. The proof of the converse result, however, remained elusive. 



Refs. 68 , 72 , on which Chapters pHq are based, clarified this question by showing that 



negative modes could arise which were not connected with the standard conditions for 
local thermodynamic stability. 

Work on Euclidean negative modes of black holes includes also higher-dimensional 
solutions [138 , 139 , the Taub-NUT and Taub-bolt instantons [140 also with cosmological 
constant [141 , and connections to Ricci-flow 142 , 143 . We shall review in Chapter [6] the 
correspondence between the classical stability of black strings or branes and the existence 
of thermodynamic negative modes. 



This Chapter is organised as follows. In Section 3.2 we review the partition function 
of Euclidean quantum gravity, paying special attention to the conformal factor problem 

we connect 



and the relevance of traceless-transverse perturbation modes. In Section 3.3 



the existence of local thermodynamic instabilities to negative modes of the action. 



3.2 The gravitational partition function 

3.2.1 Euclidean quantum gravity 

The path integral of Euclidean quantum gravity, 

Z = I d[g]e-^[^l , (3.1) 
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is a sum over imaginary time manifolds whose boundaries match a prescribed geometry, 
which is the argument of the gravitational partition function. The periodicity of the 
imaginary time coordinate on the boundary is the inverse temperature /3 = 1/T, as is 
common in thermal field theory, but the periodicity will now also include the rotation 
angles, fixing the angular velocities Oj. The action for Euclidean metrics is 

^^9] = -^ I d^x^{R-2A)-^ [ d^-^xJgiD~i)K-h. (3.2) 

The first term is the usual Einstein-Hilbert action and the second is the York-Gibbons- 
Hawking boundary term [129 ,144 , where K is the trace of the extrinsic curvature on dM.. 



This term is required for non-compact manifolds M , such as the ones we will study, in order 
that the boundary condition on dAi is a fixed induced metric, and not fixed derivatives of 
the metric normal to dM. The inclusion of such a boundary term in the action specifies 
the thermodynamic ensemble, in this case the grand-canonical (fixed temperature T and 
angular velocities fij). 

The term Iq can depend only on ^\ the induced metric on dM, and not on the 
bulk metric gab, so that it can be absorbed into the measure of the path integral. However, 
since we are interested in the partition functions of black holes, it is convenient to choose 
it so that / = for the background spacetime that the black hole solution approaches 
asymptotically. For asymptotically flat black holes [129 , the Einstein-Hilbert term is zero 
and the action becomes 

_ Jl /■ d^'-^x^^) {K - Ko) , (3.3) 
JdM ^ 

where Kq is the trace of the extrinsic curvature of the flat spacetime matching the black 
hole metric on the boundary dM at infinity. This subtraction renders the action of the 
black hole finite. For asymptotically AdS black holes 122[|145 , the boundary terms cancel 



when the background subtraction is performed, but the bulk volume integral diverges and 



requires an analogous subtraction that sets the action of AdS space to zero. See 146 - 151 



for the AdS / GET interpretation of Iq as a counterterm in the dual conformal field theory. 
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In fact, that prescription for regularising the gravitational action is preferable since it is not 
possible in general to embed an arbitrary boundary geometry in the reference spacetime, 
as the subtraction method requires. 



It should be emphasised that the gravitational path integral ( |3.1| ) is known to be 
non-renormalisable and is considered here as a low energy (and low-curvature) effective 



theory. A different issue is that the action (3.2) can be made arbitrarily negative so that 



the path integral appears to be always divergent. As we shall see, this problem can be 
addressed in the semiclassical approximation, where the path integral is dealt with by 
saddle-point methods. We consider a saddle-point gab, i-e. a non-singular solution of the 
equations of motion, 

2A , , 

Rab = ^ _ ^ gab , (3.4) 

usually referred to as a gravitational instanton. While the boundary conditions in the par- 
tition function may admit the existence of several instantons, which are different phases 
of the system, we are here insterested here in the local stability of a single black hole 
phase. A black hole instanton is defined by the analytic continuation t = —it of the 
Lorentzian solution. Regularity at the bolt (instanton horizon) requires the periodic iden- 
tifications of imaginary time and rotation angles: {T,(pi) ~ {T,(pi + 2Tr) ~ {t + P, (f)i—iQif3). 
Such an identification of the rotation angles, which fixes the angular velocities Oj, leads 
to an instanton metric for rotating black holes which is not Euclidean, having complex 
components. 

To go beyond the leading order instanton contribution, we treat as a quantum field 
the normalisable perturbations hab about the saddle-point, gab — >■ gab + hab- This leads to 
a perturbative expansion of the action, 

I[g + h]=I[g] + l2[h;g]+0{h^). (3.5) 



The first order action Ii vanishes since gab obeys the equations of motion, while the second 
order action I2, which gives the one-loop correction, is the action for the quantum field hab 
on the background geometry gab- The effective field theory requires that the background 
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geometry gab has a curvature far below the Planck scale. The partition function is 



^1— loop — 6 



(3.6) 



where (G.F.) denotes all contributions induced by fixing the gauge in the path integral. 
The second order action is given by 



1 



'\h-Gh + ^{6hf 



where • denotes the metric contraction of tensors. We have define43 



(3.7) 



hab — hab — ^Sabh^ 



(3.8) 



and 



{Gh)ab = -VVchab - ^Ra'ifh.d 



We have also introduced the following operations on tensors T: 



(3.9) 



b...c 



ab...c 1 



{aT)ab...c = ^{aTb...c) ) 



following [131 . 



(3.10a) 
(3.10b) 



3.2.2 The conformal factor problem and physical instabilities 

We mentioned that the instantons of rotating black holes are not real Euclidean geometries. 
Although the action I[g] is real, difficulties arise when dealing with quantum corrections. 
We will ignore these subtleties for now and leave them for the next Chapter, where we 
consider Kerr-AdS black holes. In the following, we assume the instanton to be real 
Euclidean. 



The operator G is related to the Lichnerowicz operator Al defined in (2.21 byG^AL — 4A/ {D — 2) 
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Consider the action for a conformal normalisable perturbation hab = fgab, 

{D-l){D-2) r 



hWg;g] = ^ / d^'x^id^y. (3.11) 



This action is always negative, so that the path integral in (3.6) seems irremediably ill- 



defined. This is the perturbative manifestation of the fact that, if the path integral is 



taken without further prescription, there are geometries for which the action (3.2) can be 
made arbitrarily negative. 

This challenge, known as the conformal factor problem, was addressed by Gibbons, 



Hawking and Perry 130 , who proposed that the integration contour over the conformal 



direction in the space of metrics should be imaginary rather than realj^ We follow here 



the more detailed procedure of 131 , straightforwardly extended to higher dimensions. 
We will decompose the second order action, applying a standard gauge fixing procedure, 
and show that the conformal divergent modes are unphysical and do not contribute to the 
one-loop partition function. 

The second order action l2[h\ g] is invariant for the diffeomorphism transformations 
hab hab + VaVb + VfeK = {h + 2aV)ab ■ (3.12) 
Following the Feynman-DeWitt-Faddeev-Popov gauge fixing method, 

(G.F.) = (det C) 5{Ca[h] - Wa) ■ (3.13) 
We consider the linear class of gauges 

Cb[h] = V" (hab - -p9abh')j > (3-14) 
where (3 is an arbitrary constant, so that the Fadeev-Popov determinant (det C) is given 



^This prescription is very suggestive because it corresponds to making the timelike direction in the 
Wheeler-DeWitt metric into a spacelike one, becoming a positive definite metric on the space of spacetime 
metrics. It is thus analogous to the imaginary time prescription in common Euclidean path integrals. 
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by the spectrum of the operator 



{CV)a = -V'VbVa - RabV + ^ " 1 j ^a^bV . (3.15) 

To study the spectrum, let us consider the Hodge-de Rham decompositiorj^ of the gauge 
vector V into harmonic (H), exact (E) and coexact (C) parts, 

V = Vn + VE + Vc. (3.16) 

This induces a decomposition of the action of C, which we denote by Ch for harmonic 
vectors, Ce for exact vectors and Cc for coexact vectors. 

The harmonic part satisfies dVn = and 6Vu = 0. We can check that 

CVH = -^^Vn. (3.17) 

The spectrum is positive for A < and zero for A = 0, with multiphcity given by the 
number of hnearly independent harmonic vector fields. For A > 0, the background solution 



satisfying (3.4) does not allow for harmonic vector fields if assumed to be compact and 



orientable [152 . Thus, the spectrum of Ch is never negative. 



The exact part is such that Ve = dx, where x is a scalar. We can show that 



spec Ce = spec ( 2 



1 \ 2A 

1 □ 



13 J D 



(3.18) 



where the operator on the right-hand side acts on scalars, and □ is the Laplacian. For A < 
0, the operator is positive for /3 > 1, being positive semi-definite for A = 0. For A > 0, the 
Lichnerowicz-Obata theorem tells us that the spectrum of the Laplacian on a compact and 



orientable manifold satisfying (3.4) is bounded from above by —2DA{D — 1)^^(D — 2)"^, 
the saturation of the bound corresponding to the sphere |152| . This implies that, for 
A > 0, the spectrum of Ce is positive for /3 > D. 

^Notice that, while the instanton background is not compact for asymptotically flat/AdS black holes, 
the off-shell perturbations hat in the path integral should be normalisable. 
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The coexact part is such that 6Vc = 0. Hence 

CVc = 26aVc (3.19) 

and the spectrum of Cc can be shown to be positive semi-definite, 

J d'^x^ [Vc ■ CVc] = 2 j [aVc • aVc] > , (3.20) 

with equality for coexact Kilhng vectors. 

The Faddeev-Popov determinant contribution to the partition function is then 

detC ~ (detC'E)(detC'c) , (3.21) 

the tilde denoting that the zero modes have been projected out. The harmonic contribution 
is not explicitly considered because, if it exists (A < 0), it is a positive factor dependent 
only on A and on the dimension of the space of harmonic vector fields, as mentioned 
above; it will not be relevant to our discussion. The contribution from the exact part is 
fundamental since it will cancel the divergent conformal modes. 

In order to make the results independent of the arbitrary vector w in the gauge fixing 



(3.13), the 't Hooft method of averaging over gauges is adopted. The arbitrariness is then 
expressed in terms of a constant 7 introduced by the weighting factor of the averaging. 
The final result will be independent of 7, as required. The unconstrained eff'ective action 
for the perturbations is given by 

If[h;g] = h[h-g] + d"" x^g C^mCaih] = 



1 

16^ 



.Gh + -{i- ^){5hY + ±\^--yh-dh-^[i--] {dh) 



2 



1- 1. x,-.9 7/ 2\ - - -if 2\,'.,n 



(3.22) 

where we denote h = h'^^. 
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We now decompose the quantum field hab into a traceless-transverse (TT) part, a 
traceless-longitudinal (TL) part, built from a vector rj, and a trace (i.e. conformal) part. 



Kb — h^b + + -r^9abh, 



1 

D' 



(3.23) 



with 



TL 

>'ab 



'^{av)ab + J^Qab^-n- 



(3.24) 



The constant /3, unspecified in the gauge condition (3.14), can be chosen so that the 



trace h and the longitudinal vector i] decouple. This requires 

D-2-i-l\~^ 



/3 = 2 1 



D 7 



(3.25) 



The effective action becomes 



^h^'^ ■ Gh^'^ - ai] ■ qAi?7 - -^(5r?n5r/+ 



+ 2(1 — 7) ( 6arj ■ 5arj + j^a6r] ■ a5r] — —5ar] ■ a5r] ) + 



where the operator F is given by 



(3.26) 



D-2 / D-2j-l\ 1. 



(3.27) 



Recalling the choice of /3 (3.25), we find that the operator on the right-hand side of the 



expression (3.18) is given by ADF/{D — 2). The contribution of the ghosts (3.21) can be 
recast as 

det C ~ (det F) (det Cc) ■ (3.28) 



For the vector rj, as we did for V in the ghost part, we perform a Hodge-de Rham 
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decomposition into harmonic, coexact and exact parts, 

r] = m + "nc + VE , (3.29) 
respectively. Using rj-^ = dx, the result for the effective action is then 



\h^^ ■ Gh^^ + -hFh+ 



4 4D 1 

+ jy^l^arjH ■ aryn - laijc • aCcVc - ^ _ 2 ^^^ " , (3.30) 

where we defined the operator 

T>ab = VaVfe - ^gabO . (3.31) 

Notice that the Hodge-de Rham decomposition of rj in harmonic, coexact and exact parts 
gives, for h'^^, a decomposition in 2arj}i, 2ar]c and 2'Dx, respectively. 

Finally, we can evaluate the Gaussian integrals in the partition function. The pre- 



scription of 130 corresponds to taking the imaginary contour for h. The dependence of 



the partition function is 

^i-ioop ~ {detC){detG)-^/^{detF)-^/^{detCc)~^/^{detF)-^/^ 

{detG)-^/\detCc)^/^ . (3.32) 

Again, the tilde on the operators denotes that the zero modes have been projected outj^ 
The Gaussian integrals are subject to (^-function regularisation |131:|. It is understood that 
the spectrum of G here is restricted to traceless-transverse normalisable modes. 

Let us summarise the treatment of the conformal factor problem. Conformal modes 
make the action for perturbations negative. However, these modes do not contribute to 
the path integral. The two factors {det F)~^/'^ arising from the Gaussian integrals in h 
(taken through an imaginary contour) and x cancel with the det F factor arising from the 



'Zero modes can be dealt with by the standard collective coordinates method 153 
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exact part of the Fadeev-Popov determinant. This makes the unphysical character of the 
divergence obvious, at least in perturbation theory. 

The relevant operators are then Cq and G. For a real metric, the operator Cq is 
positive semi-definite, as we have shown above. Once its zero modes are projected out, it 
contributes a positive factor to the final result. Physical instabilities are only possible if 
there are negative eigenvalues of the operator G, 

{Gh^'')ab = \hl^ , (3.33) 
in which case the eigenmodes h'^jf are called negative modes. Notice that, for A 7^ 0, there 



is no gauge ambiguity (3.12) because GaV = 0. 



3.3 Negative modes and local thermodynamic stability 

We saw how to identify a physical pathology of the partition function through the spec- 
trum of the operator G. The contribution of a negative mode is to lower the action I[g + h]. 
Now, we wish to connect local thermodynamic instabilities to the existence of negative 



modes. The discussion follows Ref. [35] . The idea, taken from Ref. 137 , is to consider a 
family of off-shell geometries for the Euclidean path integral for which one can show that 
the Euclidean action decreases in a certain direction if the black hole is locally thermody- 
namically unstable. It then follows that there must exist a Euclidean stationary negative 
mode. 

Let us point out first that the negative modes h'^jf considered before are regular TT 
fields on the instanton background. This is different from saying that the perturbed metric 
9ab + h'^b" is regular, since we are demanding that gab and h^^ are regular separately. 
Regularity of gab requires the association of T and fij to the periodic identification of 
imaginary time and rotation angles. A field on the background gab must respect 
the same periodic identifications in order to be regular. Therefore, the off-shell metric 
9ab + h'^if leaves T and $7j unchanged to first order. 
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Consider a black hole instanton B{x) uniquely specified by parameters x°'. Let 
r(x), il.i{x), M{x), etc. denote the temperature, angular velocities, mass, etc. of this 
solution. We can construct an off-shell generalization B{x,x), specified by parameters 
as follows 154,155]. Assume that the Killing isometries dr and 9^-, where (f)i are the 



rotation angles, of B{x) are preserved. Perform an ADM decomposition of the metric, 
using the imaginary time coordinate r. Take the spatial geometry of B{x, x) to be the 
same as that of B{x). Now choose the lapse function and shift vector so that (i) B{x,x) 
has the same asymptotics as B{x); (ii) B{x,x) is regular everywhere, in particular at the 
bolt, subject to the identifications {T,(j)i) ~ (r, + 2tt) ~ (t + I3{x),cj}i — iCli{x)f3{x)); 
(iii) B{x, x) = B{x). Note that (ii) implies that B{x, x) is a configuration in the Euclidean 
path integral defined for temperature T{x) and angular velocities ^}i{x), for which the 
saddle point is B{x). Calculating the Euclidean action of B{x,x) using the Hamiltonian 
formalism gives 

/(x, x) = fi{x)M{x) - S{x) - /3(x)a(x) Ji(x). (3.34) 

Condition (iii) implies that the geometry with x = x satisfies the equations of motion and 
hence the first derivative of the action with respect to x" must vanish for x° = x". This 
is a consequence of the fact that the black hole satisfies the first law of thermodynamics, 

dM = TdS + ^iJi . (3.35) 

The second derivative of the action, i.e. the Hessian of the action, now reduces to 
d^I \ d^M d'^S d\h \ 



dx^dxl^ J \ dx"dx(^ dx'^dx^ dx'^dx^ 

where the right-hand side is evaluated at x. If the charges M and Jj uniquely parameterize 
the solution, we can choose x" = (M, Jj). We then have 



Therefore, if —Sa/3 fails to be positive definite for a black hole (i.e. if condition (2.14) for 
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local thermodynamic stability fails) , then the Euclidean action decreases in some direction 
and hence the black hole must admit a negative modej^ Given that our off-shell geometries 
are stationary, this negative mode must also be stationary]^ In general, there must be at 
least as many negative modes as there are negative eigenvalues of this Hessian. 

We shall refer to a negative mode whose existence is predicted by this thermody- 
namic argument as a thermodynamic negative mode. We shall see in Chapters |6[|8] that 
there are some negative modes whose existence cannot be predicted in this way. The latter 



modes found in 168,72 are relevant not just for the quantum stability of the black hole 



but also for its classical stability. 



Note that we have not constructed the negative mode exphcitly by this argument: the hnearisation of 
B{x, x) around x — x will give a superposition of eigenfunctions of G. The point is that, since the action 
decreases in some direction, this must involve a negative mode. 

^The operator G commutes with dr so one can work with simultaneous eigenfunctions of these operators. 
Eigenfunctions with different eigenvalues of the latter will be orthogonal. 



Chapter 4 



Kerr-AdS negative mode 



In this Chapter, based on 156 , we wih analyse the semiclassical stabihty of Kerr-AdS black 
holes. In particular, we will see that a stationary and axisymmetric negative mode exists 
only when local thermodynamic stability fails. This test of the gravitational partition 
function is remarkable also because rotating black holes have complex instanton metrics, 
whose subtleties we briefly discuss. 

The lack of symmetry of the Kerr-AdS solution makes the problem much harder 
to solve than the spherically symmetric cases. We address this by applying a spectral 
numerical method to solve linear coupled partial differential equations. The method is 
reviewed in the Appendix at the end of the thesis. 



4.1 Quasi-Euclidean instantons 



The gravitational partition function is defined as the path integral (3.1), a sum over 
geometries with imaginary time r = it. However, while static geometries remain real for 
this analytical continuation, the same does not hold for stationary non-static geometries. 
In the canonical formalism, where 'jij is the metric on a constant time slice, N is the lapse 
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function and N'^ is the shift vector required for rotating spacetimes, we have 

ds"^ = N^dT^ + -fij{dx' - iN'dT){dx^ - INHr). (4.1) 

Regularity at the boh (instanton horizon) requires the periodic identifications of imaginary 
time and rotation angles: (r, 0i) ~ (t, + 27r) ~ (r + /3, — \.Q.ij3). These geometries 
have been called quasi- Euclidean^ This seems to pose a difficulty for the path integral 
because we expect physical quantities to be real. Notice also that the procedure applied 
in the last Chapter to decompose the metric perturbations and deal with the conformal 
factor problem assumed the instanton to be (real) Euclidean. 



Nevertheless, we share the view of 1154, 1551 that quasi-Euclidean instantons pose 



no problem of principle. The instanton action is real and gives the physical free energy 
(divided by the temperature). Notice that, while one might be tempted to take imaginary 



lapse functions, which would make the line element (4.1 ) real, the resulting geometry would 
bear no relation to the Lorentzian black hole, e.g. the bolt radius would be different from 
the event horizon radius and there would be no ergosphere. Furthermore, although this 



can be done for Myers- Perry black holes (1.5) by analytically continuing the rotation 



parameters Oj, it is impossible for black rings because not all conical singularities can be 
removed in the would-be real instanton [157 158 . 



The treatment of quantum corrections about a quasi-Euclidean instanton is more 



subtle. However, Refs. 154 155 show that the action is real not just for the instanton, but 
also for a family of off-shell geometries ~ the one used in Section |3.3| to connect negative 
modes to local thermodynamic stability. The problem of thermodynamic negative modes 
should then be posed in terms of real quantities. 



The procedure leading to the expression (3.32) for the one-loop quantum corrections 
assumed a Euclidean instanton, while we now want to consider the quasi-Euclidean case. 
However, that expression may still hold for an appropriate complex contour of integration 



^Note that, although the metric is complex, the manifold is real. Real coordinates can be defined by 
setting (jii = <j)i + iQ.iT, so the identifications are (r, (j>i) ~ (r, (f>i + 2-k) ~ (r + /3, 
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in the space of perturbations (and ghosts), specified as usual by the steepest descent 
method. This is our assumption. It would be important to construct such a contour 
explicitly. 

The numerical technique applied here differs from the analytical but perhaps sim- 



plistic first approach to the problem in 159 , which could only account for the effect of a 
single direction in the perturbation space. That single direction was provided by an easily 
constructed traceless-transverse perturbation, which kept the second order action real but 
was not an eigenmode. The difficulty with that approach, which may explain the small 
discrepancy in the final result for the Kerr-AdS negative mode, is that it is not clear if 
such a direction lies on the steepest descent path, despite the second order action being 
real. The steepest descent path is here infinite-dimensional, spanned by the normalised 
eigenmodes of G, which we can only determine numerically. 



The family of off-shell geometries related to thermodynamic stability in Section 3.3 



preserves the Killing isometrics dr and d^., where (pi are the rotation angles. Therefore, 



we will focus on stationary and axisymmetric negative modes (3.33) of the Kerr-AdS 
instanton. Since this problem is independent of the time and the rotation angle, which 
are responsible for the instanton metric being non-real, it is equivalent to the Lorentzian 
problem, and it reduces to a set of explicitly real differential equations. 

In the next Section, we review the Kerr-AdS solution and its thermodynamic stabil- 



ity. In Section 4.3, we outline our implementation of the eigenvalue problem for negative 



modes. The results are presented and discussed in Section 4.4 
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4.2 The Kerr-AdS solution 

4.2.1 Instanton geometry 

The instanton geometry is obtained through the analytical continuation of the black hole 
solution: 

2 _ A(r) ,a ._2o,,MV , ^e{0)s{ii^O fr'^ + 



ds = :^^7T7 — (dT -i — sin 9 d(f)] H — r- — — — — - — dcp + iadt 

+ (dr + ro6x sin 9 d9? + S^fr, 9) d9^ , (4.2) 

A(r) 

where £ is the curvature radius of AdS, i'^ = —3/ A, and 

A(r) = (r2 + a2)(l + r2r2)-ror, (4.3a) 

S(r,^)2 = r^ + a^cos^e, (4.3b) 

Ag{9) = l-a^r^cos^e, (4.3c) 

E = l-a^r^. (4.3d) 

The bounds on the parameter space of the instanton are given by the extremality condition 
and by the requirement that \a\ < i, since the limit \a\ — t- £, for which the 3-dimensional 
Einstein universe at infinity rotates at the speed of light, is singular (441. '^^^ avoidance 
of a conical singularity at the bolt, located at r = r^, the largest root of A, requires the 
coordinate identification (r, 0) = (r + /3, (/) + i/3(J7 — a£~'^)). Here, 

= (4 4) 

^ r+(l + a2^-2 + 3r2£-2-a2r;2) ^ > 

is the inverse temperature and 

a(l + r2r2) 

^ = 2^2 (4-5) 

is the angular velocity in a reference frame non-rotating at infinity, the quantity that 
satisfies the first law of thermodynamics [145 . It differs from the angular velocity in the 
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reference frame of (4.2), which is Q' = ^ — ai ^. The mass M, the angular momentum J 



and the Bekenstein-Hawking entropy S are given by 



M=^, J=p, 5=^(r±^. (4.6) 



4.2.2 Thermodynamic stability 

We analyse here the Kerr-AdS black hole not only because of the usual theoretical moti- 
vations (AdS/CFT correspondence) but also because there is a phase transition to local 
thermodynamic stability for large enough black holes in the grand-canonical ensemble. 
Our purpose is to verify that this criterion for stability matches the one given by the 
quantum corrections to the gravitational partition function. 

Local thermodynamic stability in the grand-canonical ensemble requires that both 
the specific heat at constant angular momentum Cj and the isothermal moment of inertia 



e are positive (2.25). These are given by 



27re^ H-i {rl + a^f [r^(3r2 + a^) + e{r\ - a' 



(J ^ ^71 M ^ + -r u, ; p +V^' + -r u ; -r c \i ^ - c "^ 

^ rXl^ - 3r^^2 _ ^2(9^6 + 23r4 £2 + _ ^4(6^4 ^ 13^2^2 + 3^4) _ ^6(^2^ + ^2) 

(4.7) 

and 

7r(r2 + a2)3 [^2 (3^2 ^ ^2) ^ ^2(^^ _ ^2^] 

r+H4[r2£2_3^4 +(^2 +^2)^2] • y -"^l 

An interesting feature of the Kerr-AdS black hole, as opposed to the black ring for instance, 
is that the specific heat at constant angular velocity is sufficient to describe the stability 
in the grand-canonical ensemble. It is positive when both Cj and e are positive, and 
negative when one of them is negative; Cj and e are never simultaneously negative and 
complement the region of instability, as we can see in the left plot of Figure [4Tj Therefore, 
the local thermodynamic stability criterium is 

87r2r+ [r\ + 0?) 

^^^ = "H(l-3rir2 + a2r;^ + a2r2) ^''^^ 
The vanishing of the denominator in the right-hand side identifies the line of critical 
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Figure 4.1: Phase diagrams of the Kerr-AdS black hole. The parameter space is bounded 
above by extremality and by the singular limit a = i. Left plot: Cj < 0, e > and 
Cfi < for the dark blue region; Cj > 0, e < and Cq < for the light blue region; local 
thermodynamical stability for the red region {Cj,e,Cn > 0). Right plot: thermodynamic 
stability in the red region; global thermodynamic instability (r_|_ < i) in the brown, blue 
and green regions; local thermodynamic instability (Cq < 0) in the blue and green regions; 
classical superradiant instability (ill > 1) in the green region. 



stability for each value of £. Across this line, Cq diverges and changes sign. In the static 
case, this occurs for r+ = £/\/3, or tq ~ 0.77^. 



In the right plot of Figure |4.1[ we present a phase diagram including information on 
thermodynamic stability (local and global) and also on classical instability. For Qi > 1, 
in the light blue region, the black hole suffers from the classical superradiant instability 
[117 -119 , the only known instability of these black holes in the absence of matter. In 
the two blue regions, the black hole is locally thermodynamically unstable. In those 
regions and also in the brown region, which all satisfy r-|- < I, the black hole is globally 
thermodynamically unstable because its Gibbs free energy with respect to thermal AdS 
space is positive, 

G = M-TS-nj=^ +^ \ + ^. 4.10 

4 r+ ^2 „ 

In the light red region, the black hole is thermodynamically stable. The transition to this 
region is the Hawking-Page phase transition [122 . The three critical curves meet at the 
point r-(_ = ro/4 = a = 
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4.3 The eigenvalue problem 
4.3.1 Ansatz 

We will look at eigenvalues of the operator G in search of stationary axisymmetric negative 
modes, Gh^^ = Xh'^'^, with A < 0. Our ansatz for the perturbed metric is given by 

2 _ A{r) Juof,_-a-2nJu;, 



ds^ = —ri-^ e*""" (dT-\- sin^ e 
2j^(r, 0) 

+ ^?!^ e^^« idr + ro5x sinO dOf + Y?ir, 9) e^f"^ dO"^ , (4.11) 
A(r) 

where the functions dvo, ^'-'i, ^fJ'O, (^/^ii ^re small perturbations which are functions 
of (r, 9). The imposition of the TT conditions on the linearised perturbation hab, 

h\ = 0, V''Kb = 0, (4.12) 

gives three constraints, which are solved by explicit relations 

T = T{drdiJ,o,drSfii,dr6x, dQ6no,do6iJ,i,dg5x, Sij,o,6ij,i,5x), (4.13) 

where T = {Suq, Si'i, 6u}}. Substituting these three relations in the metric perturbation, we 
get the general expression for respecting the isometries of the instantonj^ Notice that 
this expression includes first derivatives of the functions 6fio, 6x in the components 

t,TT uTT uTT 
"'TT ; ' ' V0 ■ 

Now we look at the eigenvalue problem E = GhJ"^ — A /i^^ = 0. There are six 
equations to be solved, corresponding to the non-zero components E^^, E^^, Err-, 
Er0, Eqq. Since G is a second order differential operator, and /i"^"^ already includes first 
derivatives, the components Err, E^^, Esd, are third order equations. The components 



The inclusion of the terms h-rr, h-re, h^r, h^e in the ansatz (4.111 leads to the conclusion that these 
must vanish when imposing the TT conditions. Hence we are working with the most general ansatz for 
stationary and axisymmetric perturbations. 
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Err, Ere, Egg are second order only, and they imply (by taking their derivatives) that the 
third order components are automatically satisfied, as expected. Hence we have the three 
coupled partial differential equations Err, Erg, Egg for the three unknown functions 5/io, 

One thing we should emphasise is that the subtlety of having a quasi-Euclidean 
metric, as opposed to Euclidean, has vanished now. The equations obtained are real 
equations for real perturbation functions. Unfortunately, it will not be possible to present 
the three equations explicitly here. Even in the asymptotically flat limit ^ — t- oo, they are 
too cumbersome. 

Let us first consider the boundary conditions for the perturbations. The region of 
integration for the differential equations is an infinite strip r+<r<oo,0<^<7r. The 
perturbations must vanish at infinity, r — )■ oo. The boundary conditions on the horizon are 
obtained by considering a regular basis. This procedure is described explicitly in Chapter[7] 
and we will not repeat it here. We find that 6x{r,9) oc (r — r+) near the horizon, while 
5fiQ{r-^-,6) and 5iii{rj^,9) must simply be finite. Regularity at the poles = 0,7r implies 
only that the functions are finite there. 

4.3.2 Implementation 

It is convenient to make rescalings such that only adimensional quantities are involved in 
the problem. 




(4.14b) 



(4.14a) 



and notice that 




(4.14c) 



We will use the coordinate 



X = cos 9 



(4.15) 
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so that the numerical integration wiU be performed on a rectangle < y < Y , —1 < x < 1, 
where y ^ y* must be sufficiently large. Let us rescale the perturbation functions too, 



Sx{r,e) = f^, (4.16b) 
1 — x^ 



6ni{r,e) = — (4.16c) 

y(l - x^) 

so that the boundary conditions discussed before are simply = (i = 1, 2, 3) along the 
edges x = ±1 and y = 0,Y. It is convenient to combine our previous equations rr, r6 and 
66, now in terms of the redefined quantities, into the form 

dlqi{y,x) + ... + X,fiy,x-y„eMy,x) = 0, (4.17) 

which turns out to be possible. Notice that the last term above gives the only dependence of 



the equations on A=k. We are now ready to implement the spectral numerical method 160 
which is briefly described in the Appendix at the end of this thesis. 



4.4 Results 

Let us first consider the Kerr case, i.e. the limit i — t- 00. The results are represented in 



Fig. 4.2 We find that there is a single negative mode for \a\ < ro/2, monotonically in- 



creasing in magnitude with the angular momentum. Surprisingly, our probe perturbation 



method 159 approached the value of the negative eigenvalue now found to within 10%. 
One way to understand the increase in magnitude is to recall the connection between the 
black hole thermodynamic negative mode and the classical Gregory-Lafiamme instabil- 



ity of the respective black string/brane 137], which we shall review in greater detail in 



Chapter [6| The threshold wavenumber k = {k ■ k)^^"^ for the Gregory-Lafiamme insta- 
bility corresponds to the four dimensional stationary solution of GK^^ = —k"^ h'^^ , with 
the appropriate boundary conditions [100 . This is exactly the problem addressed here 
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Figure 4.2: For the Kerr instanton, A* is negative, decreasing monotonically away from 
a = and evaluating to a finite value at extremality \a\ = ro/2. 

if we identify A = —k'^. The fact that we are dealing with a quasi-Euclidean geometry 
rather than a Lorentzian geometry is irrelevant since time plays no role in the solutions 



to the perturbation functions defined in (4.11). The curve in Fig. 4.2 thus implies that 



the Gregory-Laflamme instability of the Kerr string persists up to extremality. The larger 
in magnitude is the negative mode, the smaller is the threshold length scale k^^ for the 
instability. We expect on physical grounds that the centrifugal force caused by the rota- 



tion will favour the instability of ripples along the string (Ref. [110 presents a fluid dual 
analogy) thus decreasing their threshold length scale and explaining the stronger negative 



mode of the black hole. See 68 , 72 , 108 for analogous results in higher dimensions. 



We now turn to the Kerr-AdS case. If we look at Fig. 4.3, the agreement between 



the existence of a negative mode and the condition (4.9) is striking. We flnd a negative 



mode only when Cq is negative. Unfortunately the numerical method does not allow us to 
safely zoom in the line of critical stability. Nevertheless, it is clear that the gravitational 
partition function reproduces the thermodynamics of the system beyond the instanton 
approximation, even for this non-static black hole with a quasi-Euclidean instanton. 
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Figure 4.3: Phase diagram of the Kerr-AdS black hole. The points represent the parameter 
region where wc find a negative mode and the line represents the change of sign of Cq, 
which is negative in the Kerr limit ^ — >■ oo. 
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Chapter 5 

Reissner-Nordstrom negative 
mode 



In this Chapter, based on 161 , we analyse the problem of thermodynamic negative modes 
of the Reissner-Nordstrom black hole in four dimensions. We find analytically that a 
negative mode disappears when the specific heat at constant charge becomes positive. 
The sector of pertm'bations analysed here is included in the canonical partition function 
of the magnetically charged black hole. The result obeys the usual rule that the partition 
function is only well-defined when there is local thermodynamical equilibrium. 

We emphasise the challenge of quantising Einstein-Maxwell theory, even as a low 
energy effective theory, because the conformal factor problem is more intricate than in 



the vacuum case, which was discussed in Section 3.2.2 We circumvent this difficulty by 
considering a dimensional reduction from five to four dimensions. The method allows us 
to decouple the divergent gauge volume and treat the metric perturbations sector in a 
gauge-invariant way. 
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5.1 Introduction 



The inclusion of matter presents considerable challenges to the study of black hole negative 

, re- 



modes. This is because the standard decomposition of the perturbations 130 ,131 



viewed in Section 3.2.2, which singles out the unphysical conformal modes, was performed 
for the Einstein theory only (with or without cosmological constant). In particular, the 
presence of a background gauge field invalidates the procedure even for metric perturba- 
tions only, i.e. when the gauge field perturbations can be consistently set to zero. 

One way out to address this problem is suggested by the alternative procedure of 



Gratton, Lewis and Turok 162,163 for cosmological instantons. They look for well- 
behaved gauge-invariant perturbations which allow for a decoupling in the action between 



the relevant modes and the unphysical ones, which are integrated out. See Kol 164 for a 
general and systematic discussion of the same basic idea, which he calls "power of the ac- 



tion," following its initial application to the Schwarzschild black hole and black string 138 



This corresponds to the procedure we adopt in this work. There is a further complication 
though. The problem of radial perturbations requires two gauge conditions (as the trace- 
less and transverse conditions in pure gravity) for the metric perturbations. The method is 
inconsistent because the radial ansatz for the perturbation of the action only has a radial 
diffeomorphism invariance (one gauge choice), and including time does not allow for the 
construction of gauge- invariant quantities. The work of Kol |138| on the Schwarzschild 
black hole gives a solution to this problem. Considering a lift to one higher dimension, 
along which there is translational invariance, it is possible to construct gauge-invariant 
quantities, i.e. perturbations which are invariant for infinitesimal diffeomorphisms along 
the radial direction and along the extra dimension. If the action for the zero modes 
(infinite wavelength) along the extra-dimension reproduces the lower-dimensional action 
for perturbations, then the higher- dimensional action can be decomposed, and the long 
wavelength limit can be taken when a simple reduced action is available. 

Using this Kaluza-Klein method, we are able to study the four-dimensional magnetic 
Reissner-Nordstrom black hole in a sector of perturbations corresponding to the canonical 
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ensemble, i.e. fixed charge Q. We find tliat one negative mode still exists if the charge is 
small compared with the mass M, as expected, but disappears for \Q\ > ^/3M/2, exactly 
when the specific heat becomes positive. This supports the validity of the canonical 



partition function of Euclidean quantum gravity 129 even when gauge fields are present. 
One might be worried by the fact that such a theory is not renormalisable at one 



loop, which is indeed the case of Einstein-Maxwell theory 1165 , so that we cannot compute 



quantum corrections. However, we can take an effective field theory approach [166 , where 
Einstein-Maxwell theory is regarded as the low energy limit of an underlying fundamental 
theory. The effective theory is valid up to a cut-off scale, after which ultraviolet completion 
effects become important. As long as the energies of the fields involved are nowhere 
near that scale, the perturbative quantisation of the non-renormalisable effective theory is 
meaningful. It may help to recall that Einstein-Maxwell theory corresponds to the bosonic 
sector of four- dimensional J\f = 2 supergravity, which can be embedded into string theory. 

We mentioned that Einstein-Maxwell theory had not been considered for the per- 
turbative path integral quantisation on black hole backgrounds. Miyamoto and Ku- 



doh 167,168 have analysed the classical stability of magnetically charged branes and 
verified that they are stable when there is local thermodynamic stability. This result is 



related to the Gubser-Mitra conjecture mentioned in Sections 2.2.2 and 6.1.1 The issue 
we address here is whether the partition function at one loop, defined as a saddle point 
approximation to a Euclidean path integral, conforms to the usual thermodynamic stabil- 
ity criterion. It is also worth mentioning that Einstein-Maxwell theory in four dimensions 
cannot be the result of dimensionally reducing Einstein-Maxwell theory in five dimensions, 
since there is then an extra scalar field. The fact that the Reissner-Nordstrom black hole 
is not straightforwardly related to a string will force us to be careful in our Kaluza-Klein 
action method, so that we make sure the correct quantum theory is obtained. 



This Chapter is organised as follows. In Section 5.2, we discuss the problems with 
quantising Einstein-Maxwell theory, stressing the differences with the pure Einstein case, 
and recall the relation between the boundary conditions of the path integral and the cor- 
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responding thermodynamic ensemble. In Section 5.3, we explain the Kaluza-KIein action 



method to analyse the second-order action. In Section 5.4, we describe the application 
of the method to the magnetic Reissner-Nordstrom black hole. We start by presenting 
an appropriate lift to five dimensions. We then construct the gauge- invariant quantities 
and obtain the reduced action. We find that the action possesses a negative mode when 



the specific heat at constant charge is negative. Finally, in Section 5.5, we present the 
conclusions. 



5.2 The Einstein-Maxwell path integral 

We wish to consider gravity coupled to electromagnetism in asymptotically flat spacetimes. 
The partition function is given by the path integral 



Z 



I d[g]d[A]e-'^a^''\ (5.1) 



constructed from the Euclidean action 

I[g, A] = - [ d'^x^iR - FabF'^') - 2 / d^-^xVh{K - Ko), (5.2) 
Jm JdM 

twhere the first term corresponds to the usual Einstein-Hilbert action and the second term 
is the Maxwell action. As before, the third term is the Gibbons-Hawking- York boundary 
term 129, 144 , required if the configurations summed over have a prescribed induced 
metric on dM. K represents the trace of the extrinsic curvature on dM and Kq is the 
trace of the extrinsic curvature of flat spacetime, matching the black hole metric at infinity, 
necessary to render the on-shell action of asymptotically fiat black hole solutions finite. 
Recall that Fab are the components of the Maxwell field strength, the exterior derivative 
F = dA of the gauge 1-form potential A. 



As discussed in Section 3.2, the purely gravitational action can be made arbitrarily 
negative for conformal transformations that obey the boundary conditions of the path 
integral, i.e. by geometries included in the sum. In the absence of electromagnetic sources, 
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this apparent divergence of the functional integration, called the conformal factor problem, 



is circumvented by choosing an appropriate complex integration contour 130 , at least at 



the one loop level of the semiclassical quantisation. Furthermore, it was later shown in [131 



that an orthogonal decomposition of the metric perturbations into a trace, a longitudinal- 
traceless and a transverse-traceless part, complemented by an appropriate gauge choice to 
deal with the diffeomorphism invariance of the action, leads to the complete decoupling of 
these components in the second-order action. This procedure was detailed in Section [3. 2. 2 



It is a key requirement that the gauge choice kills the interaction between the trace and 
the longitudinal traceless parts. The final result confirms the prescription of [ 130] , and 
the decoupled trace can be integrated by choosing a suitable complex contour. Moreover, 
the scalar parts of the partition function, comprising the contributions from the tracelike 
(conformal) perturbations and from the scalar parts of the vector modes, both in the 
longitudinal traceless part of the metric and in the ghost vectors, cancel. This shows that 
the apparent non-positivity of the action for perturbations, and the consequent divergence 
of the path integral, are fixed by projecting out this contribution. We would like to address 
this problem in the presence of electromagnetism. 



5.2.1 The second-order action 

In this section, we will perturb the action about a saddle point {g,A), which we define 
here as a non-singular solution of the equations of motion 

Rab-^R gab = 2 (^F, ''Fbc - \ gab PP^Fp^j (5.3a) 

and 

VaF''^ = . (5.3b) 
Small perturbations {hab, cib) about this solution, 



gab Qab + hab and Ab A + flfe, 



(5.4) 
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are treated as quantum fields living on the saddle point background. We then perturb the 
action to second order, 



I[g + h,A + a] = I[g,A] + l2[h,a;g,A] + 0{h^, h^a, ha^^a^^ 



(5.5) 



so that the partition function can be approximated by a saddle point functional integral 



(5.6) 



where (G.F.) denotes all contributions induced by fixing the gauge. The first-order term is 
absent because the instanton solution satisfies the equations of motion. The second-order 
action is given by 



1 



^(abUj' - aattbR"^ + Vaa^Vba^) }, (5.7) 



where we have defined h = g°'^hab- is the Lichnerowicz operator defined in (2.2). This 



action can be checked to be invariant under the following gauge transformations 



hab hab 
fla Oa + VaX 



(5.8) 



(5.9) 



and 

hab hab + VaVfc + VfeVa 
aa^aa + V^Fba 

The first set of transformations corresponds to the U{1) invariance associated with the elec- 
tromagnetic field, and the second to the invariance under diffeomorphisms, conveniently 
mixed with the U{1) symmetry. The gauge ambiguity from both these transformations 
should be fixed in the path integral. 
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Let us point out a property of (5.7) which will be crucial in this work. If we 



focus on spherically symmetric four-dimensional instantons, the metric perturbations can 
be expanded into odd and even perturbations |82j. Substituting this general expansion in 



(5.7), in the background of the magnetic Reissner-Nordstrom solution, the cross term gives 



zero for any value of the vector field perturbation a^. If we further assume that whatever 
gauge choice that makes this problem tractable in four dimensions does not involve the 
mixture between metric perturbations and gauge potential perturbations, then the metric 
sector completely decouples from the gauge potential sector. In the electric Reissner- 
Nordstrom case, this decoupling does not occur, and a complete analysis of the problem 
requires the inclusion of such a term. 

We further decompose the metric perturbations into a traceless and a trace part: 



i>ab + J^9abh. 



(5.10) 



After this expansion, the second-order action can be written as 



2[(Pab,h,aa 



X 0' 



ab 



Ida. 



2D 



(5.11) 



The first term in the third line of (5.11) is the term that makes the study of Einstein- 



Maxwell instantons more involved. It couples the traceless part of the metric with the trace 
part. Hence one cannot simply associate the trace with the unphysical divergent modes, 
as in the pure-gravity conformal factor problem. There may be several ways of tackling 
this. One might try to make a particular gauge choice for the metric perturbations. This 
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is exactly what one does in the purely gravitational case to remove the last term in the 
second line, by choosing a gauge of the form 

Ca[h]=V' (^Kb-^9abh}j , (5.12) 



where /? is a constant to be conveniently fixed 131| . However, we were not able to 



find such a gauge when electromagnetism is introduced without making a shift in the 
gauge potential Oa- This shift results in a new coupling between the trace and the vector 
potential, invalidating once more the pure gravity interpretation of the trace. Furthermore, 
this new choice of gauge requires a new decomposition of the metric perturbations instead 



of the usual longitudinal/transverse decomposition corresponding to the choice (5.12) 



An alternative approach to the problem would be to shift the trace component of the 
metric by a term proportional to 0~^{(j)abF"''^F^^). This choice removes the problematic 
term, but also makes the second-order action dependent on the inverse box operator, 
rendering any further computations unpractical. 



A third method, based on 1381, is to attempt to solve the problem by using Kaluza- 



Klein techniques, and will be the one followed here. 



5.2.2 Boundary conditions 



The partition function is defined as the path integral (5.1) with appropriate boundary 



conditions. These boundary conditions specify the thermodynamic quantities which are 



held fixed in the ensemble. Here, we recall the discussion of Hawking and Ross 169] for 
the Reissner-Nordstrom case. 

As usual, the 3-metric on the boundary dA4 = x fixes the temperature 
T = I3~^ , where /? is the periodicity of imaginary time. The boundary condition on 
the electromagnetic field at infinity typically fixes either the charge Q or the potential 
^ = Q/rj^, as we shall discuss. Imposing a periodicity /3, the leading order approximation 
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to the path integral is the Reissner-Nordstrom instanton, 

ds" = f{r)dT^ + TTT + r^dn^ ' (5-13) 

where /(r) = 1 - 2M/r + Q'^/r'^, dn'^ = dO"^ + sin^ 6i#2, and M and Q are the black hole 
mass and charge, respectively. Since an instanton solution is required to be non-singular, 
the mass and the charge are fixed by the boundary data and by the condition of regularity 



at r = r+ = M + \/ — Q"^, the location of the outer horizon in the Lorentzian solution 
and of the "bolt" in the Euclidean solution. That condition is the formula for the Hawking 
temperature: 

However, the action of the instanton depends on the boundary terms which make 
the fixing of quantities on dAi consistent. Suppose one fixes the potential j4 at a surface of 
very large r = R, where R will be taken to infinity. In the magnetic case, this corresponds 
to specifying the charge 

which is determined by integrating the magnetic field strength over S"^, which in 
turn is determined by A on the boundary alone. But, in the electric case, the charge is 
computed using the dual of the electric field strength, 



An 



Q=±I^^^F, (5.16) 



i.e. it requires fixing F^r (the i is due to the use of imaginary time). The canonical 
ensemble, for which the charge Q is fixed, includes the configurations for which derivatives 
of A normal to the boundary are fixed. The thermodynamic quantity associated with 
specifying A on dA4, in the electric case, is the electric potential at infinity <I> = —iAj- = 
Q/r^ (the gauge choice A = —i{Q/r — ^)dT ensures that A is regular on the horizon 
r = rj^). Fixing the potential <1> corresponds to the grand-canonical ensemble. For the 



canonical ensemble, the action (5.2) must include a boundary term appropriate to the 
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variational problem in question, 

-if d^xy/hF^^UaAb, (5.17) 

JdM 

i.e. to the fixing of i<Vr and thus of the electric charge Q. This term ensures that the 
Helmholtz free energies of the electric and magnetic black holes coincide [169j : 

i^Helmh = M-TS = -rinZcanomcal(/3, Q) , (5.18) 

where S = irr'^ is the black hole entropy. 

In the present work, we look at a sector of metric perturbations about the instanton 
such that the metric is fixed at the boundary. We consider the magnetic case, leaving 
the electromagnetic potential A unperturbed, which can be done consistently for 5*0 (3) 
symmetric backgrounds, as mentioned in the previous Section. The sector is thus included 
in the canonical ensemble. 



5.3 The Kaluza-Klein action method 



We described in Section 5.2 the difficulties in quantising Einstein-Maxwell theory. The 



standard decomposition of the perturbations around the instanton 131 does not apply. 
However, a decomposition that explicitly decouples the divergent modes, the conformal 
modes in the case of pure gravity, is essential. We therefore look for a different approach. 



In 138 , Kol addresses the problem of the negative mode of the Schwarzschild black 



hole by looking at the "dynamical" part of the action. This procedure, useful if the prob- 



lem has a single non-homogeneous dimension (the radial one here), was formalised in 164 
Instead of the treatment of Gross, Perry and Yaffe |132| , which looks at the Lichnerow- 
icz operator Al acting on transverse-traceless metric perturbations, an auxiliary extra 
dimension z is added. The extended space of metric perturbations and the dependence 
on the extra dimension allow for the construction of several gauge-invariant quantities. 
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These decouple into two sectors, a "dynamical" part, which is the relevant reduced action, 
and a "non-dynamical" part, which takes away the divergent modes. The action of the 
five-dimensional zero modes, i.e. the A: = modes in a Fourier decomposition along z, 
or at least a sector of it, is the four-dimensional action. Notice that, in four dimensions, 
the radial problem has a single gauge transformation, The appropriate fixing of the 
gauge freedom requires a second condition (as in the traceless and transverse gauge), but 
including time does not allow for the construction of gauge-invariant quantities. This is 
provided by the component if we use the auxiliary extra dimension instead. 

A different way of looking at the extra dimension is to relate the black hole to a black 
string/brane, as we shall discuss in the next Chapterj^ The difficulty in the charged black 
hole case is that the standard traceless-transverse gauge does not decouple the divergent 
modes. However, the Kaluza-Klein action method ("power of the action", as Kol prefers) 
is still available. It requires: 



(i) A "maximally general ansatz" for the perturbation of the fields; i.e. the ansatz must 
reproduce all of the background field equations by variation of the metric, and must 
be closed for the relevant group of gauge transformations {S^r and ^2). 

(ii) The lift must be such that the five-dimensional action for the k = perturbation 
modes is equivalent to the action for the perturbations around the four-dimensional 
black hole instanton. Actually, it suffices that a particular sector of the five-dimen- 
sional action satisfies this. We must then restrict ourselves to that sector of the path 
integral. 



The lift we consider in this work, along a timelike direction z, is a "magnetic string" 
of a theory with electromagnetism and Chern-Simons term. We will need to restrict to a 
sector of the path integral by introducing a Delta functional on the space of perturbations. 
The restriction will ensure that we obtain the four-dimensional action when we look at 
the A: = modes. 

^ A straightforward correspondence of this type will not hold in our case because of the non-trivial factor 
in front of the second term in (5.501. 
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Before we start, let us make two clarifications. First, how do we know that the 
divergent modes are being decoupled from the path integral? The decomposition of the 
path integral into "dynamical" and "non-dynamical" parts leads to a "non-dynamical 
action" simply composed of squares of gauge-invariant quantities. One of them will have 
a minus sign, meaning that a rotation to the imaginary line is needed in order for the 



Gaussian path integral to converge. This is the exact analogy of the prescription of 130 



Second, the particular eigenvalue obtained from this method is not the same, in general, 
as the one obtained from the standard decomposition. The negative eigenvalue in [iSSj is 



quantitatively different from the one in 132 because the decomposition of the action is 
different. But the positivity properties of the action, i.e. whether a negative eigenvalue 
exists or not, must be the same. 



5.4 Magnetic Reissner-Nordstrom analysis 
5.4.1 Lift to five dimensions 



In this Section, we will apply the method described in Section 5.3 to the magnetic Reissner- 
Nordstrom black hole. The first non-trivial step is to find a five-dimensional system of 
gravity, possibly coupled to some fields, that reduces to four-dimensional Einstein-Maxwell 
theory upon a Kaluza-Klein reduction on a circle. This truncation was first discussed 
in [170 in the context of Lorentzian signature. Here we straightforwardly extend it to 



Euclidean signature spacetimes. 

We start with minimal five-dimensional supergravity with the action given by 

j(5) = _ j dx^^g[R- F-'F^f^ - ^—^F,,,,F,,,,A,^) , (5.19) 

where " represents quantities in five dimensions, g = \ det^.gl and F = dA. Our dimen- 
sional reduction ansatz will take the generic form 



dsl = A^dsl + €/\{dz + lAadx^f , 



(5.20) 
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where 7 is a constant and A is related to the exponential of the dilaton field. Here 
e = ±1, according to the signature of the five-dimensional spacetime. The signature of 
the four- dimensional space is chosen to be Euclidean. Also, dz is a Killing vector of the 
five-dimensional spacetime, meaning that A, gab and Aa do not depend on z. 

The reduction of the Maxwell field parallels that of the metric, and in particular we 
expand it as 

A = A + T.{dz + 2A) , (5.21) 

where both A and A only have components in the four-dimensional space, and again do 
not depend on z. Choosing 7 = —1/2, setting A = exp (— 4(/)/\/3) and substituting in the 



action (5.19) yields the following form of the action: 



- ^S^— {Ha,a2 - 4iai5a,S) {Ha,a, " ^Aa,da,T) , (5.22) 

V o J 

where V is the volume of the circle along which we are performing our dimensional reduc- 
tion and H = F + 2T,F. Varying this action with respect to gab, Aa, Aa, (j) and S yields 
the following equations of motion 

Rab = 2da^db(l) + eVE^da^db^ 

+2e-^^^e(^Fach' - ^FmnF'^'') + 2e'^^'^^HacH,'^ - ^HmnH^''^ , (5.23a) 

_ 2 J 4 fra.ia.2a.3b 

eV,(e-2^^*F'^^) + 2V„(Se-7i*/7-'') + Ha,a2^da,^ = , (5.23b) 

v3 V5 

Va{e ^'^F'^^) + — Ha,a2da,^ = 0, (5.23c) 

v3 \/g 

BcP + ^ePabF'^' + -^e-VE'f'HabH-' - -^e^*5'^Sa,S = , (5.23d) 
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V,(e^V'^S) - --=^—{Ha,a,-^Aa,da,^){Ha,a,-4Aa,da,^) 

v3 \/9 

+ A^—Ha,a,Aa,daJ^ " 26" ^'^//'^^F.b = . (5.23e) 

We now search for a consistent truncation in which the two fields 6 and S are set to zero. 



We are then left with two constraints on F and F coming from both Eqs. (5.23d) and 



(5.23e) 



3eF"''Fafe + = and *F + \/3F = , (5.24) 

where -k is the Hodge dual with respect to the four-dimensional geometry. If the four- 
dimensional manifold has Lorentzian signature, the second condition solves the first if e = 
1. However, if the four- dimensional manifold is assumed to have the Euclidean signature, 
one must require e = — 1. We thus from now on choose e = — 1. The equations of motion 



(5.23) then reduce to 



n,—-2(F F,^-—F F"^"-^-\-2(f F,'^-—r 



4 

VaF"^ = VaF"'' = , (5.25) 



subject to the second constraint in (5.24). In order to further simplify the equations of 
motion and to explicitly solve the constraint, we change from F and F to F^^^ and F^^^ 
given by 

FW = i*F-#F and F(2) = ^^ * F + 1f . (5.26) 



The constraint is solved by setting ★F^^) = 0, that is, F^"^^ = 0. Substitution in (5.25) 
yields 

Rab = 2(f£V(i); - ^F^Fd)'"") , V„(F(i)'^^) = 0. (5.27) 



These are precisely the equations of motion that we were searching for, i.e. the Einstein- 
Maxwell equations in four dimensions. Note that both these equations can be deduced 
from the four-dimensional action 

/(^) = -/dx^V^(i?-Fi;W). (5.28) 
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The four-dimensional quantities are related to the five-dimensional quantities via 



dsi = dsi - {dz + 2Aadx''Y , F 



and A 



yl(i) . (5.29) 



We further remark that this truncation is only valid at the level of the equations of motion. 



In fact, directly substituting the ansatz (5.29) in (5.22) gives 



(5.30) 



We can now write the five-dimensional lift of both the magnetic and electric Reissner- 
Nordstrom instantons, respectively as 



dsi = mdr' + ^ + r^dn^ - ( 



f{r) 



dz — —dr^ 
r 



A 



V3 



Q cos 



(5.31) 



and 



dsl = f(r)dT^ + ^ +r^d^}^-(dz-iQ cos edcpY , A = -i^^dT, (5.32) 

/(r) V ; 2r 

where /(r) = 1 - 2M/r + Q'^/r'^, dn'^ = dO"^ + sin^ ed<t)'^, and M and Q are the black hole 
mass and charge, respectively. We are now ready to use the "power of action". 



. \/3Q 



5.4.2 Reducing the action 

The "maximally general ansatz" for the magnetic case can be written as 

dsi = e^^^'-'^^idr - K{r, z)f + e^^^^^'^Ur^ + e^'^^^'^Un' 

-e^P^'^'^ [dz - r(r, z)dT - a(r, z)drf , (5.33) 

and A = -(^/3Q/2) cos6d(f). With this ansatz, it is consistent, in terms of gauge- 
invariance, not to perturb the electromagnetic potential. This conforms to what we have 
seen before, because in the magnetic case the metric perturbations decouple from the 
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vector potential perturbations, and the former have a positive definite action. For the 



electric case, one would have to find a form of perturbing the solution (5.32), maintaining 



all variables involved only dependent on r and z. However, the authors were not able to 
accomplish this. From now on, we will only consider the magnetic case, for which the 



metric is explicitly codimension one. We then perturb the quantities in (5.33) as 



A{r, z) = AQ{r) + a(r, z) , B{r, z) = Bo{r) + b{r, z) 



C{r, z) = Co(r) + c(r, z) , r(r, z) = ro(r) + 7(r, z) 



(5.34) 



where all lower case letters are perturbations, and thus absent in the background solution, 
which is given by 



Ac, 



-So = 2 log/, 



Cnir) = logr 



and 



Pn fr) 



Q 

r 



(5.35) 



At zero-th order, we exactly recover (5.31). We then substitute (5.34) into the action 



(5.19) and expand it to second order, 



= j drdziP^j^j dau^dbu-^ + Qfju^dau-^ + Viju^u^) , (5.36) 

where u = {a, 6, c, a, /3, 7, k}, / G {1, ...,7} and the independent non-vanishing compo- 
nents of the tensors P°'j^j , Qj"'j and Vjj are given in Appendix 5. A We can now expand 
all the fields in Fourier modes, take the limits A; = and /3 = 0, and compare the five- 
dimensional quadratic action with the four-dimensional counterpart. As we predicted in 



(5.30), the two actions are not equal. Instead, they differ by the following term: 



L 



(5) 



2,fc=0,/3=0 



I. 



(4) 



dr 



{r'^dr-fir) + Q[a{r) + b{r) - 2c(r)]}^ 
2^ 



(5.37) 



Note that here we are perturbing the four-dimensional action only in the metric sector, 
and with an ansatz equal to the / = ansatz for the Schwarzschild perturbations used 



in 1132 . We will see later how to deal with this difference in actions. 
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The second-order action (5.36) is invariant under the foUowing gauge transformations 



<5a = e-2^oer^'o + e-2^or2a,6 
6f3 = -9,6 

^7 = e-2s«err^ + 2roa.e. 



6Hi 



-2Ao 



[e.n + ro(9,e. + a,6-2e.A'o)], 



(5.38) 



where ' denotes differentiation with respect to r in zero-th-order quantities. These corre- 
spond to infinitesimal diffeomorphisms along the 1-form ^ = ^rdr + izdz. 



The quadratic action (5.36) can be cast in a different form, in which a and k only 



appear as dza and dzK, if we integrate by parts. This was expected because, in the /c = 
sector, one does not need to set a or k to zero, as one does for /3. Let us then consider 
the following gauge independent quantities 



h + 



c4 
C' 



e-2^«r2/3 



-2Bo 



\ ^0 

e -« f/3r'o - 2ro/3yl'o 



7 



-2ylo 



^0 

2(Ao+Bo) 



°92c + roa,/3 



(5.39) 



2ro/3 



Plugging in the expressions for the u^'s, as a function of the g*'s, in (|5.36l), we obtain a 



quadratic action which only depends on the g*'s. As expected, the action can be entirely 
written in terms of gauge-invariant quantities. We now proceed to integrate the "non- 
dynamical" quantities g^, and g^, which appear with no derivatives in the action. The 
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action for the perturbations can be written as 



r(5) _ .(5) 



^2;\,5 + / drdz 



(5.40) 



where i G {2,3,4}. L does not depend on or and R depends on g^, and their 
derivatives. Exphcitly, 



1 

2r2 2 2^ ■'^ 



and 



L 



3fQ 



fr' -=-fQr 



91 
2 



R 



2r 



r — rf 



r(5) 



^ + 2Qq^ + r{2rfdrq^ - Qdrq^) 



2r2 



We can thus integrate the q^ by constructing the following shifted variables 



(5.41) 



(5.42) 



(5.43) 



i = q'- {L-'fR~^ , 



(5.44) 



in terms of which the action (5.40) can be rewritten as 



r(5) _ A5) 

2 — -'q «i 



drdzR-iL-^fR-.+lND, 



(5.45) 
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where 



Ind = / drdzq L-. 



(5.46) 



The "non-dynamical" g*'s have now decoupled from the and dependent part of the 
action. The new effective action for the variables and q^ is 



drdzRiiL-^f^ R-. 



(5.47) 



This last expression is too cumbersome to be shown. However, we still have freedom to 
rotate and q^ . If one sets 



q^ = aiiq + ai^qf 
q^ = a^iq + 055^^ 



(5.48) 



where 



an 



2^/6/r2 



, «15 



^/6r 



, (5.49) 



then g"' disappears from the effective action, which reduces to the simple expression 



/ 



drdzr 



fidr 



^ + Vq'^ 



(5.50) 



where V is given by 



V 



6M2 



r2(3M-2r)2 



4 



3 V M 



(5.51) 



We are now in a position to explain why the Kaluza-Klein method works in this 
specific case. The path integral for perturbations in five dimensions Z^g) reduces to the 
one in four dimensions Z^2) k = modes if we only integrate over directions in which 



the difference (5.37) between the two actions is zero, i.e. if we tune 7. So, one has a 



particular sector of Z^^-^ , 



a + 6 - 2c) + r'^dr-iW] = Zfi) > 



(5.52) 
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where 6[.] is a Dirac delta functional. Above, Z'^^) denotes the restriction of the four- 
dimensional path integral to the sector of radial metric perturbations (we showed before 
that it is consistent, in the magnetic case, to leave the electromagnetic potential unper- 
turbed). Introducing the gauge-invariant variables gives the following form 



= lim / dimcme- 



^2f 



m 



lim 

k=0 



S[Q{q' + q^) + 



lim / d[q]d[f]d[q^]e-^^«^''^-^^^^'^'^5 



+ 



2r^ 



lim / d\q]e-^-«^ 
k=oJ 



(5.53) 



where each equality holds up to an infinite constant. It is the crucial fact that the final form 
of the action, /cffl^]) does not depend on q^ that makes the two path integrals equivalent 
when we look at the k = sector. The coefficients Ojj were chosen in such a way that the 
argument of the last Dirac delta functional depends on drq^ and not on both q^ and drq^ 
simult aneously. 



5.4.3 The negative mode 

In this section, we will compute the negative mode of the action by directly studying the 



eigenvalue associated with (5.50), when the z dependence drops out (A; = 0). The action 
reduces to 

IeB,k=o = J drr^ [fidrqf + Vq^] . (5.54) 

We first remark that, if \Q\/M > the potential V becomes positive, and as a result 

the action does not have a negative mode beyond this value of the charge. This is in exact 
agreement with the thermodynamic prediction, as we shall see. 



To proceed, we change variables in (5.54), in such a way that the resulting eigen- 



value problem reduces to a one-dimensional Schrodinger equation. A convenient change 
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of variables is 



and 



r 



log 



(5.55) 



where r± = M ± ^/M^ - 

are the locations of the black hole outer and inner horizons, 



respectively, and 7" is a constant chosen for later convenience. The action (5.54) reduces 
to 



IeS,k= 



dy 



{dyuf + Vu' 



where 



V 



r_|_ —r_ 



r2 


r 1 — r 


2 




l + ee T y 





(5.56) 



(5.57) 



and £ = (r_^ — 3r_)/(3r_^ —r_). Here we choose T = {r^ — r_)^/y^, making the potential 
negative definite. This choice is only valid for e > 0, that is, \Q\/M < \/3/2. The 
eigenvalue problem can now be formulated as 



dyU + Vu = Xu. 



(5.58) 



The change of coordinates (5.55) maps r = r , to y = +oo and r = +oc to y = 0. As 



a result, the boundary conditions are now changed to u being regular at y = +oo, and 



integrable near y = 0. Equation (5.58) can be analytically solved for these boundary 



conditions and one finds the unique normalised bound state 

{r , +r_)H 



u 



4T 



corresponding to 



('^+-'^-) (e _ 1) + (e + 1) coth 
2Vl-r?2-l 



{r+-r_)y 



2T 



(5.59) 



(5.60) 



64M2(1 - r?2)2(2/r^+ 1)2' 
where 77 = \Q\/M. As expected, there is only one negative mode, and it disappears for 



7] > (Fig. 5.1). 
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Figure 5.1: Evolution of the negative mode with \Q\/AI. 



5.4.4 Local thermodynamic stability 

Let us recah the thermodynamics of Reissner-Nordstrom black holes. The partition func- 
tions correspond to the different ensembles according to the boundary conditions of the 



path integral, as explained in Section 5.2.2 In the magnetic case, fixing the temperature 
T = by imposing the periodicity in imaginary time, and fixing the electromagnetic 
potential Aa at infinity, which gives the magnetic charge Q, corresponds to the canonical 
ensemble. 

The thermodynamic stability condition for the Reissner-Nordstrom black hole is 
stated in Section 2.2.2 if we substitute — t- <I> = Q/r^ and J ^ Q. In the canonical 



ensemble, local stability requires that the specific heat at constant charge, 

27rrl(r_|_ — r_) 



3r_ 



(5.61) 



is positive. This occurs for ^/3M/2 < \Q\ < M, which is exactly the range we found 
previously for the disappearance of the negative mode in the partition function. 

Had we studied the grand-canonical ensemble, we would expect a negative mode to 
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persist. The stability condition is the positivity of not just Cq but also the isothermal 
permittivity e, really the capacitance here, which is given by 

dQ\ _ T7f _r+(r+-3r_) 



Since Cq and e have opposite signs, the grand-canonical ensemble is always unstable and 
the negative mode is not expected to disappear. 



5.5 Conclusions 

In this Chapter, we have studied the problem of negative modes of the Euclidean section 
of the magnetic Reissner- Nordstrom black hole in four dimensions. Solving this problem 
within four dimensions seems very difficult, as the identification of the unphysical pertur- 
bations which render the partition function divergent becomes considerably more intricate 
than in the vacuum case. 



Following [138 , we devised a method to study this problem by lifting the magnetic 



Reissner-Nordstrom solution to five dimensions. The five-dimensional action is equal to 
the four-dimensional action, up to a quadratic term that can be set to zero by a suitable 
constraint, imposed by a Dirac delta functional, on the five-dimensional path integral. Fur- 
thermore, in five dimensions, the action can be solely written in terms of gauge-invariant 
variables, which in turn can be divided into two decoupled sectors: "non-dynamical" and 
"dynamical". The former is algebraic, and can thus be readily integrated out. The final 
form of the five-dimensional action depends on a single gauge-invariant variable, and the 
study of its negative modes in the long wavelength limit, corresponding to four dimensions, 
is now accessible. 

We found complete agreement between the local stability of the canonical ensemble 
and the existence of the negative mode. We analytically determined the eigenmode as a 
function of the black hole mass and magnetic charge, from which we concluded that the 
negative mode ceased to exist for \Q\/M > \/3/2. This is the range for which the specific 
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heat at constant charge Cq becomes positive. 

The Kaluza-Klein action method is a practical procedure to determine the quantum 
stabihty of gravity coupled to electromagnetism. However, a standard treatment of the 



metric and electromagnetic potential perturbations along the lines of 131 remains most 
desirable. The clarification of the divergent modes problem, analogous to the conformal 
factor problem of pure Einstein theory, would possibly lead to a better understanding of 
perturbative quantum gravity coupled to matter. 



5. A Appendix: Components of the tensors P, Q and V 



The independent components of P'^j j in the expression (5.36) are given by 



pi 1 


_ pi 1 - pi 1 - 
— -^33 — -^35 ~ 
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— ^-'^ 6 7 ' 


pi 1 _ pi 2 




pi 2 
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(5.63) 

where all the other non- vanishing components can be obtained via the symmetry P'^j = 
P^j ""j . The terms that only involve one derivative were written using the auxiliary tensor 
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Qj°'j, whose non-zero components are given by 

Ql 1 = Q5 1 = ~2Q2 1) Ql 3~ Q5 3~ "2^2 3' Ql 5 = Q5 5~ ~'^Q2 5' 

Qb 6~ ~'^Ql 6' ^1 6 ~ Q2 6' ^3 1 ~ —2(52 1' Q3 3~ ~'^Q2 3' ^3 5~ ~^Q2 5' 

Q3 6 — ~Q6 4 — ~2Q2 6' Q(i 4 ~ 6' '^2 7 — ^Ql 7' 

<5i'4 = 2^-r + r/, Q^\ = Qf, = 4r/, Qa'e = 

Qa's = -2^ + 2r + 2r/, Q2'5 = - ^ + ^ + 3r-/, Q2\ = r{l + 3f), 
^3'4 = -2(^ + 2r/), g3\ = 2Q-2^-4g/, Qg2^ = -52!, 
Q5\ = q(i-^), Q^\ = -^+r-rf. 

(5.6 

Finally, the potential Vjj takes the following form: 
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Chapter 6 



Classical stability and 
thermodynamics 



In this Chapter, we will discuss the connection between the local thermodynamic insta- 
bility and the classical stability of black branes and black holes, focusing on vacuum solu- 
tions without cosmological constant. The connection arises from the fact that stationary 
zero-modes of the Euclidean action, which mark the onset of pathologies of the partition 
function, are also classical stationary zero-modes. These may be thermodynamic in origin, 
being predicted by the thermodynamic stability Hessian, in which case they identify the 
threshold of new Gregory-Laflamme-type instabilities of black branes, but not of black 
holes. Or they may correspond to additional degrees of freedom in the partition function, 
not captured by the usual thermodynamic description, in which case they may indeed 
correspond to classical instabilities of black holes. 

When connected to the classical instability of a black brane or a black hole, a 
stationary zero-mode corresponds to the bifurcation to a new family of solutions. When 
connected to the thermodynamic instability of a black hole, the zero-mode corresponds 
simply to a change in the asymptotic charges within the same family of solutions. 



The discussion here, taken mainly from Refs. 68 , 72 , provides the explanation for 
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the numerical results in the following two Chapters. 



6.1 Negative modes and black brane stability 



We start with a review of the Gubser-Mitra conjecture. Then we will show, based on 68 
that this conjecture implies the classical instability of any black brane constructed by 
trivially extending a vacuum black hole along extra dimensions. 

6.1.1 The Gubser-Mitra conjecture 

The relevance of black branes in string theory motivated the study of the stability of these 
objects. Gregory and Laflamme showed that Schwarzschild black branes are classically 
unstable for perturbation modes with a large wavelength along the extended directions 
[100 . As we discussed in Section 2.1.2 the linearised Einstein equations for those modes 



in the traceless-transverse gauge are 

iALh)ab = -k\ab- (6.1) 

Here, k is the wavenumber of perturbations of the type e^^'^hab, and is the Lich- 
nerowicz operator on the black hole background. The instability occurs for modes with 
I A; I < A;*, where the critical wavelength k^, corresponds to the threshold stationary mode. 
The inclusion of charge may improve the stability properties. Recall Chapter [5j where 
we saw how the charge can make a thermodynamic instability disappear. Indeed, a con- 
nection between the classical and thermodynamic stability of black branes was found. It 



is known as the Gubser-Mitra conjecture 1125,126 , which states that black branes with 



a non-compact translational symmetry are classically stable if and only if they are lo- 



cally thermodynamically stable. We alluded to this conjecture in Section 2.2.2, when we 



discussed local thermodynamic stability; see the review [102 and references therein. 



The problem was clarified by Reall [137], who noticed that Eq. (6.1 ) for the station- 



ary and axisymmetric mode marking the onset of the black brane instability. A; = A;*, is 
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the same as Eq. (3.33) for traceless-transverse negative modes of the Euchdean partition 
function of black holes, if we identify A = —k"^. Recall also Section 3.3 where we show 
how negative modes arise from local thermodynamic instabilities. If there is a negative 
eigenvalue of the thermodynamic stability Hessian, there will be a negative mode hab of 
the black hole partition function. Hence, there will also be a threshold for a classical 
Gregory-Laflamme instability of the associated black branes. 

The stationary perturbation mode with k = k^, is also the perturbative signal of 
the bifurcation to a new family of non- uniform black brane solutions [104 105 . Recall in 



particular Figure 2.2 in Section 2.1.2 



6.1.2 Rotating black branes are unstable 

We shall now prove, based on [68], that the condition for local thermodynamic stability 
fails for all asymptotically flat vacuum black holes. According to the discussion above, 
this implies that all black branes of the type ( |1.11[ ) constructed with such a black hole are 
classically unstable. 

The proof goes as follows. Consider the Legendre transform of the entropy, 

W = S-/3M + pn,Ji , (6.2) 
in terms of which the first law is expressed as dW = —MdP + /3ilj(iJj. Recall from 



Section 2.2 expression (2.17), that local thermodynamic stability exists if and only if 



Wai3 = r.-. , Va = (/3, -P^i) , is positive definite. (6.3) 
A sufficient condition for Wa/^ not being positive definite is 
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Now, the Smarr relation, valid for asymptotically flat vacuum black holes, reads 

^^^M = TS -^iJi ^ M = -{D-2)WT. (6.5) 

Hence we have for such black holes, using the first law again. 

Woo = -{D- 3)MT < , (6.6) 

which implies that local thermodynamic stability fails. It then follows that any vacuum 
black hole solution must admit a thermodynamic negative mode, and that any black brane 
based on such a black hole solution must always be classically unstable. This explains our 
results for the Kerr black hole (Chapter [4]) and for Myers- Perry black holes with a single 
spin (Chapter [7]) and equal spins (Chapter [s]). We find that at least one continuous 
negative mode exists for any value of the rotation parameter, so that the black branes are 
unstable. In fact, we find for Myers-Perry black holes, which may exhibit more than one 
negative mode, that this is the most negative. 

Notice that this negative mode must be present for all values of the rotation, and 
thus there is no critical rotation for which it reduces to a zero-mode. Hence there cannot 
be a classical instability of the black hole associated with this particular mode. 



6.2 Zero-modes and black hole stability 



In this Section, we will argue that a stationary zero-mode of the black hole can correspond: 
(i) to a change in the parameters of the solution, if the zero-mode can be predicted by 
the thermodynamic Hessian, or (ii) to the threshold of a classical instability of the black 
hole, not just the black brane. We will also conjecture that classical instabilities associated 
with stationary zero-modes can only appear in a regime which we call ultraspinning. We 
analyse this regime for the Myers-Perry family of solutions. 



6.2. ZERO-MODES AND BLACK HOLE STABILITY 
6.2.1 The ultraspinning regime 
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In Ref. [71], rapidly-rotating Myers-Perry black holes with a single spin in D > 6 were 
conjectured to be unstable for Gregory-Laflamme-type modes. Such black holes have 
quasi-extended event horizons for large values of the angular momentum compared to the 
mass {\J\/M D-3 ^ i)^ acquiring some properties of black branes, which are Gregory- 
Laflamme-unstable. An order of magnitude estimate for the threshold of the black hole 



instability was given in 71 by considering the thermodynamic behaviour. For small 
rotations, the temperature decreases with the rotation for fixed mass, as it happens for 
the Kerr black hole, which is expected to be classically stable. However, after a critical 
value of the rotation, the temperature actually increases, as it happens for black branes. 
Classical stability was conjectured to fail roughly after the critical rotation. 



The critical rotation is actually a thermodynamic zero- mode of the Hessian (6.3) 



beyond which the black hole possesses two thermodynamic instabilities. To understand 



this, recall our discussion in Section 3.3 Since a traceless-transverse negative mode hab is a 
regular tensor on the black hole background Qab, it follows that T and fij are left unchanged 
by hab. A negative mode, which is an off-shell perturbation in the path integral, occurs 
when the thermodynamic Hessian has a negative eigenvalue. However, this argument also 
applies if a negative mode is continuously connected through the variation of the black 
hole parameters to a zero-mode (A = 0), which is a classical perturbation of the black hole. 
This will be the case if one of the eigenvalues of the thermodynamic Hessian changes from 
positive to negative continuously, marking the appearance of a new local thermodynamic 
instability. It is now more convenient to use the thermodynamic Hessian 



which is simply the inverse matrix of Wab (recall Appendix |2.A ). If the zero- mode hab pre- 
serves T and Jlj then, using the first law, it preserves {dS/dM)j = 1/T and {dS / dJijM = 
—Qi/T, that is, it preserves dS/dx". Hence it must correspond to an eigenvector of Sap 



106 



CHAPTER 6. CLASSICAL STABILITY AND THERMODYNAMICS 



with eigenvalue zero: 

= 6{daS) = Sx^dpd^S = S^pdxP. (6.8) 

Notice that the corresponding thermodynamic instabihty occurs only for a certain range 
of the rotation, beyond the zero-mode, while the thermodynamic instability shown to exist 
for all vacuum black holes in the last Section occurs for any value of the rotation. These 
instabilities correspond to distinct negative eigenvalues of the thermodynamic Hessian, and 
hence to distinct Gregory-Laflamme instabilities of the black branes. This is a refinement 
of the Gubser-Mitra conjecture. 



Ref. 72 confirmed numerically this picture and the conjecture of Ref. [71], showing 
that a new negative mode appears at the critical rotation and, more importantly, that 
additional negative modes which are not thermodynamic in origin occur for higher rota- 
tions. These are thresholds for classical instabilities of the black hole, as explicitly verified 
for the equal spin case in [68], on which Chapter [s] is based. Notice that the zero-modes 
which are thermodynamic in origin can be identified by the simpler Hessian matrix 



due to the identity 



det(-5„^) = det(//.,) , (6.10) 



valid for asymptotically flat vacuum black holes; for a proof of this identity, see Ap- 



pendix 6. A It follows that, for a black hole parameterised by (M, Jj), additional negative 
eigenvalues of —Sap correspond precisely to negative eigenvalues of Hij. 

In the Myers-Perry case, for fixed M, the eigenvalues of Hij are all positive for small 
enough angular momenta. However, as some or all of the angular momenta are increased, 
an eigenvalue of Hij may become negative. If we consider the space parameterised by Jj 
(for fixed M), there is some region containing the origin in which Hij is positive definite. 



We define the boundary of this region to be the ultraspinning surface. Following Ref. 72 



we shall say that a given black hole is ultraspinning if it lies outside the ultraspinning 
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surface. From the above arguments, we know that as one crosses the ultraspinning surface, 
the black hole will develop a new negative mode, and the associated black branes will 
develop a new classical instability. This is in addition to the instability already present at 
low angular momenta. Furthermore, on the ultraspinning surface, the new negative mode 
must reduce to a stationary zero- mode that corresponds to a variation of parameters within 
the Myers-Perry family of solutions, since they are identified by the Myers-Perry equation 
of state S{M, Ji). 

Ref. [72] conjectured that classical instabilities whose threshold is a stationary and 
axisymmetric zero-mode occur only for rotations higher than a thermodynamic zero-mode, 
i.e. in the ultraspinning regime. We emphasise that our conjecture gives a necessary 
condition for an instability, not a sufficient one. 

The intuition leading to the conjecture is that modes of lower symmetry are usually 
the most unstable ones. For instance, the original Gregory-Laflamme instability occurs 
for the "s-wave" of the transversal black hole. An additional classical instability will arise 
after a critical value of the rotation, and it will correspond to a "p-wave" of the transversal 
black hole. As the rotation is increased, higher order waves may become unstable. Now, 
if we consider a black hole, instead of a black brane, the "s-wave" and the "p-wave" are 
associated with the asymptotic charges, mass and angular momenta. Therefore they are 
associated to purely thermodynamic instabilities. Higher order waves, on the other hand, 
may become classically unstable as the rotation is increased, starting with the "d-wave" . 
Notice that these waves do not affect the asymptotic chargesj^ The results of Refs. 68 , 72 



described in the next two Chapters, clarify the harmonic structure underlying this problem. 

Recall also that the thresholds of the classical instabilities should be associated with 
bifurcations to different black hole families, highlighting the connection between stability 
and uniqueness. 



^Recall that, in Section 2.2.2 we discussed how the thermodynamic Hessian gave the stabihty with 
respect to mass/angular momenta exchanges with a "reservoir". If the asymptotic charges are unaltered, 
no "reservoir" is invoked, and it is the internal stability of the black hole which is being probed. 
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Silr, 3 



Figure 6.1: Parameter space for MP black holes with D = 6 (left) and D = 7 (right). The 
black hole is labelled by the horizon radius r_|_ and the spin parameters Oj which we take 
to be positive for clarity. For D = 6, the blue curve corresponds to extreme black holes. 
In the red region, both eigenvalues of Hij are positive. In the blue region, corresponding 
to ultraspinning black holes, one eigenvalue is positive and the other is negative. For 
D = 7, the blue surface corresponds to extreme black holes. The ultraspinning surface is 
the red surface near the origin. Inside this surface, Hij is positive definite. The orange 
surface is where another eigenvalue of H^j vanishes. Between the red and orange surfaces, 
two eigenvalues of Hij are positive and one is negative. Between the orange and blue 
surfaces, one eigenvalue of Hij is positive and two are negative. Ultraspinning black holes 
correspond to points between the red and blue surfaces. Note that the "cusp" where the 
red and orange surfaces meet has equal spins. 



6.2.2 Zero-modes of the Myers-Perry family 

Let us now examine the particular form of Hij in the case of Myers-Perry solutions. For 
D = 5, there is no ultraspinning region, i.e. Hij is always positive definite. For D = 6,7, 



we find a more interesting behaviour, as shown in Figure 6.1 We parameterise the black 
hole by the horizon radius r+ and the spin parameters ai ~ Ji/M (defined in Section 
1.3.1), rather than M and Jj, in order to produce clearer figures. For D = 6, Hij is a 



2x2 matrix, and the ultraspinning "surface" is a closed curve that encloses the origin. 
Inside this curve, both eigenvalues of Hij are positive but outside one is positive and one 
is negative. Notice that all extreme black holes are ultraspinning. 
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For D = 7, Hij is a 3 x 3 matrix. The ultraspinning surface encloses the origin. 
However, now there is another surface on which a second eigenvalue of Hij changes from 
positive to negative. This surface lies between the ultraspinning surface and the surface 
corresponding to extreme black holes, and touches the ultraspinning surface at a point 
with equal a^. For generic Oj, if one gradually scales up the a, then one eigenvalue of Hij 
becomes negative as the ultraspinning surface is crossed, and another eigenvalue becomes 
negative as the other surface is crossed. As each surface is crossed, the black hole should 
develop a new stationary negative mode, and the associated black branes should develop 
a new instability. 

Some D = 7 black holes are special, e.g. solutions for which one of the angu- 
lar momenta vanishes do not intersect the second surface. For instance, consider the 
singly-spinning case. All such black holes must possess the thermodynamic negative mode 



predicted in Section 6.1.2 , and a second negative mode should appear at the ultraspinning 



surface, but there should be no further thermodynamic negative modes. This is consistent 
with the results of Ref. [72] , where it was found that a stationary negative mode does 
indeed appear at the ultraspinning surface. As we shall see in the next Chapter, new and 
infinitely many stationary non-thermodynamic negative modes appear at larger angular 
momenta, the first of which marks the threshold of instability of the black hole. 



Using the expressions in Section 1.3.1 we can actually derive an explicit expression 



for the reduced thermodynamic Hessian, 

{D - 2)^ 



Hi 



+ 2 

K 



+ 



{rl + afy {rl + aj)^ 2r+ 



K 



(6.11) 



where there is no sum over i,j and we have defined Cl'^ = il?. The matrix is positive 
definite in the static case, aj = 0. 
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In the singly-spinning case, say ai / 0, we have 
2{D -2){D- 3)^r+(r2 + 



H 



11 



H 



M [{D - 3)r2 +{D- 5)a2]3 
(Z)-2)7r, 



[(Z)-3K-(Z)-5)a^], 



for (i,j)/(l,l). 



3.12) 



There is a single zero- mode in Z? > 5, as we see in Figure 6.1, which occurs for 



a 



D-3 
D-5 



3.13) 



The associated eigenvector is dn, so that the angular momenta which vanish in the back- 
ground solution are not excited by the perturbation, i.e. the zero-mode keeps the black 
hole within the singly-spinning Myers-Perry family. The ultraspinning regime occurs for 



rotations higher than (6.13). 



In the equal spins case, we have 
{D - 2)7r 



H 



Mk 



s 12^ 



2r+ 



{rl + a2)2 



1 ^ nn^ 

2r+ K 



3.14) 



where Qij = 1 yi,j- An eigenvector Vi of Hij must then be an eigenvector of Qij, which 
leaves only two options: the eigenvector is such that Vi = V^i, or is such that Vi = 0. 
In the former case, there can be no zero-mode, since this would require 



+ 2 



2r, 



1 nfi^ 

2r-i_ K 



n = 0. 



3.15) 



which can be simplified to = —{D — 3)r^/(2 — e), and thus has no real solution for a 
(recall that e = 0, 1 for odd and even D, respectively). Hence no instability occurs for 
modes that preserve the equality between the spins, which is consistent with the results 
However, the eigenvectors satisfying V^, = do change sign once at \a\ = r+j^ 



of 
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There is one such eigenvalue in D = 6, and in D = 7 there are two, as shown in Figure [6T 



^Notice that, in Chapter [sj we will use the radial variable of Ref. 120 , which is related to the variable 
presently used by f'^ = + a^, so that the ultraspinning surface is at \a\ = f+/\/2. 
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The associated eigenvectors break the symmetry between the spins, so that the perturbed 
black hole is a Myers-Perry solution which is not in the equal spins sector. Due to the 
presence of the extremality bound, it is not clear whether the rotation can be sufficiently 
high to excite classical instabilities in the ultraspinning regime. However, we shall see in 
Chapter [s] that this indeed happens for D = 9 (and, we believe, for higher D, at least in 
the odd D case for which the solutions are codimension-1). 



6. A Appendix: Thermodynamic determinants 



We wish to show that 

det(-5,^) = - p_3)^^ det(/7,,) , (6.16) 

where S^^ = d^S/dx'^dx^, with x° = (M, J^), and Hij = d'^{-S) / dJjdJj = -Stj. This 
identity follows directly from the linear algebra result ( 2.45| ), if we take Cij = Hij. The 
proportionality coefficient is given by 

a - -5oo - H, 5o.5o, " " (^J/ l^J,, A^J. ' ^''''^ 

where we have used the first law of thermodynamics, dS = j3dM — l5Q.idJi. Since 

( dm,) \-'_( dj, \ (dJA _( dj, \ ( dm,) 



(6.18) 

we get 



The identity 

dl\ ^(dl^\ _(dl\ (dJA 

dMjj \dMj^^ \dJ,) ,,\dM) ^ 
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and the ^oi = Sio "Maxwell relation", {d{pnj)/dM)j = (5/3/5 Jj)j,^, further imply that 



dl3\ 



w) 



-1 



I3n 



3.21) 



Following the steps taken between Eqs. (6.4) and (6.6), we finally obtain 



CJ = — 



(D - 3)MT 



< 0. 



3.22) 



These steps require that the Smarr relation (6.5) is valid, i.e. they apply only to asymp- 
totically flat vacuum black holes. 



Chapter 7 

Singly-spinning Myers-Perry 
instability 

In this Chapter, based on |72|, we will confirm numerically the conjecture that Myers- 
Perry black holes with a single spin become unstable at a sufficiently large value of the 
rotation, and that new black holes with pinched horizons appear at the threshold of the 
instability. We search numerically, and find, the stationary axisymmetric perturbations 
that mark the onset of the instability and the appearance of the new black hole phases. 



7.1 Introduction 

In higher-dimensional spacetimes a vast landscape of novel black holes has begun to be 



uncovered, as discussed in Section [2.1.2[ Its layout - i.e. the connections between different 
classes of black holes in the space of solutions - hinges crucially on the analysis of their 
classical stability: most novel black hole phases are conjectured to branch-off at the thresh- 
old of an instability of a known phase. Showing how this happens was an outstanding 
open problem first addressed in Ref. [72j. 

The best known class of higher-dimensional black holes, discovered by Myers and 
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Perry (MP) in 43 and described in Section 1.3.1, appears in many respects as the nat- 
ural generalisation of the Kerr solution. In particular, the event horizon is topologically 
spherical. However, the actual shape of the horizon can differ markedly from the four- 
dimensional one, which is always approximately round with a radius parametrically ~ M. 
This is not so in D > 6. Considering for simplicity the case where only one spin J is 
turned on (of the n = [{D — l)/2j independent angular momenta available), it is possible 
to have black holes with arbitrarily large J for a given mass M. These black holes fall 
into the ultraspinning regime defined in the last Section. The horizon of these black holes 
spreads along the rotation plane out to a radius a ~ J/M much larger than the thickness 
transverse to this plane, ~ {M^ / J'^y^^^~^\ 

This fact was picked out by Emparan and Myers fiT] as an indication of an instabil- 
ity and a connection to novel black hole phases. In more detail, in the limit a — )• oo with 
fixed, the geometry of the black hole in the region close to the rotation axis approaches that 
of a black brane. Black branes are known to exhibit classical instabilities [lOOj , at whose 



threshold a new branch of black branes with inhomogeneous horizons appears 1104, 105 



Ref. 71 conjectured that this same phenomenon should be present for MP black holes 
at finite but sufficiently large rotation: they should become unstable beyond a critical 
value of and the marginally stable solution should admit a stationary, axisymmet- 

ric perturbation signalling a new branch of black holes pinched along the rotation axis. 
Simple estimates suggested that in fact (a/r_|_)crit should not be much larger than one. 
As a/r+ increases, the MP solutions should admit a sequence of stationary perturbations, 
with pinches at finite latitude, giving rise to an infinite sequence of branches of 'pinched 
black holes'. Ref. |62| argued that this structure is indeed required in order to establish 
connections between MP black holes and the black ring and black Saturn solutions more 



recently discovered. See the idea in Figure 7.1, which is analogous to Figure 2.2 for the 



non- uniform branes, in Section 2.1.2 See also how this picture arises in the black hole 



phase diagram of Figure 7.2 



Our main result is a numerical analysis that verifies the conjecture of Ref. \7T\. 
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Figure 7.1: Conjectured connections between a D > 6 singly-spinning MP black hole and a 
black ring (A), a black Saturn (B) or a di-ring (C), through a family of pinched solutions. 
The bifurcation from the MP family occurs at the points A, B and C, respectively, in the 



phase diagram of Figure 7.2 



7.2 Background solution 



The solution for a MP black hole rotating in a single plane is 

ds"^ = - J^!'^\^ \dt + a sin^ 6* #1 ^ + J!!" ^, [(r^ + a^)d(/) - a dtl ^ 



+ ^^|^dr2 + S2(r, 0) + cos" 6 dO^^.^) , (7.1) 



„2 „„„2 



S(r, 0)2 = + cos^ , A(r) = + - . (7.2) 



D-3 



The parameters here are the mass-radius and the rotation-radius a, 



levrGM D-2 J 
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J 



Figure 7.2: Diagram of entropy versus angular momentum, at fixed mass, for singly- 
spinning MP black holes in D > 6 illustrating the conjecture of [7l] (see also (62j): at 
sufficiently large spin the MP solution becomes unstable, and at the threshold of the 
instability a new branch of black holes with a central pinch appears (A). As the spin 
grows new branches of black holes with further axisymmetric pinches {B, C, . . .) appear. 
We determine numerically the points where the new branches appear, but it is not yet 
known in which directions they run. We also indicate that at the inflection point (0), 
where {d'^S/d,P)M = 0, there is a stationary perturbation that should not correspond to 
an instability nor a new branch, but rather to a thermodynamic zero-mode that moves 
the solution along the curve of singly-spinning MP black holes. 



The event horizon lies at the largest real root r = r+ of A. 

The linearised perturbation theory of the Kerr black hole {D = 4) was disentangled 
171| using the Newman-Penrose formalism. Attempts to extend this formalism to 



m 



decouple a master equation for the linear perturbations of MP black holes have not yet 
been successful; see [172 for recent progress. Moreover, even though some subsectors of 
perturbations for some classes of MP black holes had been decoupled before the work 
described here [72], e.g. [Ml [Tl4lp0l[T73l[l74 ], none of them had shown signs of any 
instability and they do not contain the precise kind of perturbations we are interested in|^ 
Thus we take a more frontal numerical approach to a full set of coupled partial differential 
equations (PDEs). 



^Notice that 
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and 



68 



on which the next Chapter is based, are posterior to 



72 , 
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7.3 The eigenvalue problem 



We intend to solve for a stationary linearised perturbation h^b around the background 



(7.1 ). Choosing the traceless-transverse (TT) gauge, h^-^ = and V^hab = 0, the equations 
to solve are 

{/\Lh)ab = , (7.4) 



where A/, is the Lichnerowicz operator defined in (2.2). Actually, we solve the more 
general eigenvalue problem 



ALh)ab = -k^hab. (7.5) 



As we discussed in the previous Chapter, this problem appears in two contexts. Eq. (7.5) 
determines the stationary perturbations of a black brane obtained by adding flat directions 
z to the line element (7.1), for perturbations with a profile e^^'^hab- Thus such modes 



with k > correspond to the threshold of the Gregory-Laflamme instability of black 



branes 100 . The same equations also describe the negative modes of quadratic quantum 



corrections to the gravitational Euclidean partition function 131|132| . In Chapter [4| based 



on 



1561, we studied this problem for the Kerr black hole, and showed the existence of a 



branch of solutions extending the negative Schwarzschild mode {k^ch 7^ 0) [132] to finite 
rotation, with k growing as the rotation increases towards the extremality bound. 

If the ultraspinning instability is present for MP black holes in D > 6, then, in 
addition to the analogue of the branch studied in [156| , a new branch of negative modes 
extending to /c = must appear. The eigenvalue k = corresponds to a (perturbative) 
stationary solution with a slightly deformed horizon. In fact, as explained above, we 
expect an infinite sequence of such branches that reach k = at increasing values of the 
rotation. The solutions for A; > imply new kinds of Gregory-Laflamme instabilities and 
inhomogeneous phases of ultraspinning black branes (see also 108] ) . 



Another reason to consider (7.5) instead of trying to solve directly for = is that 
there exist powerful numerical methods for eigenvalue problems that give the eigenvalues 
k together with the eigenvectors, i.e. the metric perturbations; see the Appendix at the 
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end of this thesis. 

The modes we seek preserve the SO{D — 3) x SO{2) rotational symmetries of the 
singly-spinning MP solution and depend only on the radial and polar coordinates, r and 
6 |71|. Thus we take the ansatz 



+ 



A(r) 
~S2(r,e) 

A(r) 



dt + a sin 9 e 



5uj 



2 sin 



[r' + a^)d(l)-ae-^^ dt 



e^^« [dr + 5x sin 6 dOf + ^^(r, 9) e^^"^ dB^ + cos^ 9 e^* dnf^_^j , (7.6) 



where {Suq, Si'i, SfiQ, 5fii, 5uj, 5x, S^}, which are functions of (r, 9), describe our perturba- 
tions. These functions will be determined numerically by solving the eigenvalue problem 



(7.5) subject to appropriate boundary conditions. After imposing the TT gauge, Eq. (7.5) 



reduces to four coupled PDEs for SfiQ, 5fii, 5x and 6^ (the TT conditions then give Suq 
5vi and 6uj)^ in a procedure equivalent to the one described in Chapter |4| based on 156 
for the Kerr-AdS black hole. 



The boundary conditions are discussed in detail in the Appendix 7. A For the nu- 



merical implementation, it is convenient to use the new radial and polar variables 



y = l 



cos 9 . 



(7.7) 



which are dimensionless and take values in the interval [0, 1] . Dirichlet boundary conditions 
are simpler to implement, and so we will use the following perturbation functions: 



1 



r+ 



1 )x{l- x)6fii{y,x) , 

r J 



q2{y, x) = (1 - x) 6x{y, x) , 



Q4:{y,x) 



1 -] x(l- x) 5^(y,x) , 

r / 

(7.8) 



which vanish at the boundaries r = r+,r = oo,x = and x = 1 when the boundary 
conditions determined in the Appendix |7. A are enforced. 

It is important to ensure that the eigenmodes we find are not pure gauge, hab = 
V(^aCb)- This is trivial when A; / 0, since AiV(a^b) = 0, but not for zero-modes. We can 
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D = 7 
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D = 7 (blue), 8 (green), 9 (red) 
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Figure 7.3: Negative eigenvalues for the MP black hole in D = 7 (left) and, for comparison, 
in = 7,8,9 (right). 



prove that in the TT gauge, where the equations for the gauge parameter are separable, 
pure gauge perturbations within our ansatz necessarily diverge at the horizon or at infinity 



(the proof is along the lines of the equal spins MP case, in the Appendix 8.B.2). Thus, 
with our boundary conditions, the zero-modes are never pure gauge. 



7.4 Results and discussion 



We use again, as in Chapter [4] describing 156 , a spectral numerical method to solve the 



eigenvalue problem; see the Appendix at the end of the thesis. We have carried out the 
calculations for D = 7,8, 9. The cases D = 5 (where there is no ultraspinning regime) and 
D = Q were left out in 72 due to numerical difficulties which have since been solvedj^ 



For D = 5, there is no instability of this type. For D = 6, the results conform to our 
expectations, being similar to the ones for D = 7, 8, 9 discussed here. 

The results for D = 7, 8, 9 are displayed in Figure 7.3 The negative eigenvalue —k'^ 



is plotted as a function of the rotation parameter a. We normalise k and a relative to the 
mass-radius r^, which is equivalent to plotting their values for fixed mass (or mass per 
unit length, in the black string interpretation). As described above, the leftmost curve. 



Ref. 175 was published after the submission of this thesis. 
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which does not reach A: = 0, is the higher-dimensional counterpart of the Kerr negative 
mode, and the eigenvalues k are the wavenumbers of the Gregory-Laflamme threshold 
modes at rotation o. As we consider larger rotation, new branches of negative modes 
appear that intersect A; = at finite a/vm- We label these successive branches with an 
integer i = 1,2,3,..., and refer to them as 'harmonics'. The values of a/r^ at which the 



stationary perturbations appear are listed in Table 7.1 



D 


(aAm)k=i 


{a/rjn)\i=2 


(a/?'m)k=3 


7 


1.075 


1.714 


2.141 


8 


1.061 


1.770 


2.275 


9 


1.051 


1.792 


2.337 



Table 7.1: Values of the rotation a/vm for the first three zero-modes {k = 0). The 
estimated numerical error is ±3 x 10~^ in D = 7 and ±5 x 10~^ in I? = 8, 9. 



It is important to note that the eigenmode A: = of the harmonic i = 1 does not 
correspond to a new stationary solution. Instead it is a thermodynamic zero-mode that 
marks the onset of a new local thermodynamic instability, and takes the solution to an 
infinitesimally nearby one along the family of MP black holes. The location of this zero- 



mode is predicted by the reduced thermodynamic Hessian Hij, introduced in Section 6.2 



and corresponds to the inflection point of the curve S{J) at fixed M (point in Figure 7.2 ) 



According to the expression (6.13), this point is given by 



a 



D-3 



D 



1=1 



2(L>-4) \D-5 



. -P-5 



(7.9) 



For rotations larger than this value, the black holes are in the ultraspinning regime, as 



defined in the last Chapter. The values of (a/rm) given by (7.9) for D = 7, 8,9 agree with 
the central values of the numerically-determined rotations {a/rm) for 1 = 1 (first column 



in Table 7.1) up to the third decimal place. This is a very good check of the accuracy of 



our numerical method. 

The k = Q eigenmodes of the higher harmonics, £ > 2, do not admit this interpreta- 
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tion as perturbations along the MP family of solutions and thus correspond to genuinely 
new (perturbative) black hole solutions with deformed horizons. Their appearance con- 
forms perfectly to the predictions in [7l] and [62] . It is then natural to expect, although 
our approach does not prove it since it only captures zero-frequency perturbations, that 
the harmonic £ = 2 signals the onset of the instability conjectured in [7l]. The A: = 
eigenmodes for higher harmonics confirm the appearance of the sequence of new black 
hole phases as the rotation grows. 

To visualise the effect on the horizon of the perturbations that give new solutions, 
and provide further confirmation of our interpretation, we draw an embedding diagram 
of the unperturbed MP horizon and compare it with the deformations induced by the 



ultraspinning harmonics £ > 2. This is best done using the embedding proposed in 176 



which has the advantage of allowing one to embed the horizon along the entire range 
< 9 < tt/2 for any rotation, although at the cost of stretching the pole region, which 



acquires a conical profile. We do it for the £ = 2,3,4 ultraspinning harmonics in Figure 7.4 



In spite of the distortion created by the embedding, the effect of the perturbations is clear: 
1 = 2 modes create a pinch centred on the rotation axis 9 = 0; i = 3 modes have a pinch 
centred at finite latitude 9; i = A modes pinch the horizon twice: around the rotation 



axis and at finite latitude. These are the kind of deformations depicted in Figures 7.1 



and |7.2[ To better identify the number of times that the perturbed horizon crosses the 
unperturbed solution, in these figures we also plot the logarithmic difference between the 
two embeddings. 



Ref. 71 gave several arguments to the effect that critical values a/r^ close to 1 
were to be expected. In particular, it was pointed out that the change in the behavior of 
the black hole from 'Kerr-like' to 'black brane-like' could be pinpointed to the value of the 
spin where the temperature (i.e. surface gravity) reaches a minimum for fixed mass, which 
is the same, for solutions with a single spin, as the infiection point of S{J). As we have 
argued, the zero-mode at this solution should not signal an instability. The i = 2 mode at 
the threshold of the actual instability instead appears at larger rotation, well within the 
ultraspinning regime defined in the last Chapter. We conjecture this to be true in general: 
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Figure 7.4: Top plot: Embedding diagram at {a/rm)\E=2 of the D = 7 black hole horizon, 
unperturbed (solid), and with the first unstable harmonic perturbation (i = 2, k = 0) 
(dashed). The embedding Cartesian coordinates Z and X lie respectively along the rota- 
tion axis 9 = and the rotation plane 9 = 7r/2. We also show the logarithmic difference be- 
tween the embeddings of the perturbed (^£=2) and unperturbed (Zq) horizons. The spikes 
represent the points where the two embeddings intersect. The perturbation has two nodes, 
so the horizon squeezes around the rotation axis, then bulges out, and squeezes again at 



the equator, as in the conjectured shape A in Figures 7.1 and 7.2 Bottom-left plot: for 



£ = 3, between the first two nodes of the perturbation, the horizon has a pinch, like shape 
B. Bottom-right plot: for £ = 4, the four nodes deform the horizon into shape C. 
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instabilities of MP black holes which have stationary threshold modes can only occur in 
the ultraspinning regime. 

We have identified the points in the phase diagram where the new branches must 
appear, but we cannot determine in which direction these run. This requires calculating the 
area, mass and spin of the perturbed solutions. However, for any k ^ — and numerically 
we can never obtain an exact zero — the linear perturbations decay exponentially in the 
radial direction, and thus the mass and spin, measured at asymptotic infinity, are not 
corrected. Moreover, the modes i > 2 should not affect the asymptotic charges, which 
are usually associated with the lowest modes (we will confirm this explicitly in the next 
Chapter for MP black holes with equal spins). It seems that in order to obtain the 
directions of the new branches one has to go beyond our level of approximation or adopt 
a different approach. 

The new i > 1 branches extend to non-zero eigenvalues k. These imply a new ultra- 
spinning Gregory-Laflamme instability for black branes, in which the horizon is deformed 
not only along the extended directions of the branes, but also along the polar direction of 
the transverse black hole. Observe that, even if the i = 1 mode does not signal a classical 
instability of the MP black hole, the modes i = 1, k > are expected to correspond to 
thresholds of instabilities of MP black branes. At a given rotation, modes with larger £ 
have longer wavelength k~^ and so the branch ^ = 1 is expected to dominate the instabil- 
ity. The growth of k with a can be understood heuristically, since as a grows the horizon 
becomes thinner in directions transverse to the rotation plane and hence it can fit into a 
shorter compact circle. 



After this work was released 72 , Shibata and Yoshino [112 113 were able to find an 
instability for non-axisymmetric perturbations of singly-spinning MP black holes, using 
a non-linear numerical method. This instability seems to kick in for lower values of the 
rotation than the mode 1 = 2 found here. In particular, it also occurs in the D = 5 
case, where no instability of the type we studied occurs. Since the perturbations are not 
axisymmetric, the thermodynamics arguments on which our conjecture is based do not 
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apply. Furthermore, if the threshold mode is not stationary, it is not associated with new 
families of stationary black hole solutions bifurcating from the MP family. 

To finish, we mention that pinched phases of rotating plasma balls, dual to pinched 
black holes in Scherk-Schwarz compactifications of AdS, have been found [177 , as well 
as new kinds of deformations of rotating plasma tubes |110| and rotating plasma ball 
instabilities [178|179| . The relation of our results to these and other phenomena of rotating 
fluids is an interesting problem. 



7.A Appendix: Boundary conditions 



Boundary conditions at the event horizon 

We want to determine the conditions for the metric perturbations hab to be regular 



at the event horizon. As we discussed in Section 3.3, requiring that the traceless-transverse 
perturbation is a regular 2-tensor on the black hole background gab is more restrictive than 
requiring that the perturbed geometry is regular. In particular, such a perturbation cannot 
change the temperature T and angular velocity of the background solution. From the 
Euclidean perspective, if hab does not obey the periodic identifications of the imaginary 
time and rotation angles of the background geometry, it clearly cannot be a regular tensor 
on that geometry. 



Let us first discuss the unperturbed background geometry (7.1). Near r ~ r+, we 
can write A(r) = A'(r+)(r — r_|_) + 0[(r — r+)^], with A'(r_|_) > 0, and the near horizon 
geometry of ( |7.1[ ) reads 

{rl+a?f A'(r+)(r-r+) 

(7.10) 
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This suggests the introduction of a new azimuthal coordinate 

^ = (l)-nHt, nH = ^^, (7.11) 

with period Ac/) = 2ir. Now we perform a Wick rotation into Euchdean time r and define 
a new radial coordinate p according to 

t = -iT, r = with T ^ ^' 



r = r+ + ^^p\ (7.12) 



The near horizon geometry is then given by 

ds\^^^ S2(r+,^)[/>W + dp2]+s2(r+,e)d^2 



(7.13) 



+ S2(r+,g) # +<cos^^dO^^_4). 

The periodic identifications required to ensure regularity are 

(?,^) ~ (f + 27r,(^) ~ (?,0 + 27r). (7.14) 

These correspond to the already mentioned identifications which give the temperature 

T = I3~^ and angular velocity Q of the solution, 

(r, (/.) - (r + (/) - mp) ~ (r, + 27r) . (7.15) 

Let us now introduce Cartesian coordinates in the (r, p) plane by taking r = 
Arctan(y/x) and p = y^x^ + y^. The following 1-forms are manifestly regular: 

E'^ = p'^dr = xdy — ydx , E^ = pdp = x dx + y dy . (7.16) 
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In terms of these 1-forms, the metric perturbation reads 



6vq — 2a^ sin^ d 



't? (r+,g)A^(r+)+2r+ (r^ + a^) 
S4 (r+,0)A' (r+) 



5a; 



(7.17) 



+ (r+, 0)5^1 de'^ + 



A'(r+) 

(r^ + a^)^ sin^ i 
S2 



S2(r+,e)A'(r+) p2 



(D-4) 



Regularity then requires that 



5x 



r=r^ 



6lo\ 



0, 6i^o\ 







r=r+ 

are finite . 



and 



(7.18) 



The first and second conditions ehminate irregular contributions arising respectively from 



the terms E^dO and E'^dcj) in (7.17). The second and third conditions guarantee that there 



is no conical singularity in the (t, p) plane, since the first two terms in (7.17) then read 
simply 5fio (^p'^dr^ + dp^), which is manifestly regular. 



Boundary conditions at the 6 = equator 



We follow here the same strategy. Let us focus on = 0, and introduce the coordi- 



nate 



cos (9 = 1 - - X 



(7.19) 



The background geometry (7.1) near x = is given by 



rj Q \ ^ 



A(r) , 2 
^ ^ dt^ + 2a[l 



J.2 J_ g2 



A(r) 
+ a? 



dtd(j) 



+ (r2 + a^) {dx' + x'dcb') + '^^^dr' + ^'dnf^-,) 



(7.20) 
(7.21) 



which is manifestly regular given that (p has period 2it. 



7.A. APPENDIX: BOUNDARY CONDITIONS 



127 



Cartesian coordinates on the (x, </>) plane can be chosen by taking (p = Arctan(7//a;) 
and X = + and the fohowing 1-forms are manifestly regular: 

= xdx = X dx + y dy , E'^ = x^dcj) = xdy — ydx . (7.22) 

The metric perturbation then reads 



hab dx" dx^lg^Q ~ 2 I 2 '^'^o + A( ^ '^'^^ ^ A( \ ^^'^'^ 



+ A(r) A(r) 

+ [(r^ + + A(r)) Jo; - (r^ + a^) <^i.i + A(r)5z.o] 



_|_ 

+ (r2 + a^) [5iiidx^ + <5zvi;t2^</.2] + 5^ dnlj^_^~^ , 

(7.23) 

and regularity requires that 

" " (7.24) 

^xle^o' '^^le=o' '^^o|g=o' '^'^o|e=o' '^^l0=o are finite. 



Boundary conditions at the 9 = Tr/2 equator 



We introduce the new coordinate x = cos 6. The geometry (7.1) in the neighbour- 
hood of the rotation plane 9 = tt /2 \s given by 



ds^ 



l0~7r/2 



A(r) 



[r^ + - A(r)] + 



(^2^„2^2_^2^(^) 



dx^ + 



(7.25) 



which is manifestly regular given that is the line element of an 

Introducing the manifestly regular and smooth 1-forms, 



E"" = xdx. E^ = x'^dVL 



(d-A) , 



(7.26) 
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the metric perturbation reads 

h/iudx dx |g^^/2 ~ ^ '^^ 

[r"^ + a'^f 5ui-a'^l^{r)[5vo + 25uj) ^ 

(7.27) 



+ [(r^ + + A(r)) 5a; - (r^ + a^) 5ui + A(r)5i'o] dtdcp 



r 

^2 ^2 / 



A(r) A(r) V 2; 

Regularity requires that 

Sx\a-7L=0, = 6iJ.i\„_;^ and 

IP— 2 ''"2 ''"2 

5a;|g^2- , S^o\g^T^, Suo\g^^ , 5i^i|g^7r are finite. 



(7.28) 



Boundary conditions at the asymptotic region r — t- 00 



At spatial infinity, r — )• 00, the solutions to Eq. (7.5) behave as hat oc e^'^'''. There- 
fore, regular perturbations must vanish at infinity. 



Chapter 8 



Equal-spinning Myers-Perry 
instability 



In this Chapter, based on Ref. [68j, we present the first example of a linearised gravitational 
instability of an asymptotically flat vacuum black hole. We study perturbations of a Myers- 
Perry black hole with equal angular momenta in an odd number of dimensions. We find 
no evidence of any instability in five or seven dimensions, but in nine dimensions, for 
sufficiently rapid rotation, we find perturbations that grow exponentially in time. The 
onset of instability is associated with the appearance of time-independent perturbations 
which generically break all but one of the rotational symmetries. This is interpreted as 
evidence for the existence of a new 70-parameter family of black hole solutions with only a 
single rotational symmetry. We also present results for the Gregory-Laflamme instability 
of rotating black branes, demonstrating that rotation makes black branes more unstable. 



8.1 Introduction 

In the previous Chapter, based on Ref. [72], we presented very strong evidence for the 
existence of an ultraspinning instability of singly-spinning Myers-Perry black holes. How- 
ever, no actual instability, i.e. a perturbation growing in time, was demonstrated. In this 
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Chapter, based on Ref. [68] , we shall demonstrate that some Myers-Perry (MP) black 
holes do admit perturbations that grow exponentially in time, thereby providing the first 
example of a linearised gravitational instability of an asymptotically fiat vacuum black 
hole solutionii] 

We shall exploit the idea introduced in Ref. |120| of considering MP solutions with 
enhanced symmetry. The generic MP solution (Section 1.3.1| ) has isometry group i?x C/(l)" 
where R corresponds to time translations and n = \_{D — l)/2j . However, this is enhanced 
when some of the angular momenta coincide. In particular, for odd D, the MP solution 
with all angular momenta equal ( Jj = J/n Wi) has a much larger R x U (N + 1) isometry 
group, where D = 2N + 3. Furthermore, the solution is cohomogeneity-1, i.e. it depends 
only on a single coordinate. The metric involves a fibration over complex projective space 
CP^ . Gravitational perturbations of this solution can be decomposed into scalar, vector 
and tensor types according to how they transform under isometrics of CP^ . The tensors, 
which exist only for N > 2 {D > 7), were studied in Ref. [120 and no evidence of any 
instability was found. The special case of D = 5, for which only scalar perturbations exist, 
was studied in Ref. [114 . Again, no evidence of any instability was found. 

In this Chapter, we shall study scalar-type perturbations of these cohomogeneity- 
1 MP black holes. The symmetries enable the problem to be reduced to coupled linear 
ordinary differential equations (ODEs) which we solve numerically. We find no evidence of 
any instability for D = 5 (consistent with Ref. frT4|) or D = 7. However, for D = 9, when 
J exceeds a certain critical value Jcrit) there is a perturbation that grows exponentially in 
time, i.e. an instability. We believe that such an instability will exist for all (odd) D > 9 
although we have demonstrated this only for D = 9. 

As expected, the onset of instability is associated with the appearance of a stationary 
zero- mode of the MP solution with J = Jcrit- This zero- mode is interesting for another 
reason. It has been proven that a stationary, rotating black hole must admit a rotational 



iRef. 



112 



where an instability of singly-spinning Myers-Perry black holes for non-axisymmetric per- 



turbations was found (see also 113 ), was released shortly before [68| . However, the numerical analysis 
in 112 113 is fully non-linear, despite dealing with small perturbations. It would be very interesting to 



analyse the instability found there at the linear level, so that the threshold mode can be better studied. 
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isometry (i.e. a U{1) isometry) [55 66 . This is known as the rigidity theorem. However, 
all known higher-dimensional black holes have multiple rotational isometries (e.g. MP 
black holes have [{D — l)/2j commuting U{1) isometries), i.e. more symmetry than one 
expects on the basis of general arguments. Therefore, it has been conjectured that there 
exist solutions with less symmetry than any known solution, specifically solutions with a 
single rotational symmetry [67]. 

It was proposed in Ref. |67| that one could seek evidence for the existence of such 
solutions in the same way that the first evidence was obtained for the existence of non- 
uniform black string solutions. For black strings, the static zero-mode associated with 
the onset of the Gregory-Laflamme instability was conjectured to describe the "branching 
off' of a new family of non-uniform black string solutions from the already known branch 
of uniform solutions 101| . Perturbative [104j and numerical 105 work subsequently 
confirmed that this was correct. For rotating black holes, the idea proposed in Ref. [67] 
is to look for a stationary zero-mode of a MP solution. By analogy with the black string 
example, this could be interpreted as the branching off of a new family of solutions. If the 
zero-mode preserves only a single rotational symmetry then this would be evidence for the 
existence of new black holes with just one rotational symmetry. 

In our case, the stationary zero-mode generically preserves only a single rotational 
symmetry. Therefore we conjecture that there exists a family of stationary black hole 
solutions with just a single rotational symmetry, that bifurcates from the cohomogeneity- 
1 MP black hole solution at J = Jcrit- In fact, we find not just one stationary zero-mode, 
but a large family, corresponding to all scalar harmonics on CP^ with a certain eigenvalue. 
For D = 9, we shall argue that the new family of solutions will involve 70 independent 
parameters, considerably more than the 5 parameters required to specify the MP solution! 
If correct, this implies that any hope of specifying higher-dimensional black holes uniquely 
using just a few parameters is bound to fail. Of course, the new black hole solutions may 
turn out to be unstable themselves. 

Recently, Ref. [65 1 has provided other evidence for the existence of higher-dimensional 
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black holes with a single rotational symmetry. Using the "blackfold" approach of Refs. |63l 



64 , approximate solutions were constructed for D > 5 that describe "helical" black rings. 
The blackfold approximation is based on the fact that higher-dimensional black holes can 
have widely separated horizon scales. The results of the present Chapter concern black 
holes that lie outside the regime of validity of this approximation. Our results complement 



those of Ref. 65 because we are presenting evidence for topologically spherical black holes 



with a single rotational symmetry, whereas Ref. 65 considered black rings. 



A difference between the results to be discussed here and the results in the previ- 
ous Chapter concerns the nature of the stationary zero-mode. In the present case, the 
unstable perturbation breaks all but one rotational symmetry. However, the stationary 
zero-modes found in the previous Chapter for singly-spinning black holes preserved the 
isometrics of the background, as Ref. [7^ had conjectured. Singly-spinning black holes 
exhibit symmetry enhancement for D > 6. They are cohomogeneity-2 with isometry group 
R X C/(l) X SO{D - 3) where SO{D - 3) has 5^"^ orbits^ While these isometrics are pre- 
served by the zero-modes found previously, the results in this Chapter raise the question of 
whether singly-spinning MP black holes might admit further stationary zero-modes that 
break some of their symmetry, and provide further evidence for new black hole solutions 
with reduced symmetry. 

A corollary of our approach is the first data for the Gregory-Laflamme instability 
[100 for rotating black branes. We consider black branes obtained as the product ( 1.11 ) of 



a cohomogeneity-1 MP black hole with flats directions. We argue that such solutions are 
always classically unstable. Our numerical results demonstrate that the branes become 
more unstable as the angular momentum increases: the instability becomes stronger (i.e. 
it occurs on a shorter time scale) and the critical wavelength of unstable modes decreases, 
as the angular momentum increases. 



^One can decompose metric perturbations of these solutions into scalar, vector and tensor types using 
this SO{D—3) symmetry. The tensors, which exist only for D > 7 have been studied previously in Ref. [174| 
and show no evidence of any instability. However, tensor perturbations arise from deformations of the 5* 
part of the metric, whereas the expected ultraspinning instability should arise from perturbations of the 
metric transverse to The results of Ref. 72 (Chapter [t]) indicate that the instability should be a 

scalar-type perturbation. 
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An important feature of cohomogeneity-1 MP black holes is that they exhibit an 
upper bound on their angular momentum (Section |1.3.1 ). Solutions saturating this bound 
are extreme black holes. Because of this upper bound, it is not obvious that such black 
holes should exhibit the type of behaviour discussed in Ref. |71:j. However, we were 



motivated by the observation in Section 6.2.2 that, for D > 7, cohomogeneity-1 MP black 
holes do satisfy the ultraspinning criterion once J exceeds a critical value Juitra (for given 
M). This criterion says that a classical instability with a stationary threshold can occur 
only if the Hessian matrix Hij = {d"^ {— S) / d Jid J j) m ■, where S is the black hole entropy, 
fails to be positive definite. Hence, there might be an ultraspinning instability for J > Jcnt 
where Juitra < Jcrit < ^extreme- Our results show that there is no instability for D = 7 but 
an instability does occur for D = 9 and, we believe, for (odd) D > 9. 



For D = 5 MP black holes, Hij is always positive definite 72 so such black holes 
are never ultraspinning. However, for D > 6, singly spinning MP black holes with large 
enough angular momentum are, of course, ultraspinning and the numerical results of 
Ref. supply strong evidence that the instability appears only when Hij fails to be 
positive definite. 



This Chapter is organized as follows. Section 8.2 describes the cohomogeneity-1 
black holes that we shall study. In Section [8.3[ we explain our approach and discuss the 



results. The technical details of our work are presented in the later Sections 8.4 and 8.5 
and in the Appendices. 



8.2 Cohomogeneity-1 Myers-Perry black holes 

The Kerr solution was extended to higher dimensions by Myers and Perry [43] . The Myers- 
Perry family can be parameterized by a mass-radius parameter rM and [{D — 1) /2j angular 
momentum parameters Oj. In the particular case of equal angular momenta, ai = a, the 
solution in odd dimensions D = 2N + 3 is cohomogeneity-1. The metric can be written 
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asEl 

ds'^ = -f{rfdt^ + g{rfdr'^ + h{rf[d'4^ + Aadx" - n{r)dtf + r^abdx'^dx'' , {8.1] 

/ „2Af 2Ar 2\-l / r'^N 2\ 



f{r)= , , , f^(r) 



JIN 



g{r)h{r) ' j.2Nf^2 ' 

where is the Fubini-Study metric on CP^ with Ricci tensor Rab = 2(A^ + l)gab , and 
A = Aadx"" is related to the Kahler form J by dA = 2J. Surfaces of constant t and r 
have the geometry of a homogeneously squashed S^^"*"^, written as an fibre over CP^ . 
The fibre is parameterized by the coordinate tp, which has period 2tt. Explicit expressions 
for the metric gab and Kahler potential A of CP^ can be obtained through the iterative 
Fubini-Study construction summarized in Appendix |8.A 



The spacetime metric satisfies Rfj^u = and the solution is asymptotically flat. The 
event horizon is located at r = r_|_ (the largest real root of g~^) and it is a Killing horizon 
= dt + ^nd^ ) where the angular velocity of the horizon is given by: 



The mass M and angular momentum J, defined with respect to 5^, are 
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where A2N+1 is the area of a unit (2A'^ + l)-sphere. 

There is an extremality bound on the angular momentum which can be expressed 



^The radial coordinate used here can be related to the standard Boyer-Lindquist radial coordinate 
through + . In this Chapter, Latin indices will be used for coordinates on CP'^ , whereas 
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of 

Greek indices will be used for the black hole spacetime coordinates, unlike the previous Chapters. 
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as 



2 / s 2 



< 



N, 



or 



< 



N 



VmJ ~ \rM J (iV + l){A^+i)/^ • 



The solution saturating this bound has a regular, but degenerate, horizon. For fixed r^-, 



or tm, ultraspinning behaviour (Section 6.2.2) occurs for 



> 



ai 



or 



ai 



(-) >. 

\rM J \rM 



1 



2{N+1)/N ■ 



(8.5) 



Note that the range ai < a < Oext (for fixed r+) for which the black hole is ultraspinning 
becomes larger as N increases, and that ai = Oext if = 1, so there is no ultraspinning 
behaviour for D = 5. 



8.3 Strategy and Results 
8.3.1 Strategy 

Ref. [156| (Chapter [4]) introduced new numerical techniques for determining negative 
modes of rotating black holes. In Ref. [72] (Chapter [?]), these techniques were exploited to 
construct the stationary zero-mode expected to indicate the onset of an ultraspinning in- 



stability of a singly-rotating MP black hole. In the present Chapter, based on 68 , we shall 
determine the stationary zero- mode indicating the onset of instability for cohomogeneity-1 
black holes. However, our main achievement is to generalise these methods to demonstrate 
the existence of perturbations that grow exponentially in time. 

Our approach is explained in Chapter [6| We consider the eigenvalue problem 

{ALh)^^ = -k^hf,^, (8.6) 



where A/, is the Lichnerowicz operator (2.2) for the MP background, and h^i, are traceless- 



transverse perturbations (/i^^ = V^h^i^ = 0) of the MP black hole. This problem arises for 
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Gregory-Laflamme-type perturbations (2.3) of uniform black branes. Perturbations with 



non-zero k correspond to negative modes of Al, which, in the stationary case, may also 



correspond to negative modes (3.33) of the black hole partition function. The boundary 
conditions are that hf_iu should be regular on the future event horizon 'H'*' and vanishing 
at infinity. 

The strategy for studying perturbations of the black hole will be to seek a solution 



of (8.6), i.e. a negative mode of the black hole, and then vary the spin of the black hole 
until k vanishes, i.e. the negative mode becomes a zero-mode. This strategy is motivated 
by the availability of numerical techniques for solving eigenvalue equations of the form 



.6). Solutions with non-zero k correspond to perturbations of black branes. Therefore 
our method will yield results for the Gregory-Laflamme instability of rotating black branes 
as well as enabling us to search for black hole instabilities. 

We can Fourier analyse our perturbation in the time and directions, i.e. we assume 
that the dependence on t and ip is given by 



V oc e-'^'+'"'^ , (8.7) 



where m is an integer. As we shall explain in detail below, we can also decompose the 
perturbation into harmonics on CP^ . These can be of scalar, vector or tensor type. The 
tensors were considered in Ref. [120 . We shall restrict our attention to perturbations 



of scalar-type, which can be expanded in terms of scalar harmonics on CP^ . As usual, 
harmonics with different eigenvalue decouple from each other. The equations satisfied by 
h^y depend only on the eigenvalue of the harmonic in questionj^ Eigenvalues of the scalar 



Laplacian on CP^ are labelled by a non- negative integer k (see Section 8.4). Hence our 
perturbation is labelled by (w, m, k). 



Consider the (Lorentzian) negative mode equation (8.6). As we have explained 



above, this arises from classical perturbations of a rotating black brane. The usual ap- 



*In the last Chapter, we referred to the different negative modes as different "harmonics" . However 
the present 
codimension f 



in the present case, we do have a precise harmonic structure in CP^ , which is why the problem has 
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proach to this problem is to fix (k, m, k) and to determine u. However, our approach 
win be to fix {uj,m,n) and determine the possible eigenvalue(s) —k"^. In other words, 
we are determining the wavenumber k for which black brane perturbations with given m 
and K have time-dependence associated with the given cj. We shall fix the overall scale 
tm = 1 and determine the eigenvalue(s) —k"^ for fixed {uj,m,K) as a increases from to 
extremality. If k vanishes for some value of a then the associated black hole admits a 
zero-mode with the given values of {u, m, k) (of course it must be checked that this is not 
pure gauge). 



In searching for an instability, we are looking for solutions of (8.6) with lm(ijj) > 0. 
A problem with our approach is that we expect unstable modes to have complex u in 
general, with the real and imaginary parts of oj related in some way. In other words, 
for given m, k and a, the complex quantity oj will be a function of the real quantity k 
and hence the real and imaginary parts of u cannot be independent. If we try to follow 
the above strategy for a randomly chosen complex value of u, then this will not satisfy 
the required relation between its real and imaginary parts, and therefore our numerical 
method will not output a real value of k. In order to locate where k vanishes we would 
have to scan over both a and, say, the real part of oj. It would be difficult to do this with 
high accuracy. 

We shall circumvent this problem by restricting attention to modes with m = 0, i.e. 
modes preserving the rotational symmetry of the black hole that follows from the rigidity 
theorem. There are reasons to expect that unstable modes with m = will have purely 
imaginary oj: oj = iT, F > 0, and this is confirmed by our results. Ref. |108j has obtained 
numerically non-uniform rotating black brane solutions that do indeed bifurcate from the 
uniform branch (based on cohomogeneity-1 MP solutions) at a point corresponding to a 
stationary perturbation. Hence stationary perturbations do indeed exist. We believe that 
the reason for this is that unstable modes will have Re(u;) = if they are invariant under 



the rotational symmetry of the black hole predicted by the theorems of Refs. 55 ,66], i.e. 
the symmetry generated by the Killing field r^jmj, where are the rotational Killing 
vector fields and the associated angular velocities of the horizon. We do not have a 
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proof of this, but our results, and the results of Refs. 72|l08 , indicate that it is true. In the 
limit /c — )• A;*, this gives a stationary threshold mode that preserves this symmetry. This 
symmetry is a necessary condition for the threshold mode to correspond to a bifurcation 
into a new family of non-uniform rotating black brane solutions, since presumably this 
new family should respect the theorems of Refs. [55[|66| (although, strictly speaking, these 
theorems apply only to black holes, not black branes). 

In summary, we expect an instability of the black branes (F > 0) to appear for 
wavenumbers \k\ < k^,. If we restrict attention to modes invariant under the symmetry 
generated by ^lirm then unstable modes will have Re(a;) = 0, and the threshold unsta- 
ble mode, with k = k^,, will be stationary and invariant under the same symmetry. For 
cohomogeneity-1 black holes, rijmj is proportional to d/d'ip, so modes invariant under the 
symmetry generated by ^irrii must have m = 0, which is why we set m = above. Notice 



that only m = modes are required to obey the ultraspinning conjecture in Section 6.2 



since the conjecture only applies to perturbations which do not break the rotational sym- 
metry generated by VLirrii. 

In summary, we shall set m = 0, r^f = 1 and, for given (F, k, a) we shall determine 
the possible eigenvalues — A;^. Then we vary a until the eigenvalue vanishes. We have then 
found a black hole that admits an unstable zero- mode with the given values of F and k. 



8.3.2 Results for D = 5 

Our expectation is that, for small a, the black hole will be classically stable but the 



associated black branes will suffer a Gregory-Laflamme instability 100 . Therefore, for a 
range of F, there should exist real solutions for k (corresponding to unstable perturbations 
of the branes) but k will never vanish for F > 0, so the black hole is stable. For a static 
brane, the Gregory-Laflamme instability is an s-wave perturbation of the transverse black 
hole, which for us translates into a k = perturbation. Hence, for a rotating brane, it is 
natural to expect this instability also to have «; = 0. 

This is indeed what we find. The left plot in Fig. 8.1 shows our result for —k'^ for 
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Figure 8.1: Results in D = 5, k = 0. We represent this mode for fixed values of Fr^v/ (first 
graph) and a/rM (second graph). This unstable mode of the branes corresponds to the 
well-known Gregory-Laflamme mode. 



given a and F, with k. = 0. The plot with F = corresponds to a stationary perturbation 
of the black branes. This is the "threshold unstable mode" at the critical wavelength 
beyond which the black branes are unstable. Note that the curves exist for all values of 
a, i.e. the Gregory-Laflamme instability is always present, it does not "switch off' as a 
increases. The upper curves do not extend to /c = so there is no indication of any black 
hole instability. 



In the right plot of Fig. 8.1 we give a more familiar plot of F against k for different 
values of a. For each value of a, we have a curve that takes the usual Gregory-Laflamme 
form. The maximum value of F is 10 — 20% of rj\,f and increases with increasing a. Fur- 
thermore, the range of k for which there exists an instability increases, i.e. the instability 
persists down to shorter wavelengths as a increases. Hence rotation makes the branes 
more unstable. As usual, F — t- as /c — t- and as A; — )■ A;* > 0. The mode with T = k = 



has the usual interpretation of a gauge mode [l03j. The mode with F = and k = k^ 
is the threshold unstable mode associated with the onset of instability. This marks the 
bifurcation of a new family of non-uniform rotating black brane solutions. These are the 



solutions constructed in Ref. 1108 in the black string case. They preserve the symmetries 
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Figure 8.2: Results in D = 7, k = 0. The graphs are entirely analogous to the ones in 



Fig. 8.1 



of the cohomogeneity-1 MP black hole but break the translational symmetry along the 
string. 

Note that the slope T/k approaches a common limiting value as /c — t- 0, indepen- 
dently of the value of a. This long- wavelength limiting behaviour is captured by the 
blackfold approach: it follows from Ref. |6^ that T/k ^ l/\/D — 2 as A; — )■ 0. This is 
consistent with our numerical results. 



We find no solution of (8.6) with k = 1. Therefore our results are consistent with 



stability of these black holes, in agreement with the results of Ref. 114 . 



8.3.3 Results for = 7 



For K = 0, we have the Gregory-Laflamme instability shown in Fig. 8.2 This is qualita- 
tively the same as for D = 5. Once again, rotation makes the branes more unstable, the 
behaviour as A; — )• is consistent with T/k ^ \/\J D — 2 independently of a, and there is 
a threshold mode at a critical value of k at which we expect a bifurcation of a new family 



of non-uniform branes analogous to the strings constructed in Ref. 108 
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Figure 8.3: Results in D = 7, k = 1. We represent this mode for fixed values of T tm (first 
graph) and ajru (second graph). It corresponds to a new Gregory-Laflamme instability 
of the rotating black branes, appearing for the numerical value (ai/rA/)"^"™ = 0.5949. In 
the first graph, the vertical line to the left corresponds to the analytical prediction of 
(8.5), (ai/rj\/) = 0.5946, and the interrupted vertical line to the right corresponds to 



extremality. The second graph clearly indicates that the instability of the black branes 
does not extend to an instability of the black hole (/c = 0). 



A new feature of -D = 7 is the existence of an ultraspinning regime. This occurs 



for a > oi, where a\ was defined in equation (8.5). Just as in Ref. ^72j (Chapter [7|), 
we find that a new stationary (F = 0) negative mode of the black hole appears at this 



point. This negative mode has k = 1. Our numerical results are shown in Fig. 8.3 These 
results demonstrate that the black branes have an instability in the k = 1 sector when 
a > ai. This is a new Gregory-Laflamme instability of the black branes, distinct from 
the instability in the k = sector. The plots of F against k have the same qualitative 
shape as for the k = instability except that the slopes of the curves appear to vanish 
(for all a) as ^ 00 Once again there is a threshold unstable mode at a critical value 
of k. Presumably this corresponds to a bifurcation to a new family of non-uniform black 
brane solutions. In addition to breaking the symmetry along the branes, this mode also 
breaks some of the symmetry of the black hole (typically down to that of a generic MP 



^It would be interesting to investigate whether this behaviour can be explained using blackfold methods. 
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black hol^ so this new family has less symmetry than the non-uniform branes associated 
to the threshold unstable mode with n = 0. 

Note that the k = 1 black brane instability coexists with the k = instability. The 
latter is clearly dominant since it has much larger T and the instability exists for a larger 
range of k, i.e. down to shorter wavelengths. 

It is important to note that there is no evidence of any instability of the black hole: 
none of the curves with non-zero T extends to A; = 0. In the limit — )• 0, solutions with 
non-zero T approach a pure gauge mode, just as for k = 0. Additionally, the solution with 
r = does not approach a gauge mode as /c — >• 0. Instead, as anticipated in Chapter [6], it 
corresponds simply to a variation of parameters within the MP family of solutions. 

For D = 7, since the black hole is ultraspinning for a > oi, there is the possibility 
of an ultraspinning instability appearing at a = 02 > ai. However, we find no solution of 



equation (8.6) for k = 2 so our results are consistent with stability oi D = 7 cohomogeneity- 



1 MP black holes for axisymmetric perturbations (m = 0). 
8.3.4 Results for D = 9: Black hole instability 

For K = 0, we have the expected Gregory-Laflamme instability of the black branes. For 
K = 1, we find, as for D = 7, a new Gregory-Lafiamme instability of the black branes that 



appears at a = ai. This is shown in figure Fig. 8.4 



The new feature that appears for = 9 is an instability in the k = 2 sector, which 



appears at a = 02 > ai- This is shown in Fig. 8.5 The left plot shows a new stationary 
(r = 0) negative mode which emerges from a zero- mode at a = a^- We shall prove in 
Appendix 8.B.4| that this zero- mode cannot correspond to a variation of parameters within 



the MP family of solutions. Furthermore, in Appendix 8.B.2 we show that it cannot be 
a gauge mode. 

For a > 02, there is a new instability of the black branes, corresponding to the 



^ This is because k—\ harmonics are in one to one correspondence with Kilhng vector fields of CP^ : 



see Appendix 8. A. 2 
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Figure 8.4: Results in D = 9, k = 1. The graphs are entirely analogous to the ones in 
Fig.lO 



curves with F > in the plot. But there is a qualitative difference between the left plot of 



Fig. 8.5 and our previous plots: the curves with F > now intersect k = 0, i.e. we have 
found perturbations of the black hole that grow exponentially in time, that is, a classical 
instability of black holes with a > a2. This is our main result. 

The onset of instability is indicated by the stationary zero-mode (F = 0, = 0) 
at a = 02- This is analogous to the mode constructed for singly-spinning black holes in 
Ref. (72) (Chapter [7|). Our main achievement here is to demonstrate, for the first time, 
the existence of modes which grow exponentially with time when a > a2. 



The right plot of Fig. |8.5| shows a clear difference from our previous plots. Unlike 
the original Gregory-Laflamme instability, we find that F is maximized at k = rather 
than vanishing there. Hence, for the black branes, the most unstable k = 2 modes are 
those with A: = 0, i.e. those corresponding to the black hole instability. For larger k, the 
black brane instability "switches off' in the same way as the original Gregory-Laflamme 
instability, with a threshold mode at A; = /c* indicating a new family of non-uniform black 
brane solutions. 



Fig. 8.6 presents our result for the instability time scale of the black hole as a function 
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D = 9,K = 2 D = 9,/f = 2 




aAM k Tm 

Figure 8.5: Results in D = 9, k = 2. We represent this mode for fixed values of Tvm 
(first graph) and a/rjv/ (second graph). As opposed to the k = 1 case, the time-dependent 
mode extends all the way to A; = 0. There is not only a new Gregory-Laflamme instability 
of the black branes, but also an instability of the black hole, appearing at a = 02, where 
O'2/i'M = 0.6858 > ai/rM = 0.6300 . In the first graph, the interrupted vertical line to the 
right corresponds to extremality. In the second graph, the curve shrinks to the origin as 
a — ^ 02. 



of its spin. For a > 02, we find that F increases monotonically with a, so extreme black 
holes are the most unstable. 



Finally, we find no solutions of (8.6) with k = 3. 



8.3.5 Rotational symmetries of higher-dimensional black holes 

As explained before, the study of perturbations of higher- dimensional black holes can 
be used to investigate the possible existence of new families of black hole solutions with 



less symmetry than the known solutions. The idea proposed in Ref. 67 is to look for a 
stationary zero-mode of the black hole, which is interpreted as indicating the existence of 
a family of solutions branching off from the known solutions. 

In our case, the stationary zero-modes with k = 1 are uninteresting since these 
correspond to variations within the MP family. However, for D = 9, we found a stationary 
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Figure 8.6: Results in D = 9 for k = 2, in the limit k = 0. We represent the black hole 
instability in a plot of Fr^v/ versus a/vM- The interrupted line corresponds to extremality. 
Numerical error prevents us from extending our results all the way to extremality. 



zero-mode with k = 2 that appears at a = 02, the critical value of a beyond which the 
black hole is unstable. Therefore, we have found evidence for a new family of black hole 
solutions that bifurcates from the MP family at this pointj^ 

How much symmetry do these new solutions have? This can be inferred from the 
symmetry of the k = 2 harmonics on CP^ . There will, of course, be a family of degenerate 
scalar harmonics with k = 2. Some of these will preserve some of the symmetry of CP^ 



whereas others break it completely (this is proved in Appendix 8.A.3[ )p| For a mode of 
the latter type, the associated metric perturbation will break completely the SU{N + 1) 
subgroup of the R x U{N +1) isometry group of the background metric. It preserves only a 
RxU{l) subgroup corresponding to time-translation invariance and invariance under the 
rotations generated by d/dij). Hence the corresponding family of new black hole solutions 



We expect that this new family will have unequal angular momenta in general. However, this cannot 
be seen from our results since k = 2 modes do not change the mass or angular momenta at the linearised 



level (see Appendix 8.B.4I. A second order calculation would be required to determine these changes. 

^ It is helpful to think about the case of CP^ = S^, for which k — £, the total angular momentum 
quantum number. Scalar harmonics are labelled by £ and m. Certain modes with £ = 2 (i.e. quadrupole 
modes) may preserve some symmetry (e.g. if m = then they are axisymmetric) but generically they 
break all of the continuous symmetries of S^. 
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will possess just a single rotational symmetry. 

Very recently, Ref. [65^ has constructed approximate black ring solutions with just 
a single rotational symmetry. Our results are the first evidence for the existence of new 
black hole solutions with a single rotational symmetry and horizons of spherical topology. 

Note that if one used a k = 2 harmonic that does preserve some of the symmetry of 
CP^ then presumably this would give rise to a different family of new black hole solutions, 
with more than one rotational symmetry. Therefore, assuming that each stationary zero- 
mode corresponds to a new non-linear stationary black hole solution, there must exist 
several new black hole solutions that bifurcate from the MP family at the same point, and 
these different solutions have different numbers of rotational symmetries. So how many 
new solutions are there? 

One way of addressing this question is to determine the number of parameters in the 
most general k = 2 harmonic. If we take D = 9 then k = 2 harmonics correspond to the 
[2,0,2] representation of SU{4:), which is 84-dimensional. Hence the most general k, = 2 
harmonic is labelled by 84 parameters. Some such harmonics are related by acting with 
SU{4:), i.e. by rotations of the background spacetime. However, since Sf7(4) has dimension 
15, this can eliminate only 15 parameters, leaving 84 — 15 = 69 parameters that cannot 
be eliminated by rotations of the background. So, up to rotations of the background, we 
have a family of stationary zero-modes with 69 parameters, and presumably a family of 
new black holes with 70 parameters, the extra parameter being the mass (or r^)- This 
is considerably more parameters than the 5 that are required to specify the D = 9 MP 
solution! 

8.3.6 Expectations from the ultraspinning conjecture 



The predictions of the ultraspinning conjecture for MP black holes were discussed in 
Section 6.2.2 In the equal spins case analysed here, the reduced Hessian 
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was shown to possess an eigenvalue which is always positive (with eigenvector Vi = Vyi), 
and N degenerate eigenvalues (with eigenvectors satisfying Vi = 0) which change from 
positive to negative as the ultraspinning surface is crossed, with no further changes of 
sign at larger angular momenta. Note that only the positive eigenvalue corresponds to a 
variation which preserves the equality of the angular momenta. We expect that new ther- 
modynamic negative modes will emerge only at the unique value of the angular momentum 
corresponding to the ultraspinning surface, that there will be precisely N of these, and 
that they will break some of the symmetries of the background. 

The ultraspinning surface corresponds to a = ai (with vm = 1) and we found that 
new stationary (F = 0) negative modes do indeed emerge at this point in D = 7, 9. These 
modes correspond to k = 1 harmonics. Since these are thermodynamic negative modes, 
we know that the zero-mode at a = oi must be simply a variation of parameters within 
the MP family. 

To see that there are precisely negative modes emerging at a = oi, we use the 
fact that K = 1 harmonics are in 1-1 correspondence with Killing vector fields on CP^ , 



so there are (A^ -|- 1)^ — 1 such harmonics (see Appendix 8. A. 2). However, some of these 
are related by rotations of CP^ , so we need to determine how many parameters can be 
eliminated by rotations. The counting is the same as for SU{N + 1) gauge theory with 
an adjoint Higgs field. Generically this breaks SU{N + 1) to U{1)^ so we are left with A^ 
parameters]^ Hence there are A^ independent negative modes that emerge at a = ai, in 
agreement with the prediction from thermodynamics. 

Our numerical results confirm that the stationary negative mode that emerges at 
a = ai does indeed correspond to the onset of a new instability of the black branes in the 
K = 1 sector, in agreement with the refinement of the Gubser-Mitra conjecture discussed 
in Section [6]2l 



For D > 5, cohomogeneity-1 black holes are ultraspinning for a > ai and hence 
might exhibit an instability of the form anticipated in Section [6. 2[ Moreover, as explained 



^Thanks to David Tong for this argument. 
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at the end of Section 8.2 , as we increase D, there is more "space" between the ultraspinning 
surface and the surface of extremahty for such black holes. Therefore the likelihood of an 
instability might be expected to increase with D. This is in agreement with our numerical 
results, which show no sign of any instability for D = 7 but confirm that an instability is 
present for D = 9. Note that this instability does indeed occur inside the ultraspinning 
region. 

The onset of instability is associated with the appearance of a new stationary zero- 
mode (at a = 02). This cannot correspond to a variation of parameters within the MP 



family (we show in Appendix 8.B.4 that k = 2 modes cannot change the asymptotic 



charges). Furthermore, we prove in Appendix 8.B.2 that it cannot be a pure gauge mode. 
This zero-mode is continuously connected to a stationary (T = 0) negative mode that 
exists for a > 02. This is an example of a non-thermodynamic negative mode, i.e. one 
which is not associated to a local thermodynamic instability. The same behaviour was 
observed in Ref. |72| (Chapter [7]), i.e. the onset of a classical instability of the black hole 
is associated with the appearance of a new stationary negative mode. 



Ref. [72] found that further non-thermodynamic negative modes appear as the spin 
of the black hole is increased still further. In the present case, extremality imposes an 
upper bound on the spin of the black hole and we do not find any further negative modes 
beyond the ones associated with the instability in the k = 2 sector. However, we believe 
that, for larger D, as well as an instability in the k = 2 sector there will be further negative 
modes in sectors with larger k. These new negative modes will be associated with new 
instabilities of the black hole in sectors with larger k. The stationary zero-modes marking 
the onset of these instabilities will provide evidence for the bifurcation of new families of 
black hole solutions, involving a large number of parameters, and generically with just one 
rotational symmetry. 
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8.4 Scalar perturbations and CP^ harmonics 
8.4.1 Introduction 

The rest of this Chapter is devoted to explaining the technical details of our work. We 
shall start by explaining the decomposition of metric perturbations into harmonics on 
CP^. 

Metric perturbations can be decomposed into scalar, vector and tensor types ac- 
cording to how they transform under isometrics of CP^ . Pertubations of different type 
must decouple from each other. The decomposition is explained (for a different problem) 



in Ref. 1801. Tensor perturbations are the simplest, these were discussed in Ref. 120). 



We are interested in scalar perturbations, for which the perturbation can be expanded in 
scalar harmonics on CP^ . 

We shall assume that our perturbation has been Fourier decomposed as in equation 



7). All of our numerical results assume m = 0. However, for the sake of completeness, 
we shall derive equations that are valid for non-zero m. In order to do this, we must 
address a subtlety (already encountered in Ref. [120| ), that such a perturbation couples 



with charge m to the 1-form Aa on CP^ defined in Section 8.2 Hence we must consider 
charged scalar harmonics on CP^ . First we shall describe these harmonics and then 
explain how to construct gravitational perturbations from them. 



8.4.2 Charged scalar-derived harmonics in CP 



We define the gauge-covariant derivative acting on a charge-m tensor field on CP as 

Va = Va-imAa, (8.9) 
where V is the metric covariant derivative on CP^ . 
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Scalar s 

Charged scalar fields on CP^ can be expanded in terms of charged scalar harmonics 
defined by 

(p2 + A)Y = 0. (8.10) 



These were studied in fl81| (see summary in Appendix 8.A), where it is found that 



\ = i{e + 2N)-m^ , i = 2K + \m\, (8.11) 

with K = 0,1,2,.... The modulus sign guarantees that the eigenvalue is the same for 
positive and negative charges, the corresponding eigenfunctions being related by complex 
conjugation. 

Notice that the presence of the 'gauge field' A leads to 

[Va,Vb]Y = -i2mJabY. (8.12) 

Scalar-derived 1-forms 



Given a scalar harmonic Y, we can define 



10 



Ya = -^VaY, (8.13) 



which transforms as a charged 1-form on CP^ . This can be decomposed into its (1, 0) and 
(0, 1) parts using the complex structure on CP^ . Denote these by Y+ and Y^^ respectively, 
where 

J,''Y± = TiY±. (8.14) 



'^"Note that A = if, and only if, k = m = 0, in which case Y is uncharged and constant. In this case 
there are no scalar-derived vectors nor scalar-derived tensors. 
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We shall refer to as scalar-derived 1-form harmonics. We find that they satisfy 

p2Y± = - [A - 2(iV + 1) ^ 4 m] Y^ (8.15) 



P'^Y^ = — (A T 2 m iV) Y . (8.16) 



and 

2VX 

We shall make use of the result that Killing vectors of CP^ are in one-to-one corre- 
spondence with uncharged (m = 0) scalar harmonics with k = 1 (see e.g. [181j and 



our Appendix 8.A). Given such a harmonic Y, the corresponding Killing vector field is 
-i(Y+-Y-). 



Scalar-derived tensors 



Following Ref. [180], we decompose a symmetric tensor Yab into its Hermitian (or (1, 1)) 
and anti-Hermitian components according to the eigenvalue of the map 

iJVab = Ja'Jb'^^cd . (8.17) 

If the eigenvalue is +1 the corresponding eigentensor is called Hermitian, and if it is —1 the 
eigentensor is called anti-Hermitian. In the anti-Hermitian case, we can further distinguish 
between the (2,0) and (0,2) components of Y^f, , which are defined by J^'^Y^b = =Fi^a6 
with the upper and lower signs for the (2, 0) and (0, 2) components respectively. 

The following quantities form a basis for anti-hermitian scalar-derived tensors 

= K^t) ' ^.T = -^la^b) ■ (8-18) 

Y^ denotes the scalar-derived 1-form harmonics of the previous Section, and denotes 
the projection of Da onto its (1, 0) and (0, 1) components. Notice that the correspondence 
between m = 0, k = 1 scalar harmonics and Killing vector fields on CP^ implies that 



This follows from the scalar part of equation (39) of Ref. 180 
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vanish for such harmonics. 

Hermitian scalar-derived tensors can be written in terms of a trace and a traceless 
part, for which the foUowing quantities give a basis 



12 



9abY, Y+- = V+Y- + V-Y+ - ^ gat (V-Y). (8.19) 
These tensor harmonics satisfy 



P2Y±± = - [A - 4(iV + 3) T 8 m] Y±± , 



^.20) 



with 



5.211 



, ^ _A-4(Ar + l)T2,„(iV + 2) 

^Xa~ = -^^[(A + 2miV)Y+ + (A-2miV)Y-] . 

8.4.3 Decomposition of perturbations in scalar-derived harmonics 

Let us now consider the perturbations of the full spacetime metric. We introduce the 
orthonormal basis 



,(0) 



e^"' = fdt, e^^^=gdr, e^^'^ = h {d^j + A - n dt) , e^'^=re^'\ (8.22) 



where e'^*-' is the tetrad of the CP^ manifold. The dual basis is then 

6(0) = \{dt + ^ d^) , 6(1) = -dr, 6(2) = 7- ' Hi) = " ~ Hi)) ' 

J g 11 r 

(8.23) 

Take 6^"^^ = {e*^"-*, e^^^ e*^^-*} and a coordinate basis dx'^ on CP^ . The components Hab of 
the metric perturbation transform as scalars under isometries of CP^ and can therefore be 
decomposed using scalar harmonics on CP^ . Similarly, since we are restricting attention 



^'^ In Ref. [l80| , hermitian tensors were converted into (1, l)-forms by contracting with J°'b- The two 
quantities written here correspond to terms of the form JY and (the primitive part of) dd'^Y in equation 



(47) of Ref. 180 
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to scalar-type perturbations, components of the form hAa and hat can be decomposed 
using scalar-derived 1-forms and scalar-derived tensors on CP^ : 

hAB = Jab Y , 

hAa = r{fiYt + fA^a) , (8.24) 
/^a6 = {H+^ Y++ + H- Y-,- + H+- Y+-) + r2 9ai>Y , 



where = {M^^, X^, Z^}, and the functions multiplying the harmonics depend only on 
{t,r,ip) and not on the coordinates of CP^ . The real spacetime metric perturbation is 
given by Re [h^i^). Since dt and 9^ are Killing vectors of the background solution, we will 
Fourier expand all of these functions in t and if), i.e. we assume a dependence e~"^*+''"'^. 
It remains to determine the dependence of these functions on r. The stability problem 
will thus be reduced to a system of linear ordinary differential equations. 



8.4.4 Boundary conditions 

The metric perturbations must be regular on the future event horizon T-L^ . This boundary 
condition can be imposed by considering a basis which is regular on H.^ , since the com- 
ponents of the perturbation in that basis must be regular. Let us change to the ingoing 
Eddington-Finkelstein coordinates that are regular at T-L^: 

dt ^ dv - J dr , dtp ^ dtp -^dr , (8.25) 

and consider the basis {dv, dr, dip + A — Qdv, dx^}. Denote the components of the metric 
perturbation with respect to this new basis with a bar (e.g. /gg, 1^^). Our bound- 
ary condition is that these components should be smooth functions of {v, r, cj), x") at the 
horizon]^ 

^^In other words, we demand that the tensor field h^i, should be regular at . This is stronger than the 
statement that the metric perturbation should be regular at the horizon, e.g. it excludes the possibility that 
hfii, is singular at H'^ in a certain gauge but can be made regular by a gauge transformation. For example, 
we show in Appendix |8.B.3| that, in the traceless-transverse gauge, a perturbation with uj = m — that 
satisfies our boundary condition cannot change the temperature or angular velocity of the black hole. It 
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For this class of black holes, the horizon is located at the largest real root r = r+ 

of A = gir)^"^. For a non-extreme black hole, near the horizon, A(r) = A'(r4_)(r — r_|_) + 
0[{r — r^)^], with A'(r+) > 0. Using the relation /(r) = r/ {g{r) h(r)), we find that, near 
the horizon, the metric components in the original basis are related to the components in 
the new basis by 



hjr+f /oo f r h{r+) 

2 Mt ^ 1 /Ol — /oo — 

r^A (r+) r — r+ r+ 



/oo-2/oi + /ii« A'(r+)/„ (r-r+), 

f 1 /02 . , ^^A^ ^ 

" rWA-(r,) 7^ ' " ^^^^^ ' ^^-^'^ 

The functions ++ , , i?"' , Hl associated with the components of the metric per- 
turbation on CP^ are the same in the two bases. 

Since the components in the new basis should be regular at the horizon, the above 
expressions give us boundary conditions on the behaviour of the components in the old 
basis. In imposing these boundary conditions, it is important to remember that, near the 
horizon, 

-iusv+'mnp ^ —ioji+imip I + 



\ r+ J 

where a = h{rj^)/{rj^/S.'{r^)) is positive (for non-extreme black holes). Hence, for example, 
the radial dependence of /oo near the horizon must be 

/oo oc + F(r), (8.28) 



where F{r) is smooth at r = r+. 



follows that a variation in the parameters of the MP solution that does change Th or Q.h will not give a 
perturbation hjj,v that is regular at the horizon in this gauge. 



8.5. THE EIGENVALUE PROBLEM 



155 



When numerically solving the stability equations, it will be necessary to work with 
the combinations that maximize the information on the boundary conditions. For instance, 
one should work with /02 and /12 — fo2, instead of considering only the leading behaviour 
of /02 and /12, otherwise the information that /i2(r+) — /o2(?'+) = is lost. 

As for the behaviour of the perturbations at spatial infinity r — t- 00, we are inter- 
ested in boundary conditions that preserve the asymptotic flatness of the spacetime. For 



perturbations of the black branes, the equations of motion (8.6) then imply that all the 
functions vanish exponentially for large r. 



8.5 The eigenvalue problem 



The ansatz for the metric perturbation hf^^, is given by Eq. (8.24). Ref. [68j lists in a 



long Appendix the components of the Lichnerowicz eigenvalue equation (8.6) in the tetrad 



basis (8.22). These consist of sixteen coupled second order ordinary differential equations. 



each one being second order only in one of the perturbation functions. However, six 
of these functions can be solved for in terms of the ten remaining functions and their 
first derivatives when we impose the traceless-transverse (TT) gauge conditions, listed in 
another Appendix of [68]. The procedure is analogous to the one applied in Chapters |4] 
and [3 



Notice that the TT conditions completely fix the gauge in Eq. (8.6) when k > 0, 
since the action of the Lichnerowicz operator on a gauge mode is trivial, AiV(^^,y) = 0. As 
for k = 0, there are two distinct cases. The first is the limit A: — )• for which F — )• 0, as can 



be seen on the right plots of Figs. 8.1 , 8.2 , 8.3 and 8.4 The limiting perturbation = is 



an unphysical pure gauge mode, as happens in the original Gregory-Lafiamme case 103 



The second is the much more interesting stationary perturbation k^, = marking the onset 
of a new Gregory-Lafiamme instability when the rotation increases. In the left plots of 



Figs. 8.3 8.4 and 8.5, this is the threshold mode of the curve F = 0. In the right plots 
of the same Figures, it would correspond to the graph squeezing into the origin for a 
critical value of the rotation. These stationary perturbations are physical since there is 
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no gauge ambiguity. The TT conditions require that any gauge vector is a harmonic 



1-form, satisfying V^^^ = and V^Vp^^ = 0. In Appendix 8.B.2, we show that no regular 
harmonic 1-forms exist in this case: they lead to pure gauge metric perturbations that 
diverge either at the boundary r = r+ or at infinity r — t- oo. A proof along the same 
lines, but much more cumbersome, can be given for the modes which represent the actual 
instability of the black holes, i.e. the exponential growth of the perturbations with time 
(fe = and r > 0). 



The tetrad basis (8.22) is very convenient in the explicit derivation of the TT gauge 
conditions and the Lichnerowicz equations. However, for the actual implementation of the 
numerical problem, it is convenient to choose perturbation functions that make the final 
equations more amenable to numerics, e.g. it is useful to avoid using expressions involving 
square roots. It is also helpful to define combinations of the original perturbation functions 
which can be solved for algebraically through the TT gauge conditions. Both features are 
respected if we consider the perturbations in a related basis such that: 

foo = /oo f -2fo2fhn + f22 h^n^ , hi = foif9-fi2ghn, 

f02 = f02fh- /22 f ^, fll = fll 5^ , fl2 = fl2gh, ^22 = f22 h'^ , 

fo = r {{W+ + W-)f- {Z+ + Z')hn) , 

fo = -i 1 r {{W+ -W-)f-{Z+-Z-)hn), (8.29) 
h = -lrg{X+ + X-), ~h = -ilrgiX+ -X-), 
U = -\rh{Z+ + Z-), }^ = -i^rh{Z+-Z~), 



P = I {2H+- - H++ - H~-) , Q = i i - H-) , 
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U = I {2H+- + H++ + H—) , V = HL + i^ {H+- - \ {H++ + H—)) . 

Now, we solve for six of these functions (foo, foi h-> Q) ^) ™ terms of the ten remaining 
functions and their first derivatives by imposing the TT gauge conditions. Upon this 
substitution, the corresponding second order Lichnerowicz equations will become third 
order. The ten equations which are second order in foi, fo2, fii, fi2) f22, foi fi) fii 
will remain second order. They constitute the system of equations to be solved numerically. 
A non-trivial consistency check on the gauge choice procedure is that the ten final second 
order equations must solve the six third order equations (e.g. a third order equation is a 
derivative of a second order one). We verified explicitly that this is the case. 

The final system will be solved using a spectral numerical method, briefly described 
in the Appendix at the end of this thesis. The application of the method is simpler for 
Dirichlet boundary conditions. We then consider the following perturbation functions: 



Q3 
Q5 



(l 
(l 
(l 
(l 
(l 

(l 



\ mSlfl- )+2 



11 



\ iQ!(w— mil^f )+2 



\ ia{(jj—mQjj)+l 



01 



fl, 



P. 







r+ 


q2 = 


r 






!± 


qA = 


r 


qe = 


0- 


r+ 
r 



22 



\ mH/f )+2 



fl, 



7'_l_ \ ici(a;— mSl^f ) + l _ 



f02 + f22 + 



qio 



r 

r / 
r 



r+ 
a Q.H 

a 



roi 



+ ^Hh- — {l 

a 



1 

r+ 
r 



) [^°^ 
fii 



!± 

a 



11 



(8.30) 



which vanish linearly at the horizon location r = r_|_. This behaviour can be verified in 



the expressions (8.26) and (8.27). The particular combinations chosen for gg, ® and gio 



encode the total information about the boundary conditions imposed by regularity, as 



argued in the end of Section 8.4.4 For the numerical implementation, it is convenient to 
use the variable 

y = i-— (8.31) 
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instead of the radial coordinate r, since y is dimensionless and bounded, < y < 1 . The 
functions represented above vanish at infinity r = oo {y = 1) since the large r behaviour 



of the solutions of (8.6) is exponential, e , and regularity at infinity imposes the minus 
sign. The k = case will be obtained as the limit A; — )• 0. 

The system of ten second order ordinary differential equations is ready to be solved 



numerically. The results were presented in Section 8.3 



8. A Appendix: The geometry of CP^ 
8.A.1 The Fubini-Study construction 

We review here the Fubini-Study construction of the Einstein-Kahler metric and Kahler 



potential on CP^ 181 This construction allows us to iteratively generate the CP^ 
metric and potential from the knowledge of the metric and potential of CP^~^ . 

Take the C^^"^ manifold with complex coordinates and flat metric 

dslN+2 = dZ^ dZA , (8.32) 

where the index A runs as ^ = (0,a), with 1 < a < A^. Introduce N inhomogeneous 
coordinates = Z^/Z^ in the patch where Z^ ^ , such that 

Z° = e'^|Z°|, Z" = Z'^C = Rnu'^ , 

Z^ZA = r\ f = ! + f = i + R%. (8.33) 

Furthermore, introduce a new set of (A^ — 1) inhomogeneous coordinates {0 < i < N — 1) 
such that 

= e''^^/^\u^\, u' = u^v' with w"n" = l. (8.34) 



^*We use the coordinates {Rn,'^n} that are related to the coordinates {Ci''"} of 181 through the 
coordinate transformation sin i = R%/{1 + Rjf) and ? = *iv/2. 
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The flat metric on can then be written as 



ds2N+2 — dr +r dQ.2N+ii where dQ.2]s[j^i = [dr + A(^^^) + (8.35) 
is the metric on the unit sphere S"^^^^, and is the unit CP^ metric. Written in this 



way we see that 5'^^+^ is a Hopf fibration of over CP^ . In ( |8^ , A(^n) is the CP^ 



Kahler potential. Explicitly, the CP^ metric and Kahler potential are given by {Rn > 
and < ^^AT < 47r) 

dR^T 1 i? ^r / , ^ „ , N 2 R^r, 



d^% = g^.dx-dx" = ^ + - I (d*^ + 2^(^.,))^ + — ^ dE^_i , 



1 ii- 



in terms of the Fubini-Study metric, and Kahler potential, ^(^v-i)) the unit 



dT.l^^ = f^\dv' dv^ - f^Wv' dv'W 

A^N-i) = J i /iv-1 '^^^ - rf^O > /w-i = 1 + • (8.37) 



By definition, the Kahler form on CP^ , Jn = ^dA]\! , is covariantly conserved, 
Va J^*^^ = , and satisfies Jj'Jtc = -Qac ■ 

The lesson from this analysis is that starting from the CP^ fields we can iteratively 
construct the CP^ geometry as well as the complex coordinates that define the em- 
bedding of CP^ in C^^^. CP^ is isomorphic to the 2-sphere S^, its metric and Kahler 
potential being given by 

dS? = ^ (d6i2 + sin^ Odcj?) and ^(i) = ^ cos 6* # . (8.38) 

Examples of the embedding in C^^^ may be elucidative. For CP^, parameterized by 
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{6, (p, i?2, ^2), the map is given by 

{Z', Z\ Z^) = -l£= f 1, R, cos I \ (*^+^) , sin { \ , (8.39) 

while for CP'^, parameterized by (0, (/>, i?2, ^2,-^3, ^3), the map is 

(Z0,Z\Z2,Z3) = 

'l,^^i?2Cos'e^i(*^+^),^^i?2sin^e^|(*-^)V (8.40) 



In order to reproduce the results in Section 8.4.2 , it is useful to recall that, for CP^ , 

Rabcd = 9ac9bd — (jadElbc + JacJbd — JadJbc + "^JabJcd ■ (8-41) 

8. A. 2 Scalar harmonics and Killing vectors 

We review here the systematic way to construct all scalar harmonics and Killing vector 



fields on CP 



N 



181 



The isometry group of CP^ is SU[N + 1). Let Ta^-Aj,'^^'"^^ be a 
constant Hermitian SU{N + 1) tensor, which is symmetric in the index set {^i, . . . , Ap} 
and the index set {Bi, . . . , Bq}, and traceless in any contraction between an Ai and a Bi 
index. This defines the {p,q) representation of SU{N + 1). The charged scalar harmonics 
are then given by 

Y = r^,...^/i-^' • • • Z^^ Zb,--- Zb, , (8.42) 



and satisfy the Laplacian (8.10) for A = 2[2pq + N{p + q)]. We have k = max{p, q} and 
m = p — q. Uncharged scalar harmonics have k = p = q and A = 4k(k + N). 

The Killing vectors on an Einstein-Kahler space can be constructed from the un- 
charged scalar harmonics with k = 1 , which have eigenvalue A = 4(1 + A^). Indeed all the 
Killing vectors ^(j) of CP^ are generated by the relation 



jab o ^771=0 



(8.43) 
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where i = 1, . . . , A^(A^ + 2) . That is, setting to zero all but one of the constant components 
of the arbitrary Hermitian traceless tensor T^P we get a k = 1, m = scalar harmonic on 



CP^ through (|8.42l). Repeating the exercise for all possible combinations, we generate the 



A^(A^+2) uncharged scalar harmonics Y^^^^-^^ , and the associated (A^+1) — 1 = A^(A^+2) 



Killing vectors through (8.43). There are N linearly independent Killing vectors which 
commute with all the others, thus generating the Cartan subgroup C/(l)^ of SU{N + 1). 



8. A. 3 Symmetries of k = 2 harmonics 

The symmetry group SU{N + 1) is broken, at least partially, by any linear perturbation 
hab satisfying C^hab 7^ , where ^ is one of the N{N + 2) Kilhng vectors of CP^. For 
scalar type perturbations, if C^Y ^ for some ^ then the symmetry associated with is 
broken by the perturbation. 

We explained above how to construct the Killing vectors from the uncharged k = 1 
scalar harmonics. Consider now the most general linear combination of Killing vectors 
K = J2i Cj with i = 1, . . . , A^(A^ + 2) . The entire symmetry group SU {N + 1) is 
broken by hab if the only solution to CkY = is Cj = for all i, i.e. K = 0. There are 
(uncharged) k = 2 harmonics on C-P^ {D = 9) for which this is true. For reference, we 
present here a particular example in the coordinate system used above: a family of m = 0, 
K = 2 harmonics with three non-zero continuous parameters, (32 and /33, 



Y' 



m=0 



/3i 



R3R2 



1 + cos 6 



RiRi{l + cos( 
2(1 + R^,) 



62 



RlR2y/l-cose 
{l + Rlf{l + Rl) 



2{l + Rl) 



62 



2 



R^e 



{l + RlY{l + Rl) 



^.44) 



If one of Pi, (32, (33 vanishes then this is still a k = 2 harmonic but it preserves some 
symmetry. 
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Appendix: Properties of stationary axisymmetric modes 



In the main body of the Chapter, we have presented our numerical results for general 
axisymmetric time dependent scalar perturbations. Our numerical code has a continuous 
limit as a; = ir — 7- 0. Therefore this particular case was already included in our discussion. 
However, in this appendix we want to have a closer look at stationary axisymmetric 
modes. The reasons are: (i) we developed an independent code for this particular case 
which confirms the results from our general code with time dependence; (ii) we proved 
that the stationary zero-modes (with uj = and k = = 0) are not pure gauge modes; 
(iii) we confirmed that our stationary perturbations preserve the angular velocity and 
temperature of the background geometry; (iv) we determined which stationary zero-mode 
perturbations can change the mass and angular momenta of the background solution; 
finally (v) we want to give special attention to the stationary modes, and not just to the 
time dependent instability, since these may indicate bifurcation points to new branches of 
black hole solutions. 



8.B.1 Stationary perturbations sub-sector 



As described in Section 8.5, the stability problem of axisymmetric perturbations with time 
dependence consists of a system of 16 Lichnerowicz eigenvalue equations for 16 unknown 
functions. Choosing the TT gauge reduces the problem to a system of 6 TT gauge condi- 
tions and 10 Lichnerowicz equations. The procedure is consistent because, as explained in 
that Section, the latter 10 equations automatically imply, through the 6 gauge conditions, 
that the other 6 Lichnerowicz eigenvalue equations are satisfied. 

When we consider the axisymmetric stationary sub-sector of the perturbations, i.e. 
m = 0, w = ir = Owe find that the initial system of 16 Lichnerowicz equations decouples 
into a subsystem of 10 equations for 10 functions and another subsystem with 6 equa- 
tions involving only the remaining 6 perturbation functions. Moreover, we can use the 
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stationarity condition, w = iT = 0, to further simplify our system of equations. To see 
how this is accomphshed, let us introduce the harmonics associated with the time and 
azimuthal Killing directions as S = e-*'^*e*'"'^. We can then decompose the perturbations 
according to how they transform under the {t, ip} Killing isometries. For example, the 
scalar-derived vector perturbations satisfy hj^i ~ /a d^S for b = t,ip and the index A 
running over the radial and CP^ coordinates (the bar in b denotes that the 1-form basis 
is really {dt,dij) + A — Vtdt}). These perturbations, hAt and /i^^, must then vanish when 
we set = and m = 0. This amounts to requiring that foi = fi2 = fo = f2 = 0. The 
original 10 time dependent Lichnerowicz equations then imply that fi = when we set 
a; = 0. Similarly, the original 6 TT gauge conditions imply that Q = 0. We are then led 
to the following conditions 



foi = fl2 = fo = f2 = fl = Q = , 



^.45) 



or, using the map (8.29): /oi = /12 = 0, 
H— = H++. 



-W+, X- 



x+, z- 



and 



With (8.45) the original system reduces to a closed system of 3 TT gauge condi- 



tions and 7 Lichnerowicz equations. In this axisymmetric stationary case, the boundary 



conditions at the horizon (8.26) reduce to 



/oo - 
/02 - 



/oi> /oo + /ii ~A'(r+)/n (r-r+) 



h{r+) 



/i2\/r - r+ , /22 
W'^^r-r+, X 



1 



z 



± 



1 



r+h{r+) 



^.46) 



where we used = W (r — r-| 



and one further has W 



-W'^, X" = X+ and 



-Z~^. It is important to emphasize that (8.45) already encode the information 



that the perturbations are in the TT gauge, and that the boundary conditions (8.46) are 
compatible with the TT gauge. 
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We have done an explicit search of the stationary modes using only the subsystem 



of 3 TT gauge conditions and 7 Lichnerowicz equations described above subject to (8.46). 
We recover independently the same results that we obtain when we set w = in our time 
dependent code. 



8.B.2 Stationary zero-modes are not pure gauge 



As discussed in Section |8.5[ perturbations with k > in the TT gauge all have the gauge 
freedom fixed. However, TT perturbations with = 7^ A;* and F = are pure gauge 
modes |103| . In this subsection, we want to confirm that our stationary axisymmetric 
zero-modes with k = k^, = cannot be pure gauge modes. Given that for any residual 
gauge freedom the gauge parameter would be constrained to be a harmonic 1-form, we 
will prove that there is no regular harmonic 1-form that could generate our perturbations. 



Consider the effect of a scalar gauge transformation on the metric perturbations. 
The most general scalar-type gauge parameter can be decomposed as 

e = 6"^'^*+''"^ [Co (r ) Y e(°) + 6 (r ) Y e(i) + 6 (r ) Y e(2) + r {r)Yt + C {r)^a ) dx""] . (8.47) 

Under a gauge transformation, 



^.48) 



the tetrad components of the metric perturbations transform as 



/oo — ^ /( 



00 



/ fg 



Joi Joi + - 1 ^ - y l4o J 4i ) 



J02 — ^ 702 + -r- t,0 
n 



?i 7 ?2 



fg 



f 



Jii — Jii + - 7^ ?i , 
g or 



J 12 -5- J12 7 — 4o + ^ ?i + - ^ r ?2 

fg h g\dr h' 



/22 — ^ /22 + 2 



hg n 
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W+ - 
X^ ^ x+ - 



— Co + ^ 7 t 

r J 



r g \ or r 



r \ h r"^ 



X- X- 
z- ^ z- - 



"v/A i(a; — m il) 
— 4o H 7 ? 



r 9 \or r J 



r r 



-ii + ^=(C+(A- 2miV) + ^-(A + 2miV)) 

5 4iV^/A^ ' 



2VA 



3.49) 



Our stationary axisymmetric perturbations must satisfy (8.45) which requires that a po- 
tentially dangerous gauge parameter ^ must satisfy 



Ur) = 6(r) = , and T (0 = t{r) . 



^.50) 



We now prove that a parameter ^ obeying these conditions cannot generate a pure gauge 
metric perturbation that is regular. By regularity we mean that the gauge transformation 
cannot diverge at the horizon r = r+ nor at the asymptotic boundary r — t- oo. 

A TT gauge perturbation generated by ^ must satisfy the conditions V^^^ = 
and = 0. If we introduce the antisymmetric tensor F^i, = V[^^,y], for a Ricci flat 
background, these conditions reduce to V^^'^ = and V^F'^'^ = which read simply 



^.51) 



Using ^y—g = r'^'^^^\/g and Eq. (8.10), the first of the equations above requires that 



Xr2N. 



^.52) 



where the background function g{r) is defined in (8.1) and A is the CP eigenvalue 



(8.11). The second family of equations in (8.51) further demands that (C^)' (r) = ^i(r). 
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Introducing the new variable 



i.53) 



the solution of (8.51) must then solve 



x' {f) = V{r)x{r) , with V{r) 



N^ 



> 0. 



?.54) 



To study the regularity of the associated gauge transformation, we need the asymp- 



totic behavior of x(^) at the horizon and at infinity. The solution to Eq. (8.54) can be 



obtained in these regions by considering the dominant contributions of V{r) or by doing 



a Probenius analysis. Using (8.11), we find that 



X(?')L ~ ao(r - r+) or x{r)\ 



ao : 



X{r)\ 



bor 



i.55) 
i.56) 



Recall that 6h 



-C^g^u- We have to discard the second possibility in (8.55) because it 



would generate a dependence fu 



r+) ^, not compatible with the TT boundary 



conditions (8.46). On the other hand, we have to discard the solution with the positive 



sign in (8.56) because it generates a perturbation that grows faster than the unperturbed 



background metric as r — )• 0. The appropriate boundary conditions for (8.54) are thus 



X(r)|^ ~ ao(r - r+) , xWL ^ h,A-^^'^+^) . 



?.57) 



We can now complete our proof. Notice that 



0< / xHr? = XWX'W 

J r-i- 

fOO 



X{r)x"{r) 



r+ Jr+ 



F(r)x(r)" <0, 



^.58) 
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where we used (8.56) and (8.54). But these relations can be satisfied only for 



X(r) = 



^i(r) = 0, and ^+(r) = 0, 



^.59) 



which in addition to (8.50) proves that our regular zero- mode perturbations in the TT 



gauge cannot be pure gauge modes. 



8.B.3 Temperature and angular velocity preserved 



In Section 6.2, we discussed the connection between the classical instability of the black 
hole and its thermodynamics. This connection was built on the claim that the stationary 
and axisymmetric modes that we study preserve the temperature and the angular velocities 
of the background solution. Here we will prove that this is indeed the case. 

We will first compute the angular velocity and the temperature of the unperturbed 
background solution using standard Euclidean methods. Our strategy is then to check that 
our TT metric perturbation h^ydx^dx^ is a regular symmetric 2-tensor, when expressed 
in coordinates where the background metric is regular, which confirms that they preserve 
the angular velocity and temperature. 

We start with the computation of the background angular velocity and tempera- 
ture. This is done performing the standard three steps in the background solution: i) a 
coordinate transformation to coordinates (t, tp) that corotate with the black horizon, ii) a 
Wick rotation of the time coordinate so that we work with the Euclidean solution, and iii) 
a change to a new radial coordinate that zooms the geometry in the near-horizon region. 
That is, we perform the coordinate transformations 

A'(r4.1 

V' = V - ^Ht , t = -iT, r = r+ + — , (8.60) 



with il^H being the angular velocity of the background solution (8.2). A final coordinate 
transformation, 

f = 2^THr, with Th = ^^^, (8.61) 
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sets the period of r to be the inverse of the horizon temperature and avoids a conical sin- 
gularity at the horizon. The Euclidean sector of the near horizon region of the background 



solution (8.1) then reads 

dsl ~ p2 df 2 + dp'^ + h{r+f [d^ + Aadx^-f + rlgabdx^dx^ , (8.62) 

which is a manifestly regular geometry. Indeed the polar coordinate singularity can be 
removed by a coordinate transformation into cartesian coordinates, f = Arctan(y/3;) and 
p = a/x^ + y^. This concludes our computation of the angular velocity and temperature 
of the background black hole. 

To study the regularity of the perturbations, start by introducing the manifestly 
regular 1-forms, 

E"^ = p^df = xdy — ydx , E^ = p dp = x dx + y dy . (8.63) 
Consider now our TT metric perturbation h^ydx^dx^ . After using the boundary conditions 



(8.46), which satisfy the TT gauge conditions, we get 



h^y dx^" dx" ~ ^^^/oiY (p2 df^ + dp^) + /22Y (^dip + Aadx""^ ^ 



A^^' + W'-) dx^ + 2r+E^ (X+ + X-) dx^ (8-64) 



+ 2 (dV + Aadx^^ (Z+ + Z-) dx" 
1 



+ r'i 



dx^dx^ . 



Notice that the first term on the right-hand side ensures that there is no conical singularity 
if f has the same periodicity as the background metric. Furthermore, the remaining depen- 
dence on f and p is given by the manifestly regular 1-forms E"^ and EP. Hence, h^ydx^'dx'^ 
is a regular 2-tensor in the cartesian coordinates {x , y , ip , x"") where the background metric 
is regular, which confirms that the perturbations indeed preserve the angular velocity and 
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temperature. 



8.B.4 Perturbations of the asymptotic charges 

The stationary zero-modes can potentially change the mass and angular momenta of the 
geometry. In this subsection, we determine which of our perturbations with u = 0, m = 
and k = can change these conserved charges. 

The change on the conserved charges associated to a Killing generator ^ introduced 
by a perturbation /i^i/ can be defined via a surface boundary integral as (for TT pertur- 



The conserved charges of interest are the energy, for = —d/dt, and the angular momenta 
associated with the [{D — l)/2j U{1) Killing vectors ^ = d/d^i. The corresponding 
changes are denoted by £ and Ji, respectively. The surface integral is over a constant 
time hypersurface at asymptotic infinity, dT,. 

To compute the charges of our perturbations, we need the asymptotic behaviour of 
our solutions. This can be obtained from a Frobenius analysis of the Lichnerowicz equa- 
tions at r — 7- oo. We are interested in boundary conditions that preserve the asymptotic 
flatness of the spacetime, i.e. that decay (strictly) faster than the background geometry. 
We find that the behaviour of the regular perturbations is such that they decay at infinity 
according to 



bations) [l82l[l83 




(8.65) 





(8.66) 
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We then find the generic behaviour of the changes in the conserved charges 



S = Eor-^^l + 0{r-')) , 
J^,=Air^-^^ + Bir-^^ + 0{ 



y,-l-2K 



) 



(8.67) 



where Eo,Ai, Bi are functions of k and r+ and, in particular, Ai = Bi = for k = 0. 

We thus conclude that the k = modes are the only ones that decay sufficiently 
slowly to change the mass of the geometry. Moreover these modes cannot change the 

angular momenta. The only modes that change the angular momenta are those with 
K = 1. No modes with k > 2 can change the mass or the angular momenta. We have done 
these computations explicitly for the D = 7 case. 



Part IV 

Conclusion 
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Chapter 9 



Conclusion and outlook 



The original work described in this thesis can be divided into two parts. In the first, we 
studied the black hole partition function of Euclidean quantum gravity. We extended the 
study of negative modes, which represent pathologies in the one-loop quantum corrections, 
to black holes which are charged flGl] (Chapter [s]) or rotating |159| (Chapter |4| and 



also 68 , 72] , Chapters [7||8|) . We found that local thermodynamic instabilities are always 
signalled (individually) by the existence of a negative mode. The results strengthen the 
claim that the gravitational partition function indeed describes semiclassical quantum 
gravity at low energies, even beyond the leading order instanton approximation. In the 
charged trick based on a Kaluza-Klein reduction allowed for the decoupling of 

the unphysical divergent sector, the analogue of the conformal sector for pure-gravity. 
However, it would be convenient to have a more general procedure to deal with gravity- 
matter instantons, applicable to the charged rotating case. 



In the second part (Refs. [68,72 , Chapters [6||8]) , we explored the connection be- 
tween classical stability and local thermodynamic stability of black holes, numerically 
analysing perturbations of Myers-Perry solutions with a single spin and with equal spins 
(cohomogeneity-1 in odd D). The connection is that negative modes can have implica- 
tions for the classical stability problem. We started by refining the Gubser-Mitra conjec- 
ture, showing that one Gregory-Laflamme-type instability of uniform black branes appears 
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for each individual local thermodynamic instability, and in fact for each negative mode 
whether it is related to a thermodynamic instability or not. The onsets of these insta- 
bilities should be associated with bifurcations to new non-uniform black brane families. 
Moreover, we showed that all asymptotically vacuum black holes possess a local thermo- 
dynamic instability, which implies that the associated uniform black branes are classically 
unstable (68] . 

The main achievement in this thesis was to show that rapidly-rotating Myers-Perry 
black holes can be classically unstable, as first conjectured for the D > 6 singly-spinning 
sector in Ref. |71j. We presented in Ref. J72j (Chapter [T]) the first evidence for this in- 
stability by showing that, as high rotations are considered, additional negative modes of 
the partition function arise which are not associated with the standard local thermody- 
namic instabilities. The zero-modes marking the appearance of each new negative mode 
are instead the thresholds of classical instabilities of the black hole, and not just of the 
black branes. In Ref. [68] (Chapter [8|, we went further and, analysing cohomogeneity- 
1 MP black holes, verified explicitly the existence of such an instability in D = 9, by 
determining its timescale at the linear level. 

The non-thermodynamic zero-modes also mark the bifurcation of new families of 
stationary black holes from the MP family. In the singly-spinning MP case [72i|, there 
should be infinite new families of this type. The first may interpolate between a MP black 
hole and a black ring through a horizon topology transition, and the second may interpolate 



between a MP black hole and a black Saturn, etc., as conjectured in Refs. 62,71 . In the 
equal spins case [68] , these zero-modes indicate the existence of higher-dimensional black 
holes in D = 9 (and, we believe, higher odd D) with a single rotational symmetry. This 
is the first known example with spherical horizon topology (see Ref. [65] for the very first 
example, the helical black rings). Furthermore, we argue that this family of black holes 
which generically has a single rotational symmetry is determined by 70 parameters, while 
the MP solution has only 5 parameters. These results are a stark measure of the challenge 
of classifying higher-dimensional black holes. 



175 



The negative modes found obey a certain harmonic structure, with the thermody- 
namic negative modes corresponding to the lowest harmonics (s-wave and p-wave) There- 
fore, we conjectured that classical instabilities for perturbations preserving the rotational 
symmetry, corresponding to higher harmonics (d-wave, etc.), can only be excited for rota- 
tions higher than the first thermodynamic zero-mode (p-wave). This is the ultraspinning 



conjecture of Ref. 72 , consistent with the posterior results of Ref. |68|. The preservation 



of the rotational symmetry is required in our argument in order to make the connection 



with regular Euclideanised negative modes. It was recently found in Refs. |112, 113 that 
MP black holes can also be unstable after a critical value of the rotation for perturbations 
breaking this symmetry. This occurs even m D = 5, where there is no instability in the 
sector that we consider in this thesis. 

It would be important to focus also on different higher-dimensional solutions, such 
as black rings. What is perhaps the most important question of higher-dimensional black 
hole physics, along with the closely related classification problem, remains open: is there 
uniqueness of stable solutions? 



A major part of the work reported in this thesis, from Refs. 68,72, 159 , is based 
on a numerical analysis of coupled linear second order differential equations, ODEs or 
PDEs. The spectral method used to solve these equations is briefly described in the 
Appendix A. To the best of our knowledge, this powerful method was used in the context 
of general relativity for the first time in Ref. [159j . It would be interesting to also apply 
spectral methods to non-linear black hole problems, such as the construction of the new 
solutions whose existence we conjectured here, or of solutions relevant for applications of 
the holographic correspondence. Closely connected is the recent progress - likely to make 
an impact in the future - in extending to high energy physics problems the well-developed 



numerical techniques available for general relativity in four dimensions, e.g. 184 ,185 



Let us conclude on a different note. The plethora of higher-dimensional black holes 
and its seemingly impossible classification somewhat remind us of the situation in nu- 

^Notice that, in general, there may exist different harmonic structures corresponding to different per- 
turbation subsectors, e.g. there may exist unrelated p-waves. 
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clear/particle physics before the quark model was proposed. The challenge presented by 
higher dimensions contrasts with the simplicity of the four-dimensional case, where the 
Kerr black hole is the unique asymptotically flat vacTiTim sohition. In the same way, the 
myriad of particles detected in collision experiments contrasted with the previous situa- 
tion where only a handful of more familiar particles were known. The quark model put 
some order in that particle zoo. Our hope is that there is also an underlying structure 
connecting the several phases of black holes. This thesis is a small step towards unveiling 
that structure. Or one may ask the same question for the black ring that Rabi asked for 
the muon: "Who ordered that?" 



Appendix A 



Spectral numerical method 



In this Appendix, we briefly describe the spectral numerical method that we employed to 
solve the linear ordinary differential equations (ODEs) and the coupled partial differential 
equations (PDEs) in this thesis, namely in Chapters [4| [7] and [sj See Ref. [I60j for more 
details. 

Spectral methods can solve a system of coupled ODEs or PDEs to high accuracy 
on a finite domain, as long as the system allows for analytic solutions. To our knowledge, 



the first application of these methods to general relativity was given in 1156 , followed 



by 68,72 . We will start by considering the case of ODEs, and in the end we will discuss 



how to generalise the procedure for PDEs. 

The goal is to approximate a given function, defined on a finite domain, as a finite 
sum of algebraic polynomials p{y) = Y^'iLo'^iUi- Performing the polynomial interpolation 
in an equidistant grid with J\f + I points turns out to be catastrophic in many cases 
due to the oscillation of high-degree polynomials, the so-called Runge phenomenon. In 
general, this approximation does not converge as A/" — )• oo, and may get worse at a rate 
as large as 2-'^. The correct approach is to perform the interpolation in a non- uniform 
grid, distributed more densely near the edges of the interpolation interval. We use the 
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Chebyshev grid, whose points are the extrema of Chebyshev polynomials, 

%■ = — + — cos (^^ j , j E {0, 1, . . . , AA} , (A.l) 

for yj G [a,b]. Not only does this grid avoid the Runge phenomenon, by clustering points 
near the boundary, but it has another well-known advantage over uniform grids: it typically 
leads to an exponential accuracy of the approximantion as — )• oo. However, because 
we are approximating a function as a sum of polynomials, we must restrict to analytic 
functions. The reason for the exponential accuracy is that such functions have rapidly 
decaying Fourier transforms. Since spectral methods, as the name indicates, act in Fourier 
space in a certain sense, the common difficulties of having singular points in the equations 
can be avoided, as long as analytic solutions exist. 

The procedure to solve differential equations is in the same spirit as standard quan- 
tum mechanics. Consider an eigenvalue system of n coupled linear ODEs with variable 
coefficients, 

n n 

Y,Ha.pqf =XY.^o.pqf\ aG{0,l,...,n}, (A.2) 

/3=1 13=1 

where each i/^^ is a second order operator in y, each T^p is a scalar function and {A, q^jP} 
are the eigenvalues and eigenfunctions that we want to determine. We now perform our 
approximation and discretise the [a, 6] interval according to the grid (A.l). Each q^g'' is 



then approximated by a vector, whose entries are the values at yj of the eigenfunctions 
we want to determine. Following this procedure, one represents derivatives with respect to 
y by matrices, Dj^, that act on the vectors qj^, mixing adjacent points (see p. 53 of 
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for an explicit construction of such matrices). After this approach is complete, each 
and Tap are transformed into square matrices Hap and Tap, respectively, of dimension 
(AA + 1) X (AA + 1). 

We conveniently chose to work with Dirichlet boundary conditions by considering 



the rescaled functions (4.16), (7.8) and (8.30). The ellipticity of the perturbation operator 



is consistent with a boundary value problem. Dirichlet boundary conditions are imposed 
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by setting the first and last elements of each q"^ to zero, and by eliminating the first 
and last columns and rows of each matrix and T^p 
following system of linear algebraic equations 
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. We are then left with the 



P=l 13=1 



(A.3) 



where i^^^ and T^p are obtained from H^ip and T^p by deleting their first and last columns 



and rows, and thus are square {N — 1) x (A/" — 1) matrices. The system of equations (A.3) 
can be written as 



H 



11 



H 



In 



Hnl ■ ■ ■ H„ 



Hn 



Til ■ ■ ■ Tin 



Tnl • ■ ■ Tn 



9i 



(A.4) 



which is just a standard generalised eigenvalue problem of dimension n(J\f — 1). 

The generalisation to the cases where we have a system of PDEs is straightforward. 
Say we have two coordinates y and x. Then we can work with a two-dimensional grid 



where each coordinate is discretised according to (A.l). In our works, we have used the 



same number of points J\f for both y and x, although this is not a requirement. The vectors 
qj^^ which approximate the eigenfunctions q^^^ will have {J\f + 1)^ components instead of 
M + 1, one component for each grid point, and the matrices Hai^ and T^^ will accordingly 
become (AA + 1)^ x (A/" + 1)^ matrices, increasing the computational demands. 
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